function I=SimposnEq(a,b,n,y)
%等距离散节点积分之Simposn法
%[a,b]作n等分(有n+1个节点)
%y:n+1个节点对应的n+1维值向量
su=0;
for i=1:n-1
   if mod(i,2)==1
       su=su+4*y(i+1);
   else
       su=su+2*y(i+1);
   end
end
   I=(b-a)/(3*n)*(y(1)+y(n+1)+su);


function I=TrapeIneq(x,y)
%非等距离散节点积分之梯形法
%
%y:值向量
n=length(x);
I=0;
for i=1:n-1
   I=I+abs(y(i)+y(i+1))*abs(x(i+1)-x(i))/2;
end

非常感谢，我做做试试，现在处理数据非常需要，有不懂得地方，还会打扰的！

help trapz

还请帮忙，这两个函数我都用不了怎么回事哪？

还请帮忙，那两个函数我都用不了怎么回事啊？