a=input('a=');
b=input('b=');
x0=input('x0');
y=input('y');
n=length(x0);
h=(b-a)/n;
x(1)=y(1)+y(n+1);
x(2)=0;
x(3)=0;
for i=2:n
xx=mod(i,2);
if (xx==1)
x(2)=x(2)+y(i);
else
x(3)=x(3)+y(i);
end
end
xi=h*(x(1)+2*x(3)+4*x(2))/3;
disp(xi)
这个写的自我感觉不好,改天重写一下