Recherche:

OUI à cause de toutes les remarques ci-dessus...

Et re-OUI encore une fois à cause des indices des tableaux..

En plus de tout ce qu'on a dit au dessus, je signale que en C les tableaux commencent à 0...

Donc :



for (j=2;j<=n;j++) {
gam[j]=c[j-1]/bet;
bet=b[j]-a[j]*gam[j];
if (bet ==...

fprintf(fp,"%d %f\n",i,u[i]);


ou



fprintf(fp,"%d %g\n",i,u[i]);