#include <cmath>

 double sph_bessel(unsigned n, double x)  // C++17

Liefert Wert der sphärischen Bessel-Funktion (1. Art) $j_\nu(x) = \sqrt{\frac{\pi}{2x}} J_{\nu+1/2}(x)$ für $x \geq 0$.

Rückgabewert: $j_\nu(x)$.

sph_legendre(), sph_neumann().

sph_bessel.cpp

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