积分

#P102 上半部分
#4.1.2
#例子2 数值积分

#f(x)=9-x^2,从-3积分到3的结果为36

import numpy as np
from scipy.integrate import simps,romberg


a=-3.0;b=3.0
N=10


x=np.linspace(a,b,N)
y= 9-x*x
Yromb=lambda x:(9-x*x)
t=np.trapz(y,x)
s=simps(y,x)
r=romberg(Yromb,a,b)
aiv= (9*b-(b*b*b)/3.0)-(9*a-(a*a*a)/3.0)


print 'trapezoidal={0}({1:%}error)'.format(t,(t-aiv)/aiv)
print 'simpsons={0}({1:%}error)'.format(s,(s-aiv)/aiv)
print 'romberg={0}({1:%}error)'.format(r,(r-aiv)/aiv)
print 'actural value={0}'.format(aiv)