function f=rpp(t,x)
global omega;
f=[x(2);-x(1)*omega];
end


