fsolve为何不能得到结果？

suntht

用fsolve命令求函数 BesselJ(0,x)*BesselI(1,x)+BesselI(0,x)*BesselJ(1,x) 的零点，在x比较小的时候能得出结果，但x>58之后不能计算了，如图所示。
这是什么原因？

suntht

啊哈，我自己回答吧。是因为数据太大了，到1e24了，把公式做个变化就可以了。比如每一项都除以1e10，i写成 “BesselJ(0,x)*BesselI(1,x)/1e10+BesselI(0,x)*BesselJ(1,x)/1e10” 就能算了。