function dy = MEE380eq2(t,y)
%Pendulum differential equations
% G [m/s^2] acceleration due to gravity
% LENGTH [m] pendulum length
global G LENGTH
dy=zeros(2,1);
dy(1) = y(2);
dy(2) = -(G/LENGTH)*sin(y(1));