积分与微-余显祥.docVIP

习题一 右矩形法 xitiyi.a for n=[20 25 40];%对n赋予一个向量 h=(3-1)./n;%计算步长 x=1:h:3;%取节点 k=2; f=0; while k=n%循环计算出f f=cos(x(k))*exp(sin(x(k)))+f; k=k+1; end s=f.*h%计算体积 end s = -1.1720 s = -1.1716 s = -1.1706 梯形法 xitiyib.m for n=[20 25 40] %对n赋予一个向量 h=(3-1)./n;%计算步长 x=1:h:3;%取节点 y=cos(x).*exp(sin(x)); s=h.*(y(1)+y(n+1))/2;%计算两个端点的数值 for k=2:n%循环计算出积分s s=h.*y(k)+s; end s end s = -1.1663 s = -1.1670 s = -1.1677 三点高斯积分法 xitiyic.m t=[-0.7745 0 0.7745];%三点高斯积分点 x=t+2;%进行变量替换 f=cos(x).*exp(sin(x)); s=0.5556*f(1)+0.8889*f(2)+0.5556*f(3)%用高斯积分点积分求数值解 s =