Question

我很新见，并享受pythonic语法。但是，我遇到了一个问题，不能通过同情来解决，但可以在Maple中轻松解决。

所以，我有以下系统：

0.0165 * exp( -2.0405-0.33*b-0.5*n)+0.031 * exp(-4.164-0.62*b-0.5*n)=k*p
0.025 * exp( -2.0405-0.33*b-0.5*n) +0.025 * exp(-4.164-0.62*b-0.5*n)=5*k
2*p=p*b+5*n

我需要根据b来解决n，k和p。我可以很容易地在Maple中解决这个问题，但是使用sympy，它会花费很长时间并在最后因为耗尽RAM而崩溃。 Maple可以提供精确的符号解决方案。

我使用的情感代码是solve([eq1,eq2,eq3],[b,n,k])

感谢您的帮助！

Answer 1

使用rational = False标志：

>>> print filldedent(solve([eq1,eq2,eq3],[b,n,k], rational=False))

[(-3.44827586206897*log((-0.000335859591913345*p +
0.00110833665331404)/(2.24927535168052e-5*p - 0.000139455071804192)) +
2, 0.689655172413793*p*log((-0.000335859591913345*p +
0.00110833665331404)/(2.24927535168052e-5*p - 0.000139455071804192)),
(0.000335859591913345*((-0.000335859591913345*p +
0.00110833665331404)/(2.24927535168052e-5*p -
0.000139455071804192))**1.13793103448276 +
2.24927535168052e-5*((-0.000335859591913345*p +
0.00110833665331404)/(2.24927535168052e-5*p - 0.000139455071804192))**
2.13793103448276)*exp(-0.344827586206897*p*log((-0.000335859591913345*
p + 0.00110833665331404)/(2.24927535168052e-5*p -
0.000139455071804192))))]

非线性方程组可以在Maple中求解，但不能在同一时间解决

1 个答案: