This work focuses on a novel, to the best of our knowledge, analytical form of the phase retardation function for achieving a uniform axial intensity of Bessel beams. Traditional methods of generating Bessel beams often result in significant oscillations in the intensity along the beam’s axial path, which limits their practical applications. However, the proposed phase retardation function in this study overcomes these limitations by ensuring consistent beam creation regardless of factors such as the beam waist size, wavelength, or axicon angle. By implementing the proposed spatial phase function, a fundamental Gaussian laser beam, thereby generating a Bessel beam with an elongated and constant axial intensity profile, supports our theoretical predictions. The functionality of this new phase retardation function was further scrutinized using different wavelengths and beam waist sizes to confirm that the axial intensity remained uniform profile. Additionally, when contrasting our phase function with those from earlier researches, it was observed that our findings are consistent with both theoretical models and experimental outcomes.