salut
je dois coder l'algo de gauss-seidel sur une ti-92 plus...
en paramètre il doit avoir (a, b, x0,k)
a est une matrice carré, b un vecteur colonne, x0 une approximation inixiale (vecteur colonne), k un entier k>=1
je doit retourner le résultat de la kieme itération de cette méthode...
si me rappele bien, la formule de gauss-seidel est
X^(k+1) = (D-L)^(-1) U * X^K + (D-L)^-1 *B
c'est de la difficulté a associé les paramètre que je dois utiliser à ceux de la formule...
merci
Je vais regarder ton problème. Mais va plutôt sur le forum de:
http://tigen.org
je joins l'algo que j'ai vue sur:
http://www.unige.ch/ses/metri/gilli/.../MNE-Cours.pdf
je comprends pas comment décomposé cela
je viens de trouver:
http://phys.ucc.ie/~oll/AM2054/programs/Jacobi and gauss-Seidel.htm
la fonction en c++
au lieu de faire un do while pour savoir si on approche le taux d'erreur voulu, je ne fais que réitérer tant que la valeur n'est pas = à k....
ça donne
à l'exécutionCode:gasei(a,b,x0,k) Prgm if coldim(a) != rowdim(a) then disp "matrice non carree" return 0 elseif coldim(b)!=1 then disp "nombre de colonne !=1 return 0 elseif coldim(x0)!=1 then disp "nombre de colonne !=1 return 0 elseif rowdim(a) != rowDim(b) then disp "le vecteur colonne b n'a pas le même ordre que a" return 0 elseif rowdim(a) != rowDim(x0) then disp "le vecteur colonne x0 n'a pas le même ordre que a" return 0 elseif k<1 then disp "k est inferieur a 1" return 0 else local count, tmp for count,0, k for i,0,rowdim(a) tmp=b[i] for j,0,i tmp=tmp-a[i][j]*x0[j] endfor x0[i]=tmp/a[i][i] endfor endfor for i,0, rowdim(x0) disp x0[i] endfor endif endprgm
[10,0,0;1,15,0;-1,1,20] -> a
[2;12;17] -> b
[0;0;0] -> x0
1 -> k
gasei(a,b,x0,k)
il me retourne dimension et va à la ligne
tmp=b[i]
une idée?
