Groups | Search | Server Info | Keyboard shortcuts | Login | Register [http] [https] [nntp] [nntps]


Groups > comp.soft-sys.math.mathematica > #2085

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

Path csiph.com!x330-a1.tempe.blueboxinc.net!usenet.pasdenom.info!aioe.org!feeder.news-service.com!216.196.110.146.MISMATCH!border3.nntp.ams.giganews.com!nntp.giganews.com!ramfeed-1.ams.xsnews.nl!feed.xsnews.nl!border-1.ams.xsnews.nl!206.55.64.80.MISMATCH!newspump.sol.net!posts.news.twtelecom.net!nnrp3.twtelecom.net!not-for-mail
From DrMajorBob <btreat1@austin.rr.com>
Newsgroups comp.soft-sys.math.mathematica
Subject Re: anything faster than Solve[] for solving sets of polynomial equations symbolically?
Date Wed, 4 May 2011 10:34:24 +0000 (UTC)
Organization Steven M. Christensen and Associates, Inc and MathTensor, Inc.
Sender steve@smc.vnet.net
Approved Steven M. Christensen <steve@smc.vnet.net>, Moderator
Message-ID <ipra3g$tn$1@smc.vnet.net> (permalink)
Lines 43
NNTP-Posting-Date 04 May 2011 09:35:13 GMT
NNTP-Posting-Host 98abc764.news.twtelecom.net
X-Trace DXC=J@3m9n>SEbDBZ`KRoIR<j@C_A=>8kQj6MY;@_o827nGC0D3gFL4k0fCEFiONJ7[GoFSW0L<2j=UTG
X-Complaints-To abuse@twtelecom.net
Xref x330-a1.tempe.blueboxinc.net comp.soft-sys.math.mathematica:2085

Show key headers only | View raw


This quickly catalogs the solutions, though I'm not fond of the output  
format:

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

Bobby

On Tue, 03 May 2011 04:47:36 -0500, dantimatter <google@dantimatter.com>  
wrote:

>
> Hi Everyone,
>
> I've got a set of equations and I'd really like to determine the number  
> of real solutions in terms of the unspecified parameters.   The  
> equations are:
>
> 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};
>
> where x1, x2, x3, and x4 are the variables I'd like to solve for and k1,  
> k2, k3, k4, k5, k6, k7, and Xtot are real numbers that are greater than  
> or equal to zero.   I've tried doing
>
> Solve[eqs, {x1, x2, x3, x4}, Reals]
>
> but after days and days I still have no solution.  Perhaps there's a  
> faster or better way to do this?
>
> Thanks,
> Dan
>


-- 
DrMajorBob@yahoo.com

Back to comp.soft-sys.math.mathematica | Previous | Next | Find similar | Unroll thread


Thread

Re: anything faster than Solve[] for solving sets of polynomial equations symbolically? DrMajorBob <btreat1@austin.rr.com> - 2011-05-04 10:34 +0000

csiph-web