smpsns_fxy.m

matlab源代码,适用于开发研究,带来很好的学习效果.

function INTf=smpsns_fxy(f,x,c,d,N)
%函数f(x,y)的一维数值积分数值;
%对应的积分区域是 Ry={c<=y<=d}
if nargin<5
    N=100;
end
if abs(d-c)<eps|N<=0
    INTf=0;
    return;
end
if mod(N,2)~=0
    N=N+1;
end
h=(d-c)/N;
y=c+[0:N]*h;
fxy=feval(f,x,y);
fxy(find(fxy==inf))=realmax;
fxy(find(fxy==-inf))=-realmax;
kodd=2:2:N;
keven=3:2:N-1;
INTf=h/3*(fxy(1)+fxy(N+1)+4*sum(fxy(kodd))+2*sum(fxy(keven)));