I suggest to include the parabolic cylinder function D such as implemented in scipy: [scipy.special.pbdv](https://docs.scipy.org/doc/scipy/reference/generated/scipy.special.pbdv.html)