salut tous le monde

A quoi sert la subroutine tri dag (j'ai le code en fortran):

Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
 
subroutine tridag(if,l,a,b,c,d,rr)
      implicit real*8(A-h,o-z)
      dimension A(1),b(1),c(1),d(1),rr(1),beta(20000),gamma(20000)
      beta(if)=b(if)
      gamma(if)=d(if)/beta(if)
      ifp1=if+1
      do 1 i=ifp1,l
      beta(i)=b(i)-a(i)*c(i-1)/beta(i-1)
      gamma(i)=(d(i)-a(i)*gamma(i-1))/beta(i)
   1  continue
      rr(l)=gamma(l)
      last=l-if
      do 2 k=1,last
      i=l-k
      rr(i)=gamma(i)-c(i)*rr(i+1)/beta(i)
   2  continue
      return
      end
merci d'avance