Esope source files

   1 : C POLYNO    SOURCE    CHAT      05/01/13    02:17:21     5004           
   2 :       FUNCTION POLYNO (COEFP,LCOEFP,ICOEFP,VAR)                                 
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8 (A-H,O-Z)                                                 
   5 : ************************************************************************        
   6 : *                                                                               
   7 : *                             P O L Y N O                                       
   8 : *                             -----------                                       
   9 : *                                                                               
  10 : * FONCTION:                                                                     
  11 : * ---------                                                                     
  12 : *                                                                               
  13 : *     DONNER LA VALEUR D'UN POLYNOME POUR UN NOMBRE DONNE.                      
  14 : *                                                                               
  15 : * PARAMETRES:   (E)=ENTREE   (S)=SORTIE   (+ = CONTENU DANS UN COMMUN)          
  16 : * -----------                                                                   
  17 : *                                                                               
  18 : *     COEFP   (E)  COEFFICIENTS DU POLYNOME (EN COMMENCANT PAR LA               
  19 : *                  PUISSANCE "0").                                              
  20 : *     LCOEFP  (E)  NOMBRE DE COEFFICIENTS UTILES DU TABLEAU "COEFP".            
  21 : *     ICOEFP  (E)  INCREMENT ENTRE 2 COEFFICIENTS UTILES DU TABLEAU             
  22 : *                  "COEFP". LE 1ER COEFFICIENT SE TROUVE EN POSITION 1.         
  23 : *     VAR     (E)  NOMBRE EN LEQUEL ON CHERCHE LA VALEUR DU POLYNOME.           
  24 : *     POLYNO  (S)  VALEUR DU POLYNOME.                                          
  25 : *                  = 0  SI LE TABLEAU "COEFP" EST VIDE.                         
  26 : *                                                                               
  27 :       INTEGER LCOEFP,ICOEFP                                                     
  28 :       DIMENSION COEFP(*)                                                        
  29 : *                                                                               
  30 : * AUTEUR, DATE DE CREATION:                                                     
  31 : * -------------------------                                                     
  32 : *                                                                               
  33 : *     PASCAL MANIGOT     31 JUILLET 1986                                        
  34 : *     P.M. 19/02/87: AJOUT DE LA VARIABLE "ICOEFP".                             
  35 : *                                                                               
  36 : * LANGAGE:                                                                      
  37 : * --------                                                                      
  38 : *                                                                               
  39 : *     FORTRAN77 + EXTENSION: DECLARATION "REAL*8".                              
  40 : *                                                                               
  41 : ************************************************************************        
  42 : *                                                                               
  43 :       ZZ = 0.D0                                                                 
  44 :       INIT = (LCOEFP-1) * ICOEFP + 1                                            
  45 :       IPAS = - ICOEFP                                                           
  46 :       DO 100 IB=INIT,1,IPAS                                                     
  47 :          ZZ = ZZ * VAR + COEFP(IB)                                              
  48 :   100    CONTINUE                                                               
  49 : *     END DO                                                                    
  50 :       POLYNO = ZZ                                                               
  51 : *                                                                               
  52 :       END                                                                       
  53 :