code: def _FullEOM(y,t,pot):
      
      l2= (y[0]**2.*y[3])**2.
      return [y[1],
              l2/y[0]**3.+_evaluateRforces(pot,y[0],y[4],phi=y[2],t=t,
                                           v=[y[1],y[0]*y[3],y[5]]),
              y[3],
              1./y[0]**2.*(_evaluatephiforces(pot,y[0],y[4],phi=y[2],t=t,
                                              v=[y[1],y[0]*y[3],y[5]])
                           -2.*y[0]*y[1]*y[3]),
              y[5],
              _evaluatezforces(pot,y[0],y[4],phi=y[2],t=t,
                               v=[y[1],y[0]*y[3],y[5]])]
docs: NAME:
         _FullEOM
      PURPOSE:
         implements the EOM, i.e., the right-hand side of the differential 
         equation
      INPUT:
         y - current phase-space position
         t - current time
         pot - (list of) Potential instance(s)
      OUTPUT:
         dy/dt
      HISTORY:
         2010-04-16 - Written - Bovy (NYU)