import numpy as np
#def funcval(x):
# return np.exp**(1/X)
def funcvals(x1):
a,b=8755,6810
return (a**2*np.sin(x)**2+b**2*np.cos(x)**2)**(1/2)
#梯形公式
def trapezoid(a,b,n,f):
sum=0
sum=(f(a)+f(b))*(b-a)/2
return sum
#中矩形公式
def retangle(a,b,n,f):
sum=0
sum=f((a+b)/2)*(b-a)
return sum
#simpson
def simpson(a,b,n,f):
sum=0
sum=4*f((a+b)/2)+f(a)+f(b)
return sum*(b-a)/6
#ftrapezoid
def trapezoid(a,b,n,f):
h=(b-a)/(n-1)
x=np.zeros(n,1)
for i in range(n):
x[i]=a+i*h
sum=0
sum=f(x[0])+f(x[-1])
sum=sum+2*np.sum(f(x[1:-1]))
sum=sum*h/2
return sum
#fsimpson
def fsimpson(a,b,n,funcval1):
h=(b-a)/(n-1)
x=np.zeros(n,1)
for i in range(n):
x[i]=a+i*h
h=h*2
sum=f(x[0])+f(x[-1])
if x.shape[0]>2:
#griad
fg=f(x[2:-1:2])
fh=f(x[1:-1:2])
#half grid
sum=sum+2*np.sum(fg)+4*np.sum(fh)
sum=sum*h/6
return sum
def main():
a,b=1,2
n=20
tr=trapezoid(a,b,n,funcval)
sip=simpson(a,b,n,funcval)
ft=ftrapezoid(a,b,nfuncval)
fs=fsimpson(a,b,n,funcval)
a,b=0,np.pi/2
fs=fsimpon(a,b,n,funcval1)
print(fs*4)
if __name__=='__main__':
main