Exemple A_z_mevisu > Fichier difdiv.f |
SUBROUTINE DIFDIV (T,V,NB)
* CALCUL DES DIFFERENCES DIVISEES
* Argument d'entree :
* T (real) : VALEURS DU PARAMETRE
* V (real) : VALEURS A INTERPOLER CORRESPONDANTES
* NB (integer) : NOMBRE DE POINTS D'INTERPOLATION
* Argument de sortie :
* V (real) : VALEURS DE DIFFERENCES DIVISEES
REAL T(*),V(*)
INTEGER NB,I,J
*
DO 1 I=1,NB-1
DO 2 J=0,NB-I-1
V(NB-J)=(V(NB-J)-V(NB-J-1))/(T(NB-J)-T(NB-J-I))
2 CONTINUE
1 CONTINUE
END