Télécharger quadra.eso

Retour à la liste

Numérotation des lignes :

  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.  
  34.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales