计算二重无穷积分
这是一个计算二重无穷积分的程序,
clc;clear;
syms u v tao;
x=4; y=9; a=2;
t0=0; t=3;
p=-exp(-tao);
A=eval((1/(4*a*pi*(t-t0)))*exp(int(p,tao,t0,t)));
g1=int(exp(-((x-u).^2+(y-v).^2)/(4*a*(t-t0)))*exp(-abs(u+v)),u,-inf,inf);
g2=vpa(int(g1,-inf,inf));
s=eval(A*g2)
现在不知道问题在哪里?有谁能帮忙吗?谢谢!