function [ ydot,vdot] = twoD( y,v )

global om0 gamma


ydot = v;
vdot = -om0^2*y - 2*gamma*v;


end