function value = pendulum_nonlinear_prime ( t, y )

%*****************************************************************************80
%
%% pendulum_nonlinear_prime: right hand side of the nonlinear pendulum test.
%
  alpha = 0.05;
  alpha = 0.0;
  u = y(1);
  v = y(2);

  dudt = v;
  dvdt = - alpha * v - sin ( u );

  value = [ dudt; dvdt ];

  return
end

