/*  ATK IV NUMERICAL PROGRAMMING (2006)  */
/* Excercise 2.3a - Recurrence relations */

#include<stdio.h>
#include<math.h>

int main(void){

  float besselJ[21];
  float x = 1.0;
  int n;

  besselJ[0] = 0.7651976;
  besselJ[1] = 0.4400506;

  for(n=1;n<=20;n++){
    besselJ[n+1] = 2*n*besselJ[n]/x - besselJ[n-1];
  }

  for(n=0;n<=20;n++){
    printf("%f\n",besselJ[n]);
  }

  return 0;
}