All the classical diffeq states, cyl_bessel, sph_bessel, airy, etc could offer (say for cyl_bessel)
_Tp J_deriv2() const;
_Tp N_deriv2() const;
using the state variables.
Do this for orthogonal polynomials too.
I added deriv() to those already (mostly).
All the classical diffeq states, cyl_bessel, sph_bessel, airy, etc could offer (say for cyl_bessel) _Tp J_deriv2() const; _Tp N_deriv2() const; using the state variables.
Do this for orthogonal polynomials too. I added deriv() to those already (mostly).