24小时热门版块排行榜    

查看: 471  |  回复: 0

我爱该研究

铁虫 (初入文坛)

[求助] 求助matlab程序(高斯积分),总是有问题

求助程序(高斯积分),总是有问题
高手们帮忙看看,哪里有问题,积分结果应该是收敛的,
syms x;  
a=0;b=30;p=0.3;v=0.33;r=0.2;H=1;
n=10;
G=@grule;           %调用子函数,计算高斯积分点
[XK,AK]=G(n);        
XK1=((b-a)/2)*XK+((a+b)/2);
f=p*x(1/(4*(v-1))-(1-(1+2*H*x+2*H^2*x^2)*exp(-2*H*x))/(4*(v-1)))*...
    besselj(1,p*x)*besselj(0,r*x);
I=((b-a)/2)*sum(AK.*subs(sym(f),symvar(f),XK1));
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

相关版块跳转 我要订阅楼主 我爱该研究 的主题更新
信息提示
请填处理意见