jacobi_p expansion?



>>>>> "Barton" == Barton Willis <willisb@unk.edu> writes:

    Barton> (%i2) load("l:/orthopoly-0.94/orthopoly.lisp")$
    Barton> (%i3) jacobi_p(57,1/2,1/3,0.52);
    Barton> (%o3)          interval(0.06384839164446, 2.159480417194154E-14)
    Barton> (%i4) orthopoly_returns_intervals : false$

Interval arithmetic?  Or is that just an some other estimate of the
error?

    Barton> Orthopoly doesn't do anything like 

    Barton>     gen_laguerre(n,0,x) = laguerre(n,x)

    Barton> Someday maybe ... But in view of RJF's comments about
    Barton> automatic expansion, I should think a few things through
    Barton> before I do anything like this.

Yeah, but that one is probably pretty safe.  But there ought to be
some way to tell maxima to simplify it when asked to by the user.

Ray