double f1 ( double t, double u )
{
  double dudt;
  
  dudt = 0.05 * u;
  
  return dudt;
}