Re: [open-axiom-devel] [fricas-devel] Re: [Axiom-math] Selection of roots

2008-11-03 Thread Bill Page
On Mon, Nov 3, 2008 at 7:55 PM, Alejandro Jakubi wrote: > Billl Page wrote: >> One thing that would be very useful is if someone like you were able >> to contribute some documentation like this that could be added to >> the axiom-wiki and perhaps also the Axiom book. > > May be. But wait. I am slow

Re: [open-axiom-devel] [fricas-devel] Re: [Axiom-math] Selection of roots

2008-11-03 Thread Martin Rubey
Alejandro Jakubi <[EMAIL PROTECTED]> writes: > Fine. Now, how I recover the algebraic expressions, as produced by > 'radicalSolve', associated with such labels like %A1, %A2 ? perhaps (63) -> mainDefiningPolynomial(first allRootsOf(numer fp)$RECLOS(FRAC INT)) 4 3 3 2 (63) ?

Re: [open-axiom-devel] [fricas-devel] Re: [Axiom-math] Selection of roots

2008-11-03 Thread Bill Page
On Mon, Nov 3, 2008 at 5:00 PM, Alejandro Jakubi wrote: > > Thank you Bill for showing how to use 'RealClosure' in this example. You are welcome. > In fact, I have been looking a pair of days on how to use it, but, > because of the lack of proper documentation the best that I could > get was: > >

Re: [open-axiom-devel] [fricas-devel] Re: [Axiom-math] Selection of roots

2008-11-03 Thread Martin Rubey
Alejandro Jakubi <[EMAIL PROTECTED]> writes: > Now, my next question is about selecting the two real roots out of the four > produced by 'radicalSolve' here: > > f:=(x^3+x^2-4*x-4)/(2*x^2+7*x-4) > fp:=differentiate(f,x) > radicalSolve(fp,x) > > The problem is that > > allRootsOf(fp)$RealClosure(

Re: [open-axiom-devel] [fricas-devel] Re: [Axiom-math] Selection of roots

2008-11-03 Thread Bill Page
On Mon, Nov 3, 2008 at 1:32 PM, Martin Rubey wrote: > > "Bill Page" <[EMAIL PROTECTED]> writes: > >> As I understand it AlgebraicNumber is supposed to be able to >> properly order the roots. > > No, that's RECLOS. > > Unfortunately, there is no coercion from AN to RECLOS, and this would > actually

Re: [open-axiom-devel] [Axiom-math] Selection of roots

2008-11-03 Thread Bill Page
On Thu, Oct 30, 2008 at 8:39 AM, Alejandro Jakubi wrote: > I wonder how it is done in Axiom the selection of roots of a polynomial > with some property. As in this example, select the positive roots out > of the list of three roots generated by: > > radicalSolve(p^3-p+1/10=0,p) > Since the discrim

Re: [open-axiom-devel] [fricas-devel] Re: [Axiom-math] Selection of roots

2008-11-03 Thread Martin Rubey
"Bill Page" <[EMAIL PROTECTED]> writes: > Hmmm... So what is the meaning of < in AN? Ok, I guess that it is just > whatever is exported by Expression Integer. The fact that the positive > roots appear < 0 while the negative roots appear > 0 must be just an > accident of some strange lexical orderi

Re: [open-axiom-devel] [fricas-devel] Re: [Axiom-math] Selection of roots

2008-11-03 Thread Bill Page
On Mon, Nov 3, 2008 at 1:32 PM, Martin Rubey wrote: > > Bill Page writes: > >> As I understand it AlgebraicNumber is supposed to be able to properly >> order the roots. > > No, that's RECLOS. > > Unfortunately, there is no coercion from AN to RECLOS, and this would actually > be tricky, since sqrt(

Re: [open-axiom-devel] [fricas-devel] Re: [Axiom-math] Selection of roots

2008-11-03 Thread Martin Rubey
"Bill Page" <[EMAIL PROTECTED]> writes: > As I understand it AlgebraicNumber is supposed to be able to properly > order the roots. No, that's RECLOS. Unfortunately, there is no coercion from AN to RECLOS, and this would actually be tricky, since sqrt(-3) is not allowed in RECLOS. Martin