#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';
  }
}