function yp = pendulum(y)
%equations of motion for the pendulum

%   Copyright 2010 $Date: 2010/10/03 21:00:00 $

h1 = y(2);
h2 = -sin(y(1));
yp = [h1 h2]';
