using Thesis
using PlotlyJS: savefig

a = ∛(Sun.μ * ( year/2π )^2 )
start = [ oe_to_xyz([ a, 0., 0., 0., 0., 0.5 ], Sun.μ); 10_000. ]
path = Thesis.prop(zeros(100,3), start, bepi, year)
p = Thesis.plot([path],
                colors=["#F0F"],
                title="Single LaGuerre-Conway-Propagated Orbit",
                mode="light")
savefig(p, "LaTeX/fig/laguerre_plot.png")

t = 5year
spiral_sc = Sc("test", 1000., 2000., 0.0005, 50, 0.9)
start = [ oe_to_xyz([ 0.7a, 0.05, 0.1, 0., 0., 0.5 ], Sun.μ); 100_000. ]
spiral_thrust = spiral(0.2, 100_000, start, spiral_sc, t)
spiral_path = Thesis.prop(spiral_thrust, start, spiral_sc, t)
p = Thesis.plot([spiral_path],
                colors=["#F0F"],
                title="Spiral Orbit",
                mode="light")
savefig(p, "LaTeX/fig/spiral_plot.png")