A new iteration method is presented for achieving quantum optimal control over the expectation value of a positive definite operator. Theoretical analysis shows that this new algorithm exhibits quadratic and monotonic convergence. Numerical calculations verify that for this new algorithm, within a few steps, the optimized objective functional comes close to its converged limit.