On Mon, Aug 25, 2008 at 5:32 PM, Carl Witty <[EMAIL PROTECTED]> wrote: > > 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] >
I installed phc into the public sage notebook servers, so you can try the above at https://sage.math.washington.edu:8102 -- William --~--~---------~--~----~------------~-------~--~----~ 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 -~----------~----~----~----~------~----~------~--~---
