#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::cyl_bessel_i(0, x) 
      << '\t' << std::cyl_bessel_i(1, x) 
      << '\t' << std::cyl_bessel_i(2, x) 
      << '\n';
  }
}