void flame_exact ( int n, double t[n], double *y );
void flame_parameters ( double *t0_in, double *y0_in, double *tstop_in, 
  double *t0_out, double *y0_out, double *tstop_out  );
double lambert_w ( double x, int nb, int l );