python解非线性方程的问题

[zengxj100]

废话不多说，上代码
方法1：

from sympy import *

T=symbols('T')
eq=Eq(10**20*exp(-12000/T)/(0.01+10**14*exp(-12000/T))-18700*(T-300))
#eq=[10**20*exp(-12000/T)/(0.01+10**14*exp(-12000/T))-18700*(T-300)]#改成这样结果也是一样的
result=nonlinsolve(eq,T)
print(result)

程序运行，返回的是一个ConditionSet对象，说明不能解这个方程。

方法2：
import numpy as np
from scipy.optimize import fsolve

def f(T):
return 10**20*np.exp(-12000/T)/(0.01+10**14*np.exp(-12000/T))-18700*(T-300)

result=fsolve(f,[300])
print(result)

运行结果报错
minpack.error: Result from function call is not a proper array of floats.

方法3
import math
from scipy.optimize import fsolve

def f(T):
return 10**20*math.exp(-12000/T)/(0.01+10**14*math.exp(-12000/T))-18700*(T-300)
result=fsolve(f,[300])
print(result)

第3种方法能解出一个根 [303.3]，但是方程有3个根。如果把“result=fsolve(f,[300])”这句改成“result=fsolve(f,[320])“，可以得到第二个根 [323.7]；如果“result=fsolve(f,[300])”这句改成“result=fsolve(f,[360])“，可以得到第三个根 [349.4]。

如果把“result=fsolve(f,[300])”这句改成“result=fsolve(f,[300,320.360])“，则会报错
TypeError: only size-1 arrays can be converted to Python scalars

各位大佬，如何才能把三个根一次性解出来呢？