A versatile method for theoretical designing of laser pulses for coherent control of molecular dynamics based on evolutionary algorithm is presented. Extended normal distribution crossover (ENDX) together with minimal generation gap (MGG) is employed for laser parameter optimization. The method is applied to the selective excitation problem of isotope molecules. (C) 2004 Elsevier B.V. All rights reserved.