function g=myintegrate(f,N,a,b)
%integration using simpson method
n=N/2;
h=(b-a)/N;
x=a:h:b;
fv=f(x);
s=sum(fv(1:2:2*n-1))*4;
s=s+sum(fv(2:2:2*n))*2;
s=s+f(a)+f(b);
s=s*h/3;
g=s;
现在要做的工作就是如何把上面的f函数换成我们想要的随便什么函数。
定义函数句柄:
handle=@ functionname
之后就可以把handle作为参数传递给上面的函数(或者其他需要使用函数作为参数的函数)了。
例如在0~pi/2上对sin(x)进行积分,可以如下输入:
>>fh=@sin;
>>myintegrate(fh,100,0,pi)
ans =
函数句柄可以指定为任意的matlab函数——当然也包括我们自定义的函数,比如我们在myfun.m文件中定义了如下函数:
function f=myfun(x)
f=x;
那么,
>>fh2=@myfun;
>>myintegrate(fh2,100,0,1)
ans =