Re: anything faster than Solve[] for solving sets of polynomial equations symbolically?
- To: mathgroup at smc.vnet.net
- Subject: [mg118622] Re: anything faster than Solve[] for solving sets of polynomial equations symbolically?
- From: dantimatter <google at dantimatter.com>
- Date: Wed, 4 May 2011 19:47:36 -0400 (EDT)
Hi Bobby,
> eqs = {-k5 x1 x3 + k6 x4 == 0,
> k1 x1 - k4 x2 - 2 k2 x2^2 + 2 k3 x3 + k7 x4 == 0,
> k2 x2^2 - k3 x3 - k5 x1 x3 + k6 x4 == 0, x1 + x4 - Xtot == 0};
> nonNeg = Thread[{k1, k2, k3, k4, k5, k6, k7, Xtot} >= 0];
> real = Element[vars = {x1, x2, x3, x4}, Reals];
> problem = Flatten@{eqs, nonNeg};
> Solve[problem, Reals] // FullSimplify
I'm getting an error message that 'Xtot >= 0 is not a well-formed equation'. Is this something that you encountered as well? Also, is your Solve[] function call missing a reference to the variables, or am I just confused?
Thanks!
Dan