void f4 ( double t, double u[], double dudt[] )
{
  double alpha = 0.01, f, r;
  
  r = u[0];
  f = u[1];
  
  dudt[0] = 2 * r - alpha * r * f;
  dudt[1] =   - f + alpha * r * f;
  
  return;
}