Esope source files

   1 : C DEG3      SOURCE    CHAT      05/01/12    22:40:55     5004           
   2 :       SUBROUTINE DEG3
   3 : C
   4 : C     CALCUL DES RACINES D'UN POLYNOME DU 3-EME DEGRE
   5 : C
   6 :       IMPLICIT INTEGER(I-N)
   7 :       IMPLICIT REAL*8 (A-H,O-Z)
   8 : -INC CCOPTIO                                                                    
   9 : C                                                                               
  10 : C     LECTURE DES COEFFICIENTS DU POLYNOME
  11 : C                                                                               
  11 :       CALL LIRREE(A0,1,IRETOU)
  12 :       IF(IERR.NE.0) RETURN
  12 :       CALL LIRREE(A1,1,IRETOU)
  13 :       IF(IERR.NE.0) RETURN
  13 :       CALL LIRREE(A2,1,IRETOU)
  14 :       IF(IERR.NE.0) RETURN
  14 :       CALL LIRREE(A3,1,IRETOU)
  15 :       IF(IERR.NE.0) RETURN
  16 : *
  17 : *  MISE SOUS FORME POUR DEGRE3
  18 : *
  19 :       IF (A3.EQ.0.D0) THEN
  19 :            CALL ERREUR(616)
  20 :            RETURN
  21 :       ENDIF
  22 :       AS0=A0/A3
  23 :       AS1=A1/A3
  24 :       AS2=A2/A3
  24 :       CALL DEGRE3(AS0,AS1,AS2,XR1,XI1,XR2,XI2,XR3,XI3)
  24 :       CALL ECRREE(XI3)
  24 :       CALL ECRREE(XR3)
  24 :       CALL ECRREE(XI2)
  24 :       CALL ECRREE(XR2)
  24 :       CALL ECRREE(XI1)
  24 :       CALL ECRREE(XR1)
  25 :       RETURN                                                                    
  26 :       END
  27 :