matlab四重积分

这是在用matlab做光学衍射模拟的时候衍生出的一个问题

clear all;

f = 'x^2+y^2+m^2+n^2'

M = 10;

N = 10;

A = zeros(M,N);

for m = 1:M

for n = 1:N

ftmp = f;

ftmp = subs(ftmp,'m',num2str(m));

ftmp = subs(ftmp,'n',num2str(n));

f1 = @(x,y) eval(ftmp);

A(m,n) = dblquad(f1,0,1,0,1)

end

end

或者

clear all;

sum = 0;

f2 = 'x^2+y^2+m^2+n^2'

for n = 1:10

f = f2;

f = subs(f,'n',num2str(n));

f1 = @(x,y,m) eval(f);

sum = sum + triplequad(f1,0,1,0,1,0,1)

end