RESOLUTION NUMERIQUE D'UNE EQUATION DIFFERENTIELLE PAR LA METHODE D'EULER

#include<stdio.h>
#include<math.h>
#define g 9.81

void main(void)
{
	float tau,w,l;
	int nit,i;
	FILE *fichier;
	fichier=fopen("resultat.txt","wt");
	float x;
	float y;
	printf("\Entrer teta0:");
	scanf("%f",&x);
	y=0;
	i=0;
	printf("\nEntrer le nombre d'iterations:");
	scanf("%i",&nit);
	printf("\nEntrer le pas (valeur suggerée 0.01):");
	scanf("%f",&tau);
	printf("\Entrer la longueur du fil:");
	scanf("%f",&l);
	w=sqrt(g/l);
	printf("\nLa pulsation w0=%f",w);
	for(i=1;i<=nit;i++)
	{
		x=x+tau*y;
		y=y+tau*((-w*w)*sin(x));
		fprintf(fichier,"%f %f %f\n",i*tau,x,y);
	}
	fclose(fichier);
}