code: def Jacobi(self,*args,**kwargs):
          
          if not kwargs.get('pot',None) is None: kwargs['pot']= flatten_potential(kwargs.get('pot'))
          _check_consistent_units(self,kwargs.get('pot',None))
          out= self._orb.Jacobi(*args,**kwargs)
          if not isinstance(out,float) and len(out) == 1: return out[0]
          else: return out
docs: NAME:
  
             Jacobi
  
          PURPOSE:
  
             calculate the Jacobi integral E - Omega L
  
          INPUT:
  
             t - (optional) time at which to get the Jacobi integral (can be Quantity)
  
             OmegaP= pattern speed (can be Quantity)
             
             pot= potential instance or list of such instances
  
             vo= (Object-wide default) physical scale for velocities to use to convert (can be Quantity)
  
             use_physical= use to override Object-wide default for using a physical scale for output
  
          OUTPUT:
  
             Jacobi integral
  
          HISTORY:
  
             2011-04-18 - Written - Bovy (NYU)