
function ydot = pendnl(t,y)
% f(t,y) for nonlinear pendulum
% T. DeLillo, Math 451, Fall 2001
global wsq
ydot = [y(2); -wsq*sin(y(1))];
