 * cos(phi); }
    float Y31_im(float mu, float phi)  { return -0.21545345607610045  * P31(mu) * sin(phi); }
    float Y30_re(float mu, float phi)  { return  0.7463526651802308   * P30(mu); }
    float Y30_im(float mu, float phi)  { return 0.0; }
    float Y3m1_re(float mu, float phi) { return -Y31_re(mu, phi); }
    float Y3m1_im(float mu, float phi) { return  Y31_im(mu, phi); }
    float Y3m2_re(float mu, float phi) { return  Y32_re(mu, phi); }
    float Y3m2_im(float mu, float phi) { return -Y32_im(mu, p