function xdot = cone_eqns(t,x)

%  Integration of equations of
%  motion for particle in circular cone...
%
%  x(1)=r, x(2)=r_dot, x(3)=theta

%global r1

xdot = zeros(3,1);

xdot(1) = x(2);
xdot(2) = 0.01/(x(1))^3 - 4.2478546;
xdot(3) = 0.2/(x(1))^2;