code: def total_energy(self, x, v):
          
          warnings.warn("Use the energy methods on Orbit objects instead. In a future "
                        "release this will be removed.", DeprecationWarning)
  
          v = atleast_2d(v, insert_axis=1)
          return self.energy(x) + 0.5*np.sum(v**2, axis=0)
docs: Compute the total energy (per unit mass) of a point in phase-space
          in this potential. Assumes the last axis of the input position /
          velocity is the dimension axis, e.g., for 100 points in 3-space,
          the arrays should have shape (100,3).
  
          Parameters
          ----------
          x : array_like, numeric
              Position.
          v : array_like, numeric
              Velocity.