Numerical Method of Symplectic State Transition Matrix and Application to Fully Perturbed Earth Orbit