A Galerkin′s method based numerical procedure is used to design photonic crystal fibers with zero dispersion wavelength at 800 nm. High accuracy is maintained by using a large enough area of calculation and sufficient basis functions. The zero dispersion wavelength is designed to be the operating wavelength of the widely used Ti∶sapphire femtosecond lasers, holding much promise for the applications of femtosecond laser.