Lots of research needed but, throwing this out there to investigate:
[ ] Can we leverage ModelingToolkit so that specific known differential equations can dispatch to OrthogonalPolynomials? Is this better taken care of with ApproxFun.jl? Will the interface live in another repo?
[ ] When/where does contributing to ModelingToolkit help solve this, if at all?
Lots of research needed but, throwing this out there to investigate: