In [3]:
%matplotlib inline
import matplotlib.pyplot as plt
import seaborn as sns
import numpy as np
In [7]:
func=lambda x:x**4-2*x+1
In [8]:
N=10
a=0.0
b=2.0
h=(b-a)/N
k=np.arange(1,N)
I=h*(0.5*func(a)+0.5*func(b)+func(a+k*h).sum())
In [9]:
print(I)
In [ ]: