dblsim.m

来自「完全的二维数值积分程序 给出自变量和函数值的矩阵」· M 代码 · 共 37 行

function I=dblsim(x,y,f)

% figure;surf(x,y,f);

[n m]=size(f);
m=m-1;n=n-1;
U=ones(1,m+1);
r=1:2:m+1;
U(r)=U(r)*2/3;
r=2:2:m;
U(r)=U(r)*4/3;
U(1)=1/3;U(m+1)=1/3;

v=ones(1,n+1);
r=1:2:n+1;
v(r)=v(r)*2/3;
r=2:2:n;
v(r)=v(r)*4/3;
v(1)=1/3;v(n+1)=1/3;

t1=ones(1,m+1);
t2=ones(1,n+1);
h=x(1,2)-x(1,1);
k=y(2,1)-y(1,1);

I=h*k*t2*((v'*U).*f)*t1';