Applying this algorithm for the case discussed above, i . e . a polar orbit ( i = 90 \ text { deg } \, ) with p = 7200 \ text { km } \, ignoring all other perturbing forces then the J _ 2 \, and the J _ 3 \, forces for the numerical propagation one gets exactly the same optimal average eccentricity vector as with the " classical theory ", i . e . ( 0 \, \ 0.001036 ) \,.