for vectors
and scalars
and a cost functional
to be minimized for
and
. By considering the problem as
, singular perturbation concepts can be used to compute solutions consisting of bang-bang controls followed by singular arcs. The procedure further develops a numerical technique proposed by Jacobson, Gershwin, and Lele [18], as well as additional analytic methods developed by other authors.