大虾，看一程序！feval()如何用？

jjkboy

<>function I=trapez_n(f_name,a,b,n)
h=(b-a)/n;
x=a+(0:n)*h;f=feval(f_name,x);
I=trapez_v(f,h)</P>
<>trapez_v.m</P>
<>function I=trapez_v(f,h)
I=h*(sum(f)-(f(1)+f(length(f)))/2);</P>

itmwk

<>把你的问题帖出来</P>