#include <cmath>

 double sph_legendre(unsigned n, unsigned m, double theta)  // C++17

Liefert Wert der Kugelflächenfunktion $Y_n^m(\theta,0) = (-1)^m \sqrt{\frac{(2n-1)(n-m)}{4\pi(n+m)!}} P_n^m(\cos \theta)$ für $0 \leq m \leq n$.

Rückgabewert: $Y_n^m(\theta,0)$.

sph_bessel(), sph_neumann().

sph_legendre.cpp

#include <cmath>
#include <iostream>
 
int main()
{
  std::cout << "# x n=0,m=0 n=1,m=0 n=1,m=1 n=2,m=0 n=2,m=1 n=2,m=2\n";
 
  for (int i = 0; i <= 200; ++i)
  {
    double x = 0.01*i;
    std::cout << x 
      << '\t' << std::sph_legendre(0, 0, x) 
      << '\t' << std::sph_legendre(1, 0, x) 
      << '\t' << std::sph_legendre(1, 1, x) 
      << '\t' << std::sph_legendre(2, 0, x) 
      << '\t' << std::sph_legendre(2, 1, x) 
      << '\t' << std::sph_legendre(2, 2, x) 
      << '\n';
  }
}