MathGroup Archive 2011

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: anything faster than Solve[] for solving sets of polynomial equations symbolically?

  • To: mathgroup at smc.vnet.net
  • Subject: [mg118640] Re: anything faster than Solve[] for solving sets of polynomial equations symbolically?
  • From: DrMajorBob <btreat1 at austin.rr.com>
  • Date: Thu, 5 May 2011 05:25:12 -0400 (EDT)

The code runs here, exactly as you see it below.

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

$Version

"8.0 for Mac OS X x86 (64-bit) (February 23, 2011)"

Bobby

On Wed, 04 May 2011 18:47:36 -0500, dantimatter <google at dantimatter.com>  
wrote:

>
> 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
>
>


-- 
DrMajorBob at yahoo.com


  • Prev by Date: Re: Creating an audio
  • Next by Date: Re: Error bars
  • Previous by thread: Re: anything faster than Solve[] for solving sets of polynomial equations symbolically?
  • Next by thread: Re: anything faster than Solve[] for solving sets of polynomial equations symbolically?