Singular optimal control plays an important role in process engineering, including optimal operation of batch and semi-batch reactions. However, for many practical applications, accurate solution of singular optimal control profiles is still an open issue. In particular, numerical optimization must deal with an ill-conditioned problem that often leads to very slow convergence or failure. Starting from the nested approach in