Esope source files

   1 : C QUADRA    SOURCE    CHAT      05/01/13    02:40:43     5004           
   2 :       SUBROUTINE QUADRA (CC2, CC1, CC0, X1, X2, NRoot)
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8(A-H,O-Z)
   5 : 
   6 :       DATA EPSILO /1.D-12/
   7 : 
   8 :       X1 = 0.d0
   9 :       X2 = 0.d0
  10 :       A = CC1 * CC1 - 4.d0 * CC2 * CC0
  11 : 
  12 :       IF(ABS(A).LT.EPSILO) THEN
  13 :          B = 2.d0 * CC2
  14 :          X1 = -CC1 / B
  15 :          X2 = X1
  16 :          NRoot = 2
  17 :          RETURN
  18 :       ENDIF
  19 : 
  20 :       IF( A.GT.0.)  THEN
  21 :          B = 2.d0 * CC2
  22 :          C = SQRT(A)
  23 :          X1 = (-CC1 + C) / B
  24 :          X2 = (-CC1 - C) / B
  25 :          NRoot = 2
  26 :       ELSE
  27 :          NRoot = 0
  28 :       ENDIF
  29 :       RETURN
  30 : 
  31 :       END
  32 :  
  33 :