We determine solutions for the Lévy-Leblond operator or parabolic Dirac operator in terms of hypergeometric functions and spherical harmonics. We subsequently generalise our approach to a wider class of Dirac operators depending on 4 parameters. This is a report on our recent paper arXiv:1911.01744.