function dydt = osc_td(t,y)
% harmonic oscillator system
k=1;
m=2;
%%  we could use
% dydt(1) = y(2);
% dydt(2) = -(k/m)*y(1);
% dydt=dydt(:)  % this makes dydt a column vector as required by odexx
%% or we could just use 
dydt=[y(2); -(k/m)*y(1)];
end

