dsolve({diff(diff(x(t),t),t) = k*x(t)/(x(t)^2+y(t)^2)^(3/2), diff(diff(y(t),t),t) = k*y(t)/(x(t)^2+y(t)^2)^(3/2)})