.0*r/R; return exp(-0.5*p) * 1.0/(9.0*sqrt(3.0))   * (6.0 - 6.0*p + p*p);                }
    float R31(float r) { float p = 3.0*r/R; return exp(-0.5*p) * 1.0/(9.0*sqrt(6.0))   * (4.0 - p);                          }
    float R32(float r) { float p = 3.0*r/R; return exp(-0.5*p) * 1.0/(9.0*sqrt(30.0))  * p*p;                                }
    float R40(float r) { float p = 3.0*r/R; return exp(-0.5*p) * 1.0/(96.0)            * (24.0 - 36.0*p + 12.0*p*p - p*p*p); }
    float R41(float r) { float p = 3.0*r