How to compute derivative using Numpy?
Last Updated :
21 Apr, 2021
In this article, we will learn how to compute derivatives using NumPy. Generally, NumPy does not provide any robust function to compute the derivatives of different polynomials. However, NumPy can compute the special cases of one-dimensional polynomials using the functions numpy.poly1d() and deriv().
Functions used:
- poly1d(): It helps to define a polynomial expression or a function.
- deriv(): Calculates and gives us the derivative expression
Approach:
- At first, we need to define a polynomial function using the numpy.poly1d() function.
- Then we need to derive the derivative expression using the derive() function.
- At last, we can give the required value to x to calculate the derivative numerically.
Below are some examples where we compute the derivative of some expressions using NumPy. Here we are taking the expression in variable ‘var’ and differentiating it with respect to ‘x’.
Example 1:
Python3
import numpy as np
var = np.poly1d([ 1 , 0 , 1 ])
print ( "Polynomial function, f(x):\n" , var)
derivative = var.deriv()
print ( "Derivative, f(x)'=" , derivative)
print ( "When x=5 f(x)'=" , derivative( 5 ))
|
Output:
Example 2:
Python3
import numpy as np
var = np.poly1d([ 4 , 9 , 5 , 1 , 6 ])
print ( "Polynomial function, f(x):\n" , var)
derivative = var.deriv()
print ( "Derivative, f(x)'=\n" , derivative)
print ( "When x=3 f(x)'=" , derivative( 3 ))
|
Output:
Example 3:
Python3
import numpy as np
var = np.poly1d([ 5 , 4 , 9 , 5 , 1 , 6 ])
print ( "Polynomial function:\n" , var)
derivative = var.deriv()
print ( "Derivative, f(x)'=\n" , derivative)
print ( "When x=2 f(x)'=" , derivative( 0.2 ))
|
Output:
To calculate double derivative we can simply use the deriv() function twice.
Example 4:
Python3
import numpy as np
var = np.poly1d([ 3 , 5 , 4 , 9 , 5 , 1 , 6 ])
print ( "Polynomial function:\n" , var)
derivative = var.deriv()
print ( "Derivative, f(x)'=\n" , derivative)
print ( "When x=1 f(x)'=" , derivative( 1 ))
derivative1 = derivative.deriv()
print ( "\n\nDerivative, f(x)''=\n" , derivative1)
print ( "When x=1 f(x)'=" , derivative1( 1 ))
|
Output:
Share your thoughts in the comments
Please Login to comment...