function dy = osc(t,y)
% equations of a harmonic oscillator 
	 k = 10; m = 1;
	 dy = zeros(2,1);    % a column vector
    dy(1) = y(2);
    dy(2) = -k/m*y(1);
return;