function wprime=sys(t,w)
%
% wprime(1) is dx/dt = -y
% wprime(2) is dy/dt = x 
%
wprime = [ -w(2)
			w(1)];