% This implements the simple harmonic oscillator ODE 
% with unit angular frequency,
%     d/dt (x,y) = (-y,x).
%
% CCK 2018-Mar-10 capitalized function name to avoid namespace collision in Matlab.
function dxdt = Harmonic(x,t)

dxdt = [ -x(2), x(1) ];