double f2 ( double t0, double u )
{
  double dudt;
  
  dudt = 0.05 * u;
  
  return dudt;
}