Precise Analytical Computation of Frozen-Eccentricity, Low Earth Orbits in a Tesseral Potential

Table 2

Inclination polynomials 𝑝_{𝑘,𝑚,𝑗} in (22). We abbreviate 𝛾≡1±𝑐.

𝑘	𝑗	𝑚=0	𝑚=1
	0	0	0
2	±1	0	±18𝑐𝑠
	±2	0	∓18𝑠^2
	0	92(1−5𝑐^2)𝑠	−92(1−5𝑐^2)𝑠
3	±1	∓94(1∓𝑐)(1∓10𝑐−15𝑐^2)	∓94𝛾(1±10𝑐−15𝑐^2)
	±2	452(1∓𝑐)(1±3𝑐)𝑠	−452𝛾(1∓3𝑐)𝑠
	±3	∓1352(1∓𝑐)𝑠^2	∓1352𝛾𝑠^2
	0	0	0
	±1	±458𝑐(3−7𝑐^2)𝑠	−458𝛾(1±7𝑐−14𝑐^2)𝑠
4	±2	∓1358(1−8𝑐^2+7𝑐^4)	±1354𝛾^2(1∓7𝑐+7𝑐^2)
	±3	±9454𝑐𝑠^3	−9454𝛾^2(1∓2𝑐)𝑠
	±4	∓9454𝑠^4	∓9452𝛾^2𝑠^2
	0	−4516(1−14𝑐^2+21𝑐^4)𝑠	10532(1−9𝑐^2)𝑠^3
	±1	∓4532𝛾(1±28𝑐−42𝑐^2∓84𝑐^3+105𝑐^4)	±31564𝛾(1±6𝑐−15𝑐^2)𝑠^2
5	±2	−3158𝛾(1∓3𝑐−9𝑐^2±15𝑐^3)𝑠	315𝛾16^2(1∓12𝑐+15𝑐^2)𝑠
	±3	∓9458𝛾(1±6𝑐−15𝑐^2)𝑠^2	±315𝛾16^3(1∓3𝑐)(13∓15𝑐)
	±4	−28354𝛾(1∓5𝑐)𝑠^3	−28358𝛾^3(3∓5𝑐)𝑠
	±5	∓141754𝛾𝑠^4	∓141758𝛾^3𝑠^2