double f3 ( double t, double u )
{
  double dudt;
  
  dudt = u * cos ( t );
  
  return dudt;
}