On Aug 24, 12:16 pm, Michael <[EMAIL PROTECTED]> wrote:
> I have a polynomial system of 50 equations in 50 unknowns. I would
> like
> to numerically solve this system (I'm interested in complex zeros).
> Seems to me that if I use sage's solve, it will sttempt to solve these
> algebraically.
>
> Are there any funcitons in sage for solving a polynomial or non-linear
> system
> numerically.
Other people have suggested phcpack, but not explained how to use it.
The following may or may not work for you (it depends on your
processor architecture and operating system); but if it does work,
it's actually quite easy.
First, install phc, by typing at the command line:
sage -i phc-2.3.39.p0
Then, from within Sage, you can do things like this:
sage: from sage.interfaces.phc import phc
sage: R2.<x,y> = PolynomialRing(QQ,2)
sage: start_sys = [x^6-y^2,y^5-1]
sage: sol = phc.blackbox(start_sys, R2)
sage: len(sol.solutions())
30
sage: sol.solutions()[0] # display the first solution
[0.309016994374947 + 0.951056516295154*I,
0.104528463267653 + 0.994521895368273*I]
Carl
--~--~---------~--~----~------------~-------~--~----~
To post to this group, send email to [email protected]
To unsubscribe from this group, send email to [EMAIL PROTECTED]
For more options, visit this group at
http://groups.google.com/group/sage-support
URLs: http://www.sagemath.org
-~----------~----~----~----~------~----~------~--~---
