我试图在google中找到这些方程系统的解。我已经允许我的程序运行超过6个小时,但没有。我不知道怎么回事。
Python代码:
from sympy import *
from sympy import sqrt
init_printing(pretty_print=True)
var('x y z w')
var('a b1 s b2 o1 o2 c ', positive=True)
f=x*(1-x)-(y+z)
g=a*x*y-y*w-b1*y
h=s*x*z-z*w-b2*z
j=o1*y*w+o2*z*w-c*w
solve([f, g, h, j], (x, y, z, w)) #Also try with nonlinsolve and is the same issue.
发布于 2022-01-25 13:06:18
尝试使用manual=True
标志。它在大约一秒钟内给出了7种解决方案。
>>> solve([f, g, h, j], (x, y, z, w), manual=True)
[(0, 0, 0, 0),
(1, 0, 0, 0),
(b1/a, b1*(a - b1)/a**2, 0, 0),
(b2/s, 0, b2*(-b2 + s)/s**2, 0),
((b1 - b2)/(a - s),
(c*(a - s)**2*(o1 - o2) + o2*(-a*b1*o1 + a*b2*o1 + b1**2*o1 - 2*b1*b2*o1 +
b1*o1*s + b2**2*o1 - b2*o1*s + c*(a - s)**2))/(o1*(a - s)**2*(o1 - o2)),
(a*b1*o1 - a*b2*o1 - b1**2*o1 + 2*b1*b2*o1 - b1*o1*s - b2**2*o1 + b2*o1*s -
c*(a - s)**2)/((a - s)**2*(o1 - o2)),
(a*b2 - b1*s)/(-a + s)),
(1/2 - sqrt(-4*c + o2)/(2*sqrt(o2)), 0, c/o2,
-b2 + s/2 - s*sqrt(-4*c + o2)/(2*sqrt(o2))),
(1/2 + sqrt(-4*c + o2)/(2*sqrt(o2)), 0, c/o2,
-b2 + s/2 + s*sqrt(-4*c + o2)/(2*sqrt(o2)))]
https://stackoverflow.com/questions/70837238
复制相似问题