salut tous le monde
A quoi sert la subroutine tri dag (j'ai le code en fortran):
merci d'avance
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
Partager