Sources Esope | Cast3M

quadra
C QUADRA    SOURCE    CHAT      05/01/13    02:40:43     5004
      SUBROUTINE QUADRA (CC2, CC1, CC0, X1, X2, NRoot)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
 
      DATA EPSILO /1.D-12/
 
      X1 = 0.d0
      X2 = 0.d0
      A = CC1 * CC1 - 4.d0 * CC2 * CC0
 
      IF(ABS(A).LT.EPSILO) THEN
         B = 2.d0 * CC2
         X1 = -CC1 / B
         X2 = X1
         NRoot = 2
         RETURN
      ENDIF
 
      IF( A.GT.0.)  THEN
         B = 2.d0 * CC2
         C = SQRT(A)
         X1 = (-CC1 + C) / B
         X2 = (-CC1 - C) / B
         NRoot = 2
      ELSE
         NRoot = 0
      ENDIF
      RETURN
 
      END