Решение СЛАУ методом Гауса ( немного подправить программу)

[pe3a4ok]

Всем привет!
Подскажите что сделать) Программа считает все матрицы, кроме тех у, которых на диагонали стоит 0.
Насколько я знаю, нужно сделать обмен строк матрицы...
но как притворить это в жизнь - пока с трудом догадываюсь
помогите пожалуйста!!!
Вот тут весь код на языке C:

# include <stdio.h>
# include <conio.h>
int main()
{
 int i, j, k, n ;
 float a[20][20], x[20] ;
 double s, p ;
 
 printf("Enter the size of the matrix : ") ;
 scanf("%d", &n) ;
 printf("\nEnter the co-efficients of the matrix :\n\n") ;
 for(i = 0 ; i < n ; i++)
 {
  for(j = 0 ; j < n ; j++)
  {
   printf("a[%d][%d] = ", i + 1, j + 1) ;
   scanf("%f", &a[i][j]) ;
  }
  printf("b[%d] = ", i + 1) ;
  scanf("%f", &a[i][n]) ;
 }
 for(k = 0 ; k < n - 1 ; k++)
 {
  for(i = k + 1 ; i < n ; i++)
  {
   p = a[i][k] / a[k][k] ;
   for(j = k ; j < n + 1 ; j++)
    a[i][j] = a[i][j] - p * a[k][j] ;
  }
 }
 x[n-1] = a[n-1][n] / a[n-1][n-1] ;
 for(i = n - 2 ; i >= 0 ; i--)
 {
  s = 0 ;
  for(j = i + 1 ; j < n ; j++)
  {
   s += (a[i][j] * x[j]) ;
   x[i] = (a[i][n] - s) / a[i][i] ;
  }
 }
 printf("\nThe result is :\n") ;
 for(i = 0 ; i < n ; i++)
  printf("\nx[%d] = %.2f", i + 1, x[i]) ;
 getch() ;
}

[SaLoKiN]

с прямым ходом напутано было у вас...

# include <stdio.h>

int main()
{
 int i, j, k, n ;
 float a[20][20], x[20] ;
 double s, p ;
 
 printf("Enter the size of the matrix : ") ;
 scanf("%d", &n) ;
 printf("\nEnter the co-efficients of the matrix :\n\n") ;
 for(i = 0 ; i < n ; i++)
 {
  for(j = 0 ; j < n ; j++)
  {
   printf("a[%d][%d] = ", i + 1, j + 1) ;
   scanf("%f", &a[i][j]) ;
  }
  printf("b[%d] = ", i + 1) ;
  scanf("%f", &a[i][n]) ;
 }
 for(k = 0 ; k < n - 1 ; k++)
 {
  for(i = k + 1 ; i < n ; i++)
  {
   p = a[k][k] / a[i][k] ;
   for(j = 0 ; j <= n ; j++)
   a[i][j] = a[i][j] *p - a[k][j] ;
   }
 }
 x[n-1] = a[n-1][n] / a[n-1][n-1] ;
 for(i = n - 2 ; i >= 0 ; i--)
 {
  s = 0 ;
  for(j = i + 1 ; j < n ; j++)
  {
   s += (a[i][j] * x[j]) ;
   x[i] = (a[i][n] - s) / a[i][i] ;
  }
 }
 printf("\nThe result is :\n") ;
 for(i = 0 ; i < n ; i++)
  printf("\nx[%d] = %.2f", i + 1, x[i]) ;
 getch() ;

}

что касается нулей, вот тут есть проверочка на нули.
кста, судя по реализации обратного метода, код был взят отсюда....