代码:
from scipy.optimize import fsolve
import math
m = 1.2764
def func(x):
return (4 * x * math.gamma(x) ** 2 / math.pi / math.gamma(x + 0.5) ** 2 - m)
v21 = fsolve(func, 0.01)
print(v21)
问题:解非线性方程,当m>=1.2765时候,可以正确解出x值,但是当m<=1.2764时候,程序报错 OverflowError: (34, ‘Result too large’)
这个程序我是按照matlab程序翻译成python的,在matlab里面,m<=1.2764时候,都可以解出答案,在python中就报错,请大神指导原因并解决方法,多谢!!