n=100时梯形法运行结果
n=1000时梯形法运行结果
n=5000时梯形法运行结果
3、辛普生法： 程序编写如下：
real a,b,s integer n real sinpson
print*,'输入a,b和n的值' read*,a,b,n
s=sinpson(a,b,n) print 10,a,b,n print 20,s
10 format('a=',f5.2,3x,'b=',f5.2,3x,'n=',i4) 20 format('s=',f15.8) end
real function sinpson(a,b,n) implicit none real a,b,h,f2,f4,x integer i,n real func
h=(b-a)/(2.0*n) x=a+h f2=0
f4=func(x) do i=1,n-1 x=x+h
f2=f2+func(x) x=x+h
f4=f4+func(x) enddo
sinpson=(func(a)+func(b)+4.0*f4+2.0*f2)*h/3.0 end
real function func(x) real x
func=1+sin(x) end
n=10时辛普生法运行结果
n=100时辛普生法运行结果
n=1000时辛普生法运行结果
n=5000时辛普生法运行结果 三、利用Gauss-Jordan法求联立方程组：
{ x+4y+7z=12 2x+5y+8z=15 3x+6y+8z=17