def _from_derivatives(xi, yi, x, order=None, der=<NUM_LIT:0>, extrapolate=False):
from scipy import interpolate<EOL>method = interpolate.BPoly.from_derivatives<EOL>m = method(xi, yi.reshape(-<NUM_LIT:1>, <NUM_LIT:1>),<EOL>orders=order, extrapolate=extrapolate)<EOL>return m(x)<EOL>
Convenience function for interpolate.BPoly.from_derivatives.

Construct a piecewise polynomial in the Bernstein basis, compatible
with the specified values and derivatives at breakpoints.

Parameters
----------
xi : array_like
    sorted 1D array of x-coordinates
yi : array_like or list of array-likes
    yi[i][j] is the j-th derivative known at xi[i]
order: None or int or array_like of ints. Default: None.
    Specifies the degree of local polynomials. If not None, some
    derivatives are ignored.
der : int or list
    How many derivatives to extract; None for all potentially nonzero
    derivatives (that is a number equal to the number of points), or a
    list of derivatives to extract. This numberincludes the function
    value as 0th derivative.
 extrapolate : bool, optional
    Whether to extrapolate to ouf-of-bounds points based on first and last
    intervals, or to return NaNs. Default: True.

See Also
--------
scipy.interpolate.BPoly.from_derivatives

Returns
-------
y : scalar or array_like
    The result, of length R or length M or M by R.