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