《利用MATLAB软件编写辛普森求积公式程序.docx》由会员分享，可在线阅读，更多相关《利用MATLAB软件编写辛普森求积公式程序.docx（2页珍藏版）》请在三一办公上搜索。
1、利用MATLAB软件编写辛普森求积公式程序a=input(a=?); b=input(b=?); x0=(a+b)/2; f=0 0 0; f(1)=a2*cos(a); f(2)=x02*cos(x0); f(3)=b2*cos(b); s=(b-a)/6)*(f(1)+4*f(2)+f(3) a=input(a=?); b=input(b=?); x=a:(b-a)/2:b f=(x.2).*cos(x) s=(b-a)/6)*(f(1)+4*f(2)+f(3) a=?1 b=?2 x = 1.0000 1.5000 2.0000 f = 0.5403 0.1592 -1.6646 s =
2、 -0.0813 function C,D=newpoly(X,Y) X=0.4 0.55 0.65 0.8 0.9 1.05; Y=0.41075 0.51875 0.69675 0.88811 1.02652 1.25382; %X=1 2 3; Y%=4 5 7; n=length(X); d=zeros(n,n); D(:,1)=Y for j=2:n for k=j:n D(k,j)=(D(k,j-1)-D(k-1,j-1)/(X(k)-X(k-j+1); end end C=D(n,n) for k=(n-1):-1:1 C=conv(C,poly(X(k); m=length(C) C(m)= C(m)+D(k,k); C end
