Курсовая работа: Поиск решений системы линейных уравнений методом Гаусса

double A[n][n]; // матрицакоэффицентовГаусса

double B[n]; // рабочий массив свободных членов

int i,

int j,

int k; // рабочие переменные

for( i = 0; i < n; i++ ) // копирование в рабочую матрицу A

{

for( j = 0; j < n; j++ )

A[i][j] = Z[i][j];

B[i] = Y[i]; // копирование свободных членов

X[i] = 0;

}

for( k = 0; k < n-1; k++ )

for( i = k+1; i < n; i++ ) // преобразованиестрок

{

double r = A[i][k] / A[k][k];

for( j = k; j < n; j++ )

A[i][j] -= A[k][j] * r;

B[i] -= B[k] * r;

}

X[n-1] = B[n-1] / A[n-1][n-1];

for( i = n-2; i >= 0; i-- )

{

double s = 0;

for( j = i+1; j < n; j++ )

s += A[i][j] * X[j];

X[i] = ( B[i] - s ) / A[i][i];

}

printf( "\n" );