function dydx = pendulum_ode(x,y)
% dydx = pendulum_ode(x,y)
% comments including meanings of variables,
% math methods, your name and the date

% mms 1/3/16

dydx=[y(2);-3*sin(y(1))];
