function output=rhs(t,y,dummy,omega)

output = [y(2);-omega^2*y(1)];