python - 如何在 Python 中求解 (x+1)e^x=c(constant)？

Question

我想用 Python 求解 (x+1)e^x=c 方程。

如下图所示，已使用 lambert w 函数手动成功求解该方程:

使用相同的步骤，我想以编程方式解决 (x+1)e^x。我按照上图中所示的步骤使用模块 SymPy 对其进行了编码，但没有成功。

有没有Python可以解决这类方程式的？

import numpy as np
from sympy import *
n = symbols('n')
sigmao=0.06866
sigmas=0.142038295
theta=38.9
rad=(np.pi/180)*38.9076
cos=np.cos(rad)
sec=1/np.cos(rad)
out = (0.06*0.7781598455*n*(1-exp(-2*0.42*sec*n))+exp(-2*0.42*n*sec)*sigmas)/sigmao
#Apply diff for the above expression. 
fin=diff(out, n)
print(solve(fin,n))

Answer 1

from scipy.optimize import fsolve
import numpy as np

const = 20
def func(x):
    return [(x[0]+1) * np.exp(x[0]) - const]

result = fsolve(func, [1])[0]
print('constant: ', const, ', solution: ', result)
#check
print('check: ', (result+1) * np.exp(result))


#Output[]:
    constant:  20.0 , solution:  1.9230907433218063
    check:  20.0

预览:https://onlinegdb.com/By8Z2Jwgw