function yprime = spring(t,y)
%
m = 1;
k = 1;

yp(1) = y(2);
yp(2) = -k*y(1)/m;
yprime = yp';