SciPy

numpy.polyder

numpy.polyder(p, m=1)[source]

Return the derivative of the specified order of a polynomial.

Parameters:

p : poly1d or sequence

Polynomial to differentiate. A sequence is interpreted as polynomial coefficients, see poly1d.

m : int, optional

Order of differentiation (default: 1)

Returns:

der : poly1d

A new polynomial representing the derivative.

See also

polyint
Anti-derivative of a polynomial.
poly1d
Class for one-dimensional polynomials.

Examples

The derivative of the polynomial

System Message: WARNING/2 (x^3 + x^2 + x^1 + 1)

latex exited with error [stdout] This is pdfTeX, Version 3.14159265-2.6-1.40.16 (TeX Live 2015/TeX Live for SUSE Linux) (preloaded format=latex) restricted \write18 enabled. entering extended mode (./math.tex LaTeX2e <2015/01/01> patch level 2 Babel <3.9m> and hyphenation patterns for 79 languages loaded. (/usr/share/texmf/tex/latex/base/article.cls Document Class: article 2014/09/29 v1.4h Standard LaTeX document class (/usr/share/texmf/tex/latex/base/size12.clo)) (/usr/share/texmf/tex/latex/base/inputenc.sty ! LaTeX Error: File `utf8x.def’ not found. Type X to quit or <RETURN> to proceed, or enter new name. (Default extension: def) Enter file name: ! Emergency stop. <read *> l.173 \endinput ^^M No pages of output. Transcript written on math.log.
is:

>>> p = np.poly1d([1,1,1,1])
>>> p2 = np.polyder(p)
>>> p2
poly1d([3, 2, 1])

which evaluates to:

>>> p2(2.)
17.0

We can verify this, approximating the derivative with (f(x + h) - f(x))/h:

>>> (p(2. + 0.001) - p(2.)) / 0.001
17.007000999997857

The fourth-order derivative of a 3rd-order polynomial is zero:

>>> np.polyder(p, 2)
poly1d([6, 2])
>>> np.polyder(p, 3)
poly1d([6])
>>> np.polyder(p, 4)
poly1d([ 0.])

Previous topic

numpy.polyfit

Next topic

numpy.polyint