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.  
  13. -INC PPARAM
  14. -INC CCOPTIO
  15. DIMENSION PRT(1),PRF(2)
  16. DATA EPSI/1.D-10/
  17. N1=1
  18. N2=LON
  19. 1 CONTINUE
  20. N2N1=N2-N1
  21. IF(N2N1.EQ.1) GOTO 2
  22. NC=(N1+N2)/2
  23. IF(T0.LT.PRT(NC)) THEN
  24. N2=NC
  25. ELSE
  26. N1=NC
  27. ENDIF
  28. GOTO 1
  29. 2 CONTINUE
  30. ALPHA=(T0-PRT(N1))/(PRT(N2)-PRT(N1))
  31. FT0=ALPHA*(PRF(N2)-PRF(N1))+PRF(N1)
  32. BETA=1.D0-ALPHA
  33. IF(BETA.LT.EPSI)FT0=PRF(N2)
  34. RETURN
  35. END
  36.  
  37.  

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