查看: 1318  |  回复: 1

carter191

木虫 (小有名气)

[求助] 我想用quad2dggen函数计算一个二重积分

我想用quad2dggen函数求解二重积分总出错,请会的大神帮忙看看。
子函数:
function z=eff(x,y)
z=A^2/(besselk(m,W)*besselj(m,U))*besselk(m,W*sqrt((d+x)^2)/a2)*besselj(m,U*sqrt(x^2+y^2)/a1)*cos(m*atan(y/(d+x)))*cos(m*acos(x/sqrt(x^2+y^2)))%其中A,m,W,U,a1,a2,d都是已知量
主程序:
clear all
clc
clf
syms x y;
C=quad2dggen('eff',-sqrt(160^2-y^2),sqrt(160^2-y^2),-160,160,0.001,1);
运行的时候总出错,错误如附图。matlab初学,没用过quad2dggen函数请会的大神指点下,这个程序该怎么改。

我想用quad2dggen函数计算一个二重积分
新建位图图像.jpg
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

carter191

木虫 (小有名气)

A,m,W,U,a1,a2,d 赋值分别为 1,0,50,40,200,50,10
2楼2014-04-24 10:08:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 carter191 的主题更新
信息提示
请填处理意见