Télécharger inte22.eso

Retour à la liste

Numérotation des lignes :

  1. C INTE22 SOURCE CHAT 05/01/13 00:39:53 5004
  2. SUBROUTINE INTE22(PRT,PRF,T0,FT0,LON)
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8(A-H,O-Z)
  5. C
  6. C ====================================================================
  7. C = INTERPOLATION LINEAIRE RECHERCHE PAR DICHOTOMIE =
  8. C = CREATION 14/01/86 =
  9. C = PROGRAMMEUR GUILBAUD =
  10. C ====================================================================
  11. C
  12. -INC CCOPTIO
  13. DIMENSION PRT(1),PRF(2)
  14. DATA EPSI/1.D-10/
  15. N1=1
  16. N2=LON
  17. 1 CONTINUE
  18. N2N1=N2-N1
  19. IF(N2N1.EQ.1) GOTO 2
  20. NC=(N1+N2)/2
  21. IF(T0.LT.PRT(NC)) THEN
  22. N2=NC
  23. ELSE
  24. N1=NC
  25. ENDIF
  26. GOTO 1
  27. 2 CONTINUE
  28. ALPHA=(T0-PRT(N1))/(PRT(N2)-PRT(N1))
  29. FT0=ALPHA*(PRF(N2)-PRF(N1))+PRF(N1)
  30. BETA=1.D0-ALPHA
  31. IF(BETA.LT.EPSI)FT0=PRF(N2)
  32. RETURN
  33. END
  34.  
  35.  

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