In several optical systems, a specific Point Spread Function (PSF) needs to be generated. This can be achieved by shaping the complex field at the pupil. The Extended Nijboer-Zernike (ENZ) theory relates complex Zernike modes on the pupil directly to functions in the focal region. In this paper, we introduce a method to engineer a PSF using the ENZ theory. In particular, we present an optimization algorithm to design an extended depth of focus with high lateral resolution, while keeping the transmission of light high (over 60%). We also have demonstrated three outcomes of the algorithm using a Spatial Light Modulator (SLM).