This work describes an implementation of an iterative algorithm to design two-zone binary phase diffractive optical elements (DOEs) which produce a specific irradiance distribution in both the axial and transverse directions under high numerical aperture focusing. The full vector character of a polarized beam is rigorously accounted for using the vector diffraction integrals. The axial and tran verse irradiance distributions are calculated for applications involving multiphoton direct la er writing, specifically, a novel method known as resolution augmentation through photo induced deactivation or, RAPID. In the RAPID method two beams are used, one for excitation and the other for de-excitation. By having both the excitation and de-excitation beams, the resolution achievable using the RAPID method is maximized.