code: def angleR(self,**kwargs): #pragma: no cover
          
          if hasattr(self,'_angleR'):
              return self._angleR
          (rperi,rap)= self.calcRapRperi(**kwargs)
          if rap == rperi:
              return 0.
          TR= self.TR(**kwargs)[0]
          EL= self.calcEL(**kwargs)
          E, L= EL
          Rmean= m.exp((m.log(rperi)+m.log(rap))/2.)
          if self._R < Rmean:
              if self._R > rperi:
                  wR= (2.*m.pi/TR*
                       nu.array(integrate.quadrature(_TRAxiIntegrandSmall,
                                                     0.,m.sqrt(self._R-rperi),
                                                     args=(E,L,self._pot,rperi),
                                                     **kwargs)))\
                                                     +nu.array([m.pi,0.])
              else:
                  wR= nu.array([m.pi,0.])
          else:
              if self._R < rap:
                  wR= -(2.*m.pi/TR*
                        nu.array(integrate.quadrature(_TRAxiIntegrandLarge,
                                                      0.,m.sqrt(rap-self._R),
                                                      args=(E,L,self._pot,rap),
                                                      **kwargs)))
              else:
                  wR= nu.array([0.,0.])
          if self._vR < 0.:
              wR[0]+= m.pi
          self._angleR= nu.array([wR[0] % (2.*m.pi),wR[1]])
          return self._angleR
docs: NAME:
             AngleR
          PURPOSE:
             Calculate the radial angle
          INPUT:
             scipy.integrate.quadrature keywords
          OUTPUT:
             w_R(R,vT,vT) in radians + 
             estimate of the error (does not include TR error)
          HISTORY:
             2010-12-01 - Written - Bovy (NYU)