double euler ( double t0, double u0, double dt, double f ( double t, double u ) )
{
  double u1;
  
  u1 = u0 + dt * f ( t0, u0 );
  
  return u1;
}