Files
thesis/julia/lamberts_fig.jl
2022-03-13 23:34:38 -06:00

18 lines
512 B
Julia

using PlotlyJS: savefig
r1 = oe_to_xyz([Earth.a, 0.6, 0.001, 0., 0., deg2rad(10.)], Sun.μ)[1:3]
r2 = oe_to_xyz([Earth.a, 0.6, 0.001, 0., 0., deg2rad(140.)], Sun.μ)[1:3]
tof = 0.5year
v1 = Thesis.lamberts1(r1,r2,tof)[1]
v2 = Thesis.lamberts2(r1,r2,tof)[1]
state1 = [ r1; v1; 1000. ]
state2 = [ r1; v2; 1000. ]
path1 = prop(state1, tof)
path2 = prop(state2, tof)
plot([path2, path1],
title="Lambert's Problem Solutions",
colors=["#F00", "#0FF"],
labels=["Type I", "Type II"],
mode="light")