docs: def phi(self):
          """Polar angle / inclination of this vector in radians
  
          Based on sperical coordinate space
          returns angle between this vector and the positive z-azis
          range: (0 <= phi <= pi)
          """
          return np.arctan2(np.sqrt(self.x**2 + self.y**2), self.z)