Esope source files

   1 : C FLA017    SOURCE    TTMF3     12/07/05    21:15:40     7425           
   2 :       SUBROUTINE FLA017(X,Y,YP,R,A,B,D,E)
   3 : C------------------------------------------------------------------------
   4 : C Calcul de f(x) et f'(x) 
   5 : C (voir fla014.eso)
   6 : C------------------------------------------------------------------------
   7 : C
   8 : C---------------------------
   9 : C Parametres Entree/Sortie :
  10 : C---------------------------
  11 : C
  12 : C E/  X      : flottant : Variable 
  13 : C  /S Y      : flottant : Valeur de la fonction f(X)
  14 : C  /S YP     : flottant : Valeur de f'(X)
  15 : C E/  R      : flottant : Coefficient multiplicateur
  16 : C E/  A      : flottant : Coefficient devant X^2
  17 : C E/  B      : flottant : Coefficient devant X
  18 : C E/  D      : flottant : Coefficient
  19 : C E/  E      : flottant : Coefficient 
  20 : C
  21 : C------------------------------------------------------------------------
  22 : C
  23 : C Langage : FORTRAN 77
  24 : C
  25 : C Mise en oeuvre : H. Paillère (1997, TTMF)
  26 : C
  27 : C------------------------------------------------------------------------
  28 :       IMPLICIT INTEGER(I-N)
  29 :       IMPLICIT REAL*8 (A-H,O-Z)
  30 : C
  31 :       X2  = X*X
  32 :       R2  = R*R
  33 :       E2  = E*E
  34 :       XD2 = (X+D)*(X+D)
  35 :       Y   = A*X2 + B*X - R2*(1-X2/(E2*XD2))
  36 :       YP  = 2.D0*A*X + B + 2*D*R2*X/(E2*(X+D)*XD2)
  37 : C
  38 :       RETURN 
  39 :       END
  40 : 
  41 :  
  42 :  
  43 :  
  44 :