复合梯形公式求积分c语言,用复合梯形公式和复合辛普森公式求函数积分..doc

用复合梯形公式和复合辛普森公式求函数积分.

附录一：

《数值分析》实验报告(模板)

学号 班级 姓名

【实验课题】

【实验目标】

【理论概述与算法描述】，将区间【a,b】划分为n等份，分点x(k)=a+kh,h=(b-a)/n,k=0,1,2,3, ……，在每个区间【x(k),x(k+1)】(k=0,1,2……n-1)上采用梯形公式，得

，

记

，则此公式Tn为复合梯形公式。

2．根据辛普森公式，将区间【a,b】分为n等分，在每个区间【x(k),x(k+1)】上采用辛普森公式，记x(k+1/2)=x(k)+k/2,则得到

，

记

，此公式为复合辛普森求积公式。【实验问题】，(精确解：I=1)

2 ，(取，)

【实验过程与结果】

【结果分析、讨论与结论】

h=b-a;

k=0;

T=((f(a)+f(b))*h)/2;

P=T+1;

while abs(P-T)>tol

P=T;

m=0; h=h/2;

for i=1:2^k

m=m+f(a+(2*i-1)*h);

end

T=0.5*P+m*h; k=k+1;

end

复合辛普森公式

function S=comsinp(a,b,tol)

h=b-a;

k=1;

S=((f(a)+f(b)+4*f((a+b)/2))*h)/6;

P=S;

while abs(P-S)>tol

P=S;

m=0;

n=0;

h=h/2;

for i=1:2^k-1

m=m+f(a+((2*i+1)/2)*h);

end

for j=1:2^k-1

n=n+f(a+j*h);

end

S=0.5*P+(h*(4*m+2*n))/6;

end

主程序

main

a=0;

b=1;

tol=10^-6;

T=comptra(a,b,tol)

s=comsinp(a,b,tol)

令y=f(x) 分别令f(x)=sinx.sin(x+eps)/(x+eps)

数值分析实验指导

第1页