这是一个计算二重无穷积分的程序，
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)

现在不知道问题在哪里？有谁能帮忙吗？谢谢！