In this work, by introducing Darboux operator in evolutionary computing frame, we propose a novel analytic evolutionary algorithm to obtain exact higher-order iteration solutions for model equation. We construct the first-, second- and third-order solutions of a generalized Schrödinger equation by applying this algorithm. The higher-order solutions not only retain the basic features of the lower-order cases, but also become more abundant than the lower-order cases.