function dx=simple(t, x)
%SIMPLE  simple function for testing phase plots

% Harmonic oscillator
dx = [x(2); -x(1)];