/* ATK IV NUMERICAL PROGRAMMING (2006) */ /* Excercise 2.3a - Recurrence relations */ #include #include 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; }