function dx = harmonic_oscillator(t,X) %input parameters: time and the states
%times has to be a scalar
%state has to be a Nx1 vector

%output: the derivatives Nx1 vector


%implements the differential equation of a harmonic oscillator