list of functions
- Subject: list of functions
- From: Barton Willis
- Date: Sat, 16 Jun 2007 21:31:44 -0500
-----maxima-bounces at math.utexas.edu wrote: -----
>To: "'Martin DRUON'" <martin.druon at wanadoo.fr>, <maxima at math.utexas.edu>
>From: "Richard Fateman" <fateman at cs.berkeley.edu>
>Sent by: maxima-bounces at math.utexas.edu
>Date: 06/16/2007 04:11PM
>Subject: Re: [Maxima] list of functions
>
>If you want to evaluate legendre polynomials at a floating-point number
>point, you are doing it a very poor way.
To illustrate:
(%i1) load(orthopoly)$
(%i2) orthopoly_returns_intervals : false$
(%i3) legendre_p(100,0.707);
(%o3) -0.08793237492018
(%i4) subst(x = 0.707,legendre_p(100,x));
(%o4) 3.5557176931913412*10^+14
(%i5) subst(x = 0.707,horner(legendre_p(100,x)));
(%o5) -6.5042014223259163*10^+10
(%i6) subst(x = 0.707, expand(legendre_p(100,x)));
(%o6) -7.3593884484921875*10^+10
(%i7) float(subst(x = rationalize(0.707),legendre_p(100,x)));
(%o7) -0.087932374920181
Suggestion: If you use orthopoly, compile it:
(%i2) load(orthopoly);
(%o2) C:/PROGRA~1/MAXIMA~3.0/share/maxima/5.12.0
/share/orthopoly/orthopoly.lisp
(%i3) compile_file(%)$
(%i4) load(orthopoly);
(%o4) C:/PROGRA~1/MAXIMA~3.0/share/maxima/5.12.0
/share/orthopoly/orthopoly.o
Barton