function I=simpson(fun,a,b,k)n=2.*k+1;h=(b-a)/(n-1);x=a:h:b;f=feval(fun,x);I=(h/3).*(f(1)+4.*sum(f(2:2:n-1))+2.*sum(f(3:2:n-2))+f(n));
你说的是simpson求积公式吗？quad这个函数就是这个方法求定积分
