wusemm 发表于 2010-6-13 19:15

双重积分问题

请教一个双重积分问题, 我的代码为
function =integ(x,y,gama)
gama=linspace(0,1,100); m=1; n=1; a=1; b=1;
k=((m*pi/a)^2+(n*pi/b)^2).*gama;
alpha=k.*a.*sin(x).*cos(y); beta=k.*b.*sin(x).*sin(y);
outnum=cos(alpha./2).*cos(beta./2);
outdem=((alpha/m/pi).^2-1).*((beta/n/pi).^2-1);
out=(outnum./outdem).^2.*sin(x);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
result=dblquad(@integ,0,pi/2,0,pi/2);

为什么就不能通过运行呢,gama要如何修改呢,请教

[ 本帖最后由 ChaChing 于 2010-6-13 21:24 编辑 ]

ChaChing 发表于 2010-6-20 00:37

参考下
关于带参数的积分问题 by xjzuo
http://forum.vibunion.com/forum/viewthread.php?tid=42369&extra=page%3D1%26amp%3Bfilter%3Ddigest
页: [1]
查看完整版本: 双重积分问题