function dydt = stiff_prime ( t, y )

  dydt = 50.0 * ( cos ( t ) - y );

  return
end

