Matlab7.0中二重积分的数值计算问题symsxy;symsnnxnny;nx=int(nnx);ny=int(nny);nx=0;ny=0;maple('with(orthopoly)');Hx=maple('H(nx,x)');maple('with(orthopoly)');Hy=maple('H(ny,y)');a=1.6627;Nx=(a/sqrt(pi)*2^nx*factorial(nx))^0.5;Ny=(a/sqrt

Question

Matlab7.0 中二重积分的数值计算问题
syms x y;
syms nnx nny;
nx=int(nnx);
ny=int(nny);
nx=0;
ny=0;
maple('with(orthopoly)');
Hx=maple('H(nx,x)');
maple('with(orthopoly)');
Hy=maple('H(ny,y)');
a=1.6627;
Nx=(a/sqrt(pi)*2^nx*factorial(nx))^0.5;
Ny=(a/sqrt(pi)*2^ny*factorial(ny))^0.5;
fx=Nx*exp(-0.5*a^2*x^2)*Hx;
fy=Ny*exp(-0.5*a^2*y^2)*Hy;
Fxy=fx*fy;
Goal1=diff(diff(Fxy,x),y);
Goal=Goal1*Fxy;
Final=@(x,y)Goal;
result= dblquad(Final,-inf,inf,-inf,inf)
为什么
rror in ==> Try1 at 20
result= dblquad(Final,-inf,inf,-inf,inf)

Answer 1

不能