A numerical method for the optimisation of dual-wavelength Bragg mirro
rs is presented. It takes into account the dispersion of the refractiv
e index, and can be applied to generate two highly separated stopbands
. The method is tested in a configuration useful for doubly-resonant s
econd harmonic generation (stopbands at lambda and 2 lambda), a situat
ion for which simple analytical methods fail.