Курсовая работа: Линейные электрические цепи 2

// Прямой ход метода Гаусса: приводим матрицу к диагональному виду.

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

{

// Поиск максимального по модулю элементов в k-м столбце.

max=fabs (a[k] [k]);

r=k;

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

if (fabs (a[i] [k])>max)

{

max=fabs (a[i] [k]);

r=i;

}

// Меняем строки местами k-ю и r-ю.

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

{

c=a[k] [j];

a[k] [j]=a[r] [j];

a[r] [j]=c;

}

c=b[k];

b[k]=b[r];

b[r]=c;

// Приведение матрицы к диагональному виду.

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

{

for (M=a[i] [k]/a[k] [k], j=k; j<n; j++)

a[i] [j] – =M*a[k] [j];

b[i] – =M*b[k];

}