@@ -245,13 +245,13 @@ for p0 ∈ p0s # plot for each p₀ in p0s
245245 flow_p0 = ocp_flow((t0, tf_), x0, p0; saveat=tspan, callback=cbt)
246246
247247 T = tspan
248- X = flow_p0. state.(T)
249- P = flow_p0. costate.(T)
248+ X = state(flow_p0) .(T)
249+ P = costate(flow_p0) .(T)
250250
251- plot!(plt_x, T, X; color=:blue)
252- plot!(plt_p, T, P; color=:blue)
253- plot!(plt_u, T, u.(X, P); color=:blue)
254- plot!(plt_phase, X, P; color=:blue)
251+ plot!(plt_x, T, X; color=:blue)
252+ plot!(plt_p, T, P; color=:blue)
253+ plot!(plt_u, T, u.(X, P); color=:blue)
254+ plot!(plt_phase, X, P; color=:blue)
255255
256256end
257257
@@ -318,13 +318,13 @@ for (p0, label) ∈ zip(p0s, labels) # plot for each p₀ in p0s
318318 flow_p0 = ocp_flow((t0, tf), x0, p0; saveat=tspan)
319319
320320 T = tspan
321- X = flow_p0. state.(T)
322- P = flow_p0. costate.(T)
321+ X = state(flow_p0) .(T)
322+ P = costate(flow_p0) .(T)
323323
324- plot!(plt2_x, T, X; label=label)
325- plot!(plt2_p, T, P; label=label)
326- plot!(plt2_u, T, u.(X, P); label=label)
327- plot!(plt2_phase, X, P; label=label)
324+ plot!(plt2_x, T, X; label=label)
325+ plot!(plt2_p, T, P; label=label)
326+ plot!(plt2_u, T, u.(X, P); label=label)
327+ plot!(plt2_phase, X, P; label=label)
328328
329329end
330330
0 commit comments