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

#define DELTAT 0.001

int main() {
  FILE *f;
  float y, t=0., deltay, tmax;
  printf ("Inserisci il valore iniziale ed il tempo di integrazione:\n");
  scanf ("%f %f",&y,&tmax);
  if ( !(f=fopen("output.dat","w")) ) return -1;
  while (t<tmax)
  {
    deltay = sin(y)*DELTAT;
    y+=deltay;
    fprintf (f,"t=%f, y=%f\n",t,y);
    t+=DELTAT; 
  }
  fclose(f);
  return 0;
}