2.3 初等模型-.pptVIP

数值积分及其在数模中的应用 %辛普森公式 z1=quad(1./(1-sin(x)),0,pi/4) %与精确值比较 dz1=z1-sqrt(2) %Gauss-Lobatto公式 z2=quadl(1./(1-sin(x)),0,pi/4) %与精确值比较 dz2=z2-sqrt(2) %梯形公式 x=0:pi/400:pi/4; y=1./(1-sin(x)); z3=trapz(y)*pi/400 z4=trapz(x,y) %与精确值比较 dz3=z3-sqrt(2) dz4=z4-sqrt(2) function y=weixing(t) a=7782.5;b=7721.5; y=sqrt(a^2*sin(t).^2+b^2*cos(t).^2); t=0:pi/10:pi/2; y1=weixing(t); 11=4*trapz(t,y1) 12=4*quad(‘weixing’,o,pi/2,1e-6) 用MATLAB 作数值积分 复化辛普森公式 quad(‘fun’,a,b,tol) 用辛普森(2阶)公式计算以 fun.m 命名的函数在 (a, b) 上的积分 tol为绝对误差，缺省时为10-6 quad8(‘fun’,a,b,tol) 用辛普森(8阶)公式计算 Gauss-Lobatto公式