“…In (10), V and V are minimum and maximum flight speeds; in (11), 𝛼 is the angle of attack at which vehicle starts to stall; in (12), 𝜙 is the maximum bank angle; in (13), 𝜓 is the maximal change of yaw angle within a route planning cycle (such as t ′ + 1 → t ′ + 2) and it reflects the limit by vehicle's side overload. In (14), T(⋅) is the available maximal thrust associated with z and Mach number Ma (t), which is subject to (15); in (15), 𝜅 is adiabatic exponent, ℜ is real gas constant for air, Te (z) is the air temperature at altitude z. At altitude z, T(⋅) is represented as a quadratic function of Ma (t) in ( 16) 33 ; c 1 T (z), c 2 T (z) and c 3 T (z) are parameters with z. Initial/final condition for state/control variables: The components of state variable s (t) include x (t), y (t), V (t) and 𝜓(t); the components of control variable u (t) include the rest of variables except T t f .…”