|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
请教一个双重积分问题, 我的代码为
function [out]=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 编辑 ] |
|