[未答复] 三重积分积分不成功,怎么办

[复制链接]
造梦_lHwW0 发表于 7 天前
本帖最后由 造梦_lHwW0 于 2021-5-3 11:55 编辑

G=6.67*10^-11;
M=1.8982*10^27;
T=86400;
P=7.98*10^3;
h=4.5*10^6;
d=1*10^5;
r=2.8*10^8;
func1=@(x,y,z)G.*P.*(y+r)./((y+r).^2+x.^2+z.^2).^1.5;
I1=integral3(func1,-(r+d),(r+d),@(x)-sqrt((r+d).^2-x.^2),@(x)sqrt((r+d).^2-x.^2),-h/2,h/2)-integral3(func1,-(r),r,@(x)-sqrt((r).^2-x.^2),@(x)sqrt((r).^2-x.^2),-h/2,h/2)
////
求不出来I1显示如下
////
警告: 已达到函数计算的最大数目(10000)。结果未通过全局误差测试。
> 位置:integral2Calc>integral2t (第 129 行)
位置: integral2Calc (第 9 行)
位置: integral3/innerintegral (第 146 行)
位置: integralCalc/iterateScalarValued (第 314 行)
位置: integralCalc/vadapt (第 132 行)
位置: integralCalc (第 75 行)
位置: integral3 (第 121 行)
警告: 积分未成功。
////
想问下这个程序应该怎么改


您需要登录后才可以回帖 登录 | 注册

本版积分规则

相关帖子
热门教程
站长推荐
快速回复 返回顶部 返回列表