首页 文章

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

提问于
浏览
0

我很新见,并享受pythonic语法 . 但是,我遇到了一个问题,无法通过sympy解决,但可以在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

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

我用过的用语是 solve([eq1,eq2,eq3],[b,n,k])

谢谢你的帮助!

1 回答

  • 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))))]
    

相关问题