There is an old open PR: http://github.com/sympy/sympy/pull/1026
which tries to make orthogonal polynomials work for symbolic order 'n'. An open question is when to apply "simplifications" like: L_n(-x) ---> (-1)**n * L_n(x) L_{-n}(x) ---> L_{n-1}(x) Do we want to have this: a) at construction time b) only when calling expand_func Both approaches have their obvious advantages. Please see the commits: a) a958cb78c2bfd6d6e94dd05123f43195f2cdbe83 and b) 659afc3474b146814279d39840c61f4c6589071c to get an implementation of both approaches. -- You received this message because you are subscribed to the Google Groups "sympy" group. To view this discussion on the web visit https://groups.google.com/d/msg/sympy/-/eDJFpg7fz5cJ. To post to this group, send email to sympy@googlegroups.com. To unsubscribe from this group, send email to sympy+unsubscr...@googlegroups.com. For more options, visit this group at http://groups.google.com/group/sympy?hl=en.