Computing feasible trajectories for control-constrained systems: the PVTOL example