void eulero (double& t, double h, int N, double *x, double* (funzione)(double*, double)) {
  double *d = (*funzione)(x,t);
  for (int i=0; i<N; i++) x[i]+=(d[i]*h);
  t+=h;
  delete d;
}