#### Previous topic

scipy.interpolate.PPoly.construct_fast

#### Next topic

scipy.interpolate.BPoly.__call__

# scipy.interpolate.BPoly

class `scipy.interpolate.``BPoly`(c, x, extrapolate=None, axis=0)

  Piecewise polynomial in terms of coefficients and breakpoints.

  The polynomial between `x[i]` and `x[i + 1]` is written in the Bernstein polynomial basis:

      S = sum(c[a, i] * b(a, k; x) for a in range(k+1)),

  where `k` is the degree of the polynomial, and:

      b(a, k; x) = binom(k, a) * t**a * (1 - t)**(k - a),

  with `t = (x - x[i]) / (x[i+1] - x[i])` and `binom` is the binomial coefficient.

  Parameters:

  c : ndarray, shape (k, m, …)

  > Polynomial coefficients, order k and m intervals

  x : ndarray, shape (m+1,)

  > Polynomial breakpoints. Must be sorted in either increasing or decreasing order.

  extrapolate : bool, optional

  > If bool, determines whether to extrapolate to out-of-bounds points based on first and last intervals, or to return NaNs. If ‘periodic’, periodic extrapolation is used. Default is True.

  axis : int, optional

  > Interpolation axis. Default is zero.

  See also

  `PPoly`
    piecewise polynomials in the power basis

  Notes

  Properties of Bernstein polynomials are well documented in the literature. Here’s a non-exhaustive list:

  [R52] | http://en.wikipedia.org/wiki/Bernstein_polynomial

  [R53] | Kenneth I. Joy, Bernstein polynomials, http://www.idav.ucdavis.edu/education/CAGDNotes/Bernstein-Polynomials.pdf

  [R54] | E. H. Doha, A. H. Bhrawy, and M. A. Saker, Boundary Value Problems, vol 2011, article ID 829546, DOI:10.1155/2011/829543.

  Examples

      >>> from scipy.interpolate import BPoly
      >>> x = [0, 1]
      >>> c = [[1], [2], [3]]
      >>> bp = BPoly(c, x)

  This creates a 2nd order polynomial

  \[\begin{split}B(x) = 1 \times b_{0, 2}(x) + 2 \times b_{1, 2}(x) + 3 \times b_{2, 2}(x) \\ = 1 \times (1-x)^2 + 2 \times 2 x (1 - x) + 3 \times x^2\end{split}\]

  Attributes

  x | (ndarray) Breakpoints.
  c | (ndarray) Coefficients of the polynomials. They are reshaped to a 3-dimensional array with the last dimension representing the trailing dimensions of the original coefficient array.
  axis | (int) Interpolation axis.

  Methods

  `__call__`(x[, nu, extrapolate]) | Evaluate the piecewise polynomial or its derivative.
  `extend`(c, x[, right]) | Add additional breakpoints and coefficients to the polynomial.
  `derivative`([nu]) | Construct a new piecewise polynomial representing the derivative.
  `antiderivative`([nu]) | Construct a new piecewise polynomial representing the antiderivative.
  `integrate`(a, b[, extrapolate]) | Compute a definite integral over a piecewise polynomial.
  `construct_fast`(c, x[, extrapolate, axis]) | Construct the piecewise polynomial without making checks.
  `from_power_basis`(pp[, extrapolate]) | Construct a piecewise polynomial in Bernstein basis from a power basis polynomial.
  `from_derivatives`(xi, yi[, orders, extrapolate]) | Construct a piecewise polynomial in the Bernstein basis, compatible with the specified values and derivatives at breakpoints.