function ddt=RHS(y)

y1=y(1);
y2=y(2);

ddt=[y2; -y1];