end
if nargout==2 out1=djf; out2=n; end case 'mjx'
x=linspace(a,b,n+1)'; %n等分点记入数组中 h=(b-a)/n;
djf=0; %累加器清零
for i=1:1:n
djf=djf+h*f((x(i)+x(i+1))/2); end
if nargout==0 out1=djf; end
if nargout==1 out1=djf; end
if nargout==2 out1=djf; out2=n; end case 'tx'
x=linspace(a,b,n+1)'; %n等分点记入数组中
h=(b-a)/n;
djf=0; %累加器清零
for i=1:1:n
djf=djf+(f(x(i))+f(x(i+1)))*h/2; end
if nargout==0
17
out1=djf; end
if nargout==1 out1=djf; end
if nargout==2 out1=djf; out2=n; end case 'sps'
x=linspace(a,b,n+1)'; %n等分点记入数组中
h=(b-a)/n;
djf=0; %累加器清零
for i=1:1:n
djf=djf+(h/6)*(f(x(i))+4*f((x(i)+x(i+1))/2)+f(x(i+1))); end
if nargout==0 out1=djf; end
if nargout==1 out1=djf; end
if nargout==2 out1=djf; out2=n; end case 'nc'
x=linspace(a,b,n+1)'; %n等分点记入数组中
h=(b-a)/n;
18
djf=0; %累加器清零
for i=1:1:n
djf=djf+(h/90)*(7*f(x(i))+32*f(x(i)+h/4)+12*f(x(i)+2*h/4)+32*f(x(i)+3*h/4)+7*f(x(i)+h));
end
if nargout==0 out1=djf; end
if nargout==1 out1=djf; end
if nargout==2 out1=djf; out2=n; end case 'zsy'
n=1;s1=(f(a)+4*f((a+b)/2)+f(b))*(b-a)/n/6; n=2; s2=s1+1;
while abs(s2-s1)>wc s1=s2; n=n*2; s=0; h=(b-a)/n; for i=1:n
s=s+4*f(a+(i-1/2)*(h)); end; for i=1:n-1 s=s+2*f(a+i*h);
19
end;
s=(f(a)+f(b)+s)*h/6; s2=s; end; if nargout==0 out1=s2; end
if nargout==1 out1=s2; end
if nargout==2 out1=s2; out2=n; end end
20
致 谢
这篇论文是在段敏老师细心的指导下完成的。导师渊博的专业知识,严谨的治学态度,精益求精的工作作风,诲人不倦的高尚师德,严以律己、宽以待人的崇高风范,朴实无华、平易近人的人格魅力对我影响深远。不仅使我树了远大的学术目标、掌握了基本的研究方法,还使我明白了许多待人接物与为人处世的道理。本论文从选题到完成,每一步都是在导师的指导下完成的,倾注了导师的心血。至此,我要向我的导师表示衷心的感谢!四年学业即将完成,人生学习的道路才刚刚开始,以后的道路上我将努力的学习,好好地工作,给老师争光,给学校争气!
黄 凯
2015年5月20日于合肥学院
21