“…The motion equations were integrated via Lagrange planetary equations, Equation 1, 2, 3 and 4, which are a system of nonlinear differential equations in Keplerian terms (orbital elements) that describe the motion of bodies orbiting a central body (in this case Earth). In [8], the Lagrange planetary equations restricted to the models are Where M is the average anomaly and R is the perturbing function given in [5], modeled to orbital parameters considering elliptical and inclined orbits, and with reduced degree of freedom, where is the sum of the main forces that perturb the orbital motion of the debris-sail system, represented as In Maple, line graph simulations lasted around 5 hours for plotting individually, using a single GPU capable of 2.3 GHz average processing time over four independent cores (AMD Ryzen 7 3700U). For numerical simulations, values such as mass, equatorial radius, zonal coefficients and mean radius of the Earth, made available by the online platform of NASA "Jet Propulsion Laboratory".…”