We use penalty methods coupled with unitary exponential operator metho
ds to solve the optimal control problem for molecular time-dependent S
chrodinger equations involving laser pulse excitations. A stable numer
ical algorithm is presented which propagates directly from initial sta
tes to given final states. Results are reported for an analytically so
lvable model for the complete inversion of a three-state system.