function xdot = odetest()

g = 9.81;
mu = 0.25;
r0 = 5.0;
omega = 1.0;
t0 = 0.0;
tf = 10.0;
x = [r0 0 0 r0*omega];
xdot(1) = x(3);
xdot(2) = x(4);
xdot(3) = -mu*g*(x(3)-omega*x(2))/sqrt((x(3)-omega*x(2))^2 + (x(4)+omega*x(1))^2);
xdot(4) = -mu*g*(x(4)+omega*x(1))/sqrt((x(3)-omega*x(2))^2 + (x(4)+omega*x(1))^2);

[t,s] = ode45(@odetest,[t0,tf],x);
x = s(:,1);
y = s(:,2);
vx = s(:,3);
vy = s(:,4);