@@ -55,7 +55,7 @@ const m0 = 1 # initial mass
5555const vmax = 0.1 # maximal authorized speed
5656const mf = 0.6 # final mass to target
5757
58- ocp = @def begin # definition of the optimal control problem
58+ goddard = @def begin # definition of the optimal control problem
5959
6060 tf ∈ R, variable
6161 t ∈ [t0, tf], time
@@ -98,7 +98,7 @@ nothing # hide
9898We then solve it
9999
100100``` @example main-goddard
101- direct_sol = solve(ocp ; grid_size=250)
101+ direct_sol = solve(goddard ; grid_size=250)
102102nothing # hide
103103```
104104
@@ -192,10 +192,10 @@ ub(x) = -(F0⋅g)(x) / (F1⋅g)(x) # boundary control
192192μ(x, p) = H01(x, p) / (F1⋅g)(x) # multiplier associated to the state constraint g
193193
194194# Flows
195- f0 = Flow(ocp , (x, p, tf) -> u0)
196- f1 = Flow(ocp , (x, p, tf) -> u1)
197- fs = Flow(ocp , (x, p, tf) -> us(x, p))
198- fb = Flow(ocp , (x, p, tf) -> ub(x), (x, u, tf) -> g(x), (x, p, tf) -> μ(x, p))
195+ f0 = Flow(goddard , (x, p, tf) -> u0)
196+ f1 = Flow(goddard , (x, p, tf) -> u1)
197+ fs = Flow(goddard , (x, p, tf) -> us(x, p))
198+ fb = Flow(goddard , (x, p, tf) -> ub(x), (x, u, tf) -> g(x), (x, p, tf) -> μ(x, p))
199199nothing # hide
200200```
201201
0 commit comments