code: def calcEL(self,**kwargs):
                                     
          E,L= calcELAxi(self._R,self._vR,self._vT,self._pot)
          if self._gamma != 0.:
              #Adjust E
              E-= self._vT**2./2.
              L= m.fabs(L)+self._gamma*self.Jz(**kwargs)
              E+= L**2./2./self._R**2.
          return (E,L)
docs: NAME:
             calcEL
          PURPOSE:
             calculate the energy and angular momentum
          INPUT:
             scipy.integrate.quadrature keywords
          OUTPUT:
             (E,L)
          HISTORY:
             2012-07-26 - Written - Bovy (IAS)