#ODE for BS11 ## note that Dz = (d/dz) operator (16*z^3-z)*Dz^2+(48*z^2-1)*Dz+16*z