Imaging depth of field is shallow in applications with high magnification and high numerical aperture, such as microscopy, resulting in images with in- and out-of-focus regions. Therefore, methods to extend depth of field are of particular interest. Researchers have previously shown the advantages of using freeform components to extend depth of field, with each optical system requiring a specially designed phase plate. In this paper we present a method to enable extended depth-of-field imaging for a range of numerical apertures using freeform phase plates to create variable cubic wavefronts. The concept is similar to an Alvarez lens which creates variable spherical wavefronts through the relative translation of two transmissive elements with XY polynomial surfaces. We discuss design and optimization methods to enable extended depth of field for lenses with different numerical aperture values by considering through-focus variation of the point spread function and compare on- and off-axis performance through multiple metrics.