我使用Mathematica和Python计算了以下内容。你知道吗
Mathematica使用以下代码
f[x_] = a*b/(a - b)^2*Exp[-r*x] (Exp[-b*x] - Exp[-a*x]) (a*Exp[-b*x] - b*Exp[-a*x])
Assuming[{a > 0, b > 0, r > 0}, Integrate[f[x], {x, 0, \[Infinity]}]]
给出了一个相当不错的结果:
但是,下面的Python(带SymPy)代码
from sympy import *
init_printing()
x = symbols('x')
a, b, r = symbols('a b r', positive=True)
fun = a*b/((a-b)**2) * exp(-r*x) * (exp(-b*x) - exp(-a*x)) * (a*exp(-b*x) - b*exp(-a*x))
simplify(integrate(fun, (x, 0, oo)))
产生一个相当混乱的结果:
为了在Mathematica中获得相同的结果,Python代码中缺少了什么?还是有可能?你知道吗
函数^{} 可用于取消分数:
给予
相关问题 更多 >
编程相关推荐