gauss_seidel.rar_To the Letter

#include <stdio.h> #include <math.h> #include <conio.h> /*fonction de lecture de la matrice a et du vecteur b*/ void lire(float a[19][19],float b[19],float x0[19],int &n) { int i,j; for(i=0;i<n;i++) for(j=0;j<n;j++) {printf("a[%d][%d]",i,j); scanf("%f",&a[i][j]);} for(j=0;j<n;j++) { printf("donnez l element b[%d]= ",j); scanf("%f",&b[j]);} for(j=0;j<n;j++) { printf("donnez l elements x0[%d]= \n", j); scanf("%f",&x0[j]); } } /*fonction d affichage*/ void affiche_systeme(float a[19][19],float b[19],float x0[19],int& n){ int i,j; for (i=0;i<n;i++){ for(j=0;j<n;j++) printf("%.10f ",a[i][j]); printf("%.10f",b[i]); printf("\n");} printf("\n"); } /*fonction de la decomposition lu*/ void iteration_jacobie(float a[19][19],float b[19],float x[19],float x0[19],int& n){ int i,j,k; float s; for (i=0;i<n;i++){ s=0; for (j=0;j<i;j++) s=s+a[i][j]*x[j]; for (j=i+1;j<n;j++) s=s+a[i][j]*x0[j]; x[i]=(1/a[i][i])*(b[i]-s);} for (i=0;i<n;i++) x0[i]=x[i]; } /*fonction d affichage de la solution*/ void afficher_solution(float x[19],int &n){ int i; for(i=0;i<n;i++) printf("%.10f\n",x[i]); } /***************************************************************************/ /* debut de la fonction main du programme*/ /********************************************************************************/ main(){ float a[19][19], b[19],x[19],x0[19], s; int n,i,iteration; /*lecture de la dimension n*/ printf("donner la dimension de votre matrice et de votre vecteur \n\n"); scanf("%d",&n); /*lecture de nombre d iteration de jacobie*/ printf("donner le nombre d iteration \n\n"); scanf("%d",&iteration); /*lecture des donnees*/ lire(a,b,x0,n); printf("\n"); /*affichage des elements systeme*/ printf("Voici les elemnts du systeme \n\n"); affiche_systeme(a,b,x0,n); /*decomposition lu*/ for(i=1;i<=iteration;i++){ iteration_jacobie(a,b,x,x0,n); printf("Voici les elemnts de la solution d iteration %d \n\n",i); afficher_solution(x,n); printf("Appyer sur une touche pour voir la solution d iteration suivant\n\n"); getch(); } getch(); } /********************************************************************************/ /* fin de la fonction main du programme*/ /********************************************************************************/