stability.f

Go to the documentation of this file.

    real*8 function stability(mm1,mm2,mm3,ein,eout,inc)

*       three-body stability test(mardling & aarseth 1999)
*
        implicit real*8 (a-h,m,o-z)
        REAL*8 inc
    
*       employ the semi-analytical stability criterion(ma 1999).
      q = mm3/(mm1 + mm2)
      IF (eout.LT.1.0) THEN
          xfac = (1.0 + q)*(1.0 + eout)/sqrt(1.0 - eout)
      ELSE
          xfac = 40.0*(1.0 + q)
      END IF
      pcrit = 2.8*xfac**0.4
*
*       include the inclination fudge factor.
      yfac = 1.0 - 0.3*inc/3.1415
      pcrit = yfac*pcrit
*
        stability=pcrit

    end