.0*x*x - 3.0)*x/2.0;                }

    // l=0 wavefunction (s-orbital)
    float Y00_re(float mu, float phi) { return 0.28209479177387814 * P00(mu); }
    float Y00_im(float mu, float phi) { return 0.0; }

    // l=1 wavefunctions (p-orbitals)
    float Y11_re(float mu, float phi)  { return -0.3454941494713355 * P11(mu) * cos(phi); }
    float Y11_im(float mu, float phi)  { return -0.3454941494713355 * P11(mu) * sin(phi); }
    float Y10_re(float mu, float phi)  { return  0.4886025119029199 * P10(mu); }