Hi all,

In SymPy is there any way of stopping orthogonal polynomials from
automatically simplifying down to concrete expressions?  So for:

  sy.jacobi_normalized(2, 0, 1, x)

to give:

  sqrt(6)/3*jacobi(2, 0, 1, x)

as opposed to:

  sqrt(6)*(5*x**2/2 - x - 1/2)/2

I have a use-case where it is highly desirable for polynomials to be
kept in the above form (as it permits more efficient/numerically stable
evaluations later on).  Other than keeping n as a Symbol is this
currently possible?  If not---and I suspect, from a quick glance over
the source code that it isn't!)---is there wider interest in such a feature?

Regards, Freddie.

Attachment: signature.asc
Description: OpenPGP digital signature

Reply via email to