-----maxima-bounces at math.utexas.edu wrote: -----
>is there any functionality in which Maxima is better than its rival?
For lots of reasons (good reasons, I believe), I'm uncomfortable with most
such comparisons. But you might try some experiements similar to:
Maxima:
(%i1) load(orthopoly)$
(%i2) pochhammer_max_index : 1000$
(%i3) jacobi_p(290,5,28,0.99);
(%o3) interval(-51721.80913981711,2.4071480684073597*10^-5)
(%i4) float(jacobi_p(290,5,28,99/100));
(%o4) -51721.80913981827
With Maple 10, try:
with(orthopoly):
Digits := 15:
convert(P(290,5,28,99/100),float);
P(290,5,28,0.99);
To get 6 digits of P(290,5,28,0.99) correct, you'll need to set Digits to
about 25. Maxima gets about 13 digits correct with just double floats.
Note: Maxima doesn't do big float evaluation of the Jacobi polynomials,
Maple does.
It's likely there are other significant differences between the way Maple
and Maxima
do floating point evaluation of the orthogonal polynomials, but I can't
think of them
right now.
Barton