Matlab数值积分程序集合.docVIP

Matlab数值积分程序集合[图书馆+网络收集] 近来学习数值积分,手头积累了不少程序,也拿来和各位朋友分享一下。。。主要是来自数值积分教材和网络，基本的原理也就不打算多说了，随便搜索一下就可以得到，那就开始上代码了，呵呵，非原创，但是全部验证过，有疑问可以给我e-mail： 1 梯形数值积分的MATLAB主程序 function T=rctrap(fun,a,b,m)??????? %fun 函数，a 积分上限 b积分下限 m 递归次数 n=1;h=b-a; T=zeros(1,m+1); x=a; T(1)=h*(feval(fun,a)+feval(fun,b))/2; for i=1:m ?????????? h=h/2; n=2*n; s=0; ????????? for k=1:n/2 ?????????? x=a+h*(2*k-1); s=s+feval(fun,x); end T(i+1)=T(i)/2+h*s; end T=T(1:m);e.g 运行后屏幕显示 精确值Fs，用rctrap计算的递归值T和T与精确值Fs的绝对误差wT ) exp((-x^.2./2)./(sqrt(2*pi)))T=rctrap(fun,0,pi/2,14), syms tfi=int(exp((-t^2)/2)/(sqrt(2*pi)),t,0, pi/2);Fs