A family of monotonically convergent algorithms is presented for solving a wide class of quantum optimal control problems satisfying an inhomogeneous integrodifferential equation of motion. The convergence behavior is examined using a four-level model system under the influence of non-Markovian relaxation. The results show that high quality solutions can be obtained over a wide range of parameters that characterize the algorithms, independent of the presence or absence of relaxation.