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

Table 3

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

𝑘	𝑗	𝑚=0	𝑚=1
	0	0	0
2	±1	0	−18𝑐𝑠
	±2	0	−18𝑠^2
	0	0	0
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	−916(3−30𝑐^2+35𝑐^4)	458(1−8𝑐^2+7𝑐^4)
	±1	−458𝑐(3−7𝑐^2)𝑠	±458𝛾(1±7𝑐−14𝑐^2)𝑠
4	±2	−1358(1−7𝑐^2)𝑠^2	1354𝛾^2(1∓7𝑐+7𝑐^2)
	±3	−9454𝑐𝑠^3	±9454𝛾^2(1∓2𝑐)𝑠
	±4	−9454𝑠^4	−9452𝛾^2𝑠^2
	0	0	0
	±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