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.
signature.asc
Description: OpenPGP digital signature