利用MATLAB软件编写辛普森求积公式程序

a=input('a=?');

b=input('b=?');

x0=(a+b)/2;

f=[0 0 0];

f(1)=a^2*cos(a);

f(2)=x0^2*cos(x0);

f(3)=b^2*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 =

-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