我想用python得到不同温度下的振动配分函数。我有两个一维数组:freq=[100,150,200]
(这实际上是从.dat文件中读取的)和Temp=np.arange(273.,973.,25)
。现在,对于一个特定的温度(1/1-exp(-h*nu/kb/T))
,通过频率qvib=product。我的for循环是:
for i in range(len(T)):
for j in range(len(freq)):
q[i][j]=1/(1-np.exp(-h*freq[j]/kb/T[i]))
Q[i]=np.prod(q[i][:])
我只用了两到三天的python。我甚至不知道python的所有内联操作。我有两个问题: 1)不使用for循环(比如只使用sum或no.prod)我能做到吗? 2)我的for循环不工作了。我在互联网上搜索,但没有找到任何好的for循环引用这种类型的数学运算。提前感谢您的帮助。
https://stackoverflow.com/questions/27735182
复制相似问题