在python3中使用fsolve

Question

我正在尝试使用Python3找到以下方程的解： 83.045269**(8)np.exp(-3.045269a) + 63.516372**(8) np.exp(-[a+1.042817]3.516372) - 42.8578**8 * np.exp(-2.8578a) - 4*3.3**(8) *exp(-[a+1.042817]*3.30) = 0 我做了以下事情：

>>> import numpy as np
>>> from scipy.optimize import fsolve
>>> def equation(a):
...     term1 = 8 * 3.045269**8 * np.exp(-3.045269)
...     term2 = 6 * 3.516372**8 * np.exp(-(a + 1.042817) * 3.516372)
...     term3 = -4 * 2.8578**8 * np.exp(-2.8578 * a)
...     term4 = -4 * 3.3**8 * np.exp(-(a + 1.042817) *3.3)
...     return term1 + term2 +term3 + term4
... initial_guess = 0.5
File "<stdin>", line 7
initial_guess = 0.5
^^^^^^^^^^^^^
SyntaxError: invalid syntax
 >>> initial guess = 0.5
File "<stdin>", line 1
initial guess = 0.5
        ^^^^^
 SyntaxError: invalid syntax
>>> 

以上有什么问题吗？

Answer 1

这不是关于 fsolve 的问题，而是关于正确使用 Python REPL 的问题。使用 Python REPL 时，您需要注意称为主提示的三个大于号和称为辅助提示的三个点之间的区别。

主提示 (>>>) 表明口译员已准备好接受输入。当您输入复合语句或续行时，会出现辅助提示 (...)。 在您的例子中，您定义了一个带有参数

equation

的函数

a

。按回车键后，会出现后续代码行的辅助提示，表示它们位于功能块中。在 return 语句之后，您需要退出功能块，这意味着您应该再次按 Enter 键，直到看到主提示。本质上，仅当您看到

initial_guess = 0.5

提示时才应输入

>>>

，即：

>>> import numpy as np
>>> from scipy.optimize import fsolve
>>> def equation(a):
...       term1 = 8 * 3.045269**8 *   np.exp(-3.045269)
...       term2 = 6 * 3.516372**8 * np.exp(-(a + 1.042817) * 3.516372)
...       term3 = -4 * 2.8578**8 * np.exp(-2.8578 * a)
...       term4 = -4 * 3.3**8 * np.exp(-(a + 1.042817) *3.3)
...       return term1 + term2 +term3 + term4
... 
>>> initial_guess = 0.5