'The Gouy phase.
Examples
>>> from sympy.physics.optics import BeamParameter
>>> p = BeamParameter(530e-9, 1, w=1e-3)
>>> p.gouy
atan(0.53/pi)'

@property
def gouy(self):

    return atan2(self.z, self.z_r)
