High-Order Expansion of the Solution of Preliminary Orbit Determination Problem