二重积分结果出来是N​aN,上下限相同结果​本应是0。

版本：matlab2018b；

命令：

R=1;Ba=0;a=1;x=a;

fun3=@(r,theta)(1-theta.*cot(theta))./sqrt(1-sin(theta).*(r./R).^2);

thetamin=@(r)asin(a./r);

Nx3=integral2(fun3,x,1,thetamin,pi/2);

运行后命令窗口显示：

警告: 非有限结果。积分未成功。可能具有奇异性。

> In integral2Calc>integral2t (line 121)

In integral2Calc (line 9)

In integral2 (line 106)

Nx3 =

NaN

因为x=1，其上下限相同结果本应是0，不知道为什么结果出现奇异性。

Answer 1

[hxtgkmy]

但是有数值精度误差呀

你用符号积分int试试