Optimization of action integrals is suggested as an alternative to initial value algorithms for the study of time-evolving quantum systems. The finite element method provides a convenient and rapidly convergent scheme for some simple cases where analytical solutions are available. Two level systems are given particular attention.