In [1]:

    

%matplotlib inline
import matplotlib.pyplot as plt
import seaborn as sns
import numpy as np


    

In [2]:

    

func=lambda x:x**4-2*x+1


    

In [4]:

    

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())
print(I)


    

    




4.50656

In [6]:

    

import scipy.integrate as integrate


    

In [8]:

    

fun=lambda x:np.exp(-x)*np.sin(x)
result,error=integrate.quad(fun,0,2*np.pi)
print(result,error)


    

    




0.49906627863414593 6.023731631928322e-15

In [9]:

    

I=integrate.quad(fun,0,np.infty)
print(I)


    

    




(0.5000000000000002, 1.4875911973447031e-08)

In [10]:

    

print(abs(I[0]-0.5))


    

    




2.220446049250313e-16

In [ ]: