High-order expansion of the solution of preliminary orbit determination problem