Trajectory optimization is a field of study in the aeronautical domain that has experienced a significant surge in interest over the past decades, both in civil and military applications. It represents an optimal control problem aimed at identifying the controls and the corresponding trajectory to achieve an objective while satisfying certain constraints. The objective of this thesis is to develop a bilevel algorithm capable of optimizing the trajectories of multiple aircraft F104 in formation flight, each represented by a 3 DOF point mass model. The identified trajectories must meet predefined waypoints for each aircraft, which are established by the user prior to the start of the optimization process.