ПРЯМІ МЕТОДИ РОЗВ'ЯЗУВАННЯ СИСТЕМИ ЛІНІЙНИХ АЛГЕБРАЇЧНИХ РІВНЯНЬ

[dizzy69]

Здравствуйте!
Для виконання роботи необхiдно записати систему рiвнянь у матричнiй формi i скласти програму її розв'язування одним методом оптимального виключення

Вивести вектор розв'язку та нев'язки. Вектор нев'язки обчислюють за початковою системою.
Ось моя матриця
8 10 14 1 = 1
1000 0 2 5 = 2
-1 500 29 13=34
900 5 6 47 =6


це все що я зумів зробити
# include <stdlib.h>
# include <conio.h>
# include <stdio.h>
#include <conio.h>
void main()
{
clrscr();
const n=4;
int i=0,j=0,l=0;
float z=0;
float a[n][n+1]={{8,10,14,1,1},{1000,0,2,5,2},{-1,500,29,13,34},{900,5,6,47,6}};
for (i=0;i<n;i++) {
for(j=0;j<n+1;j++)
printf("%i | ",a[i][j]);
printf("\n");
}
for(i=0;i<n;i++)
{
z=a[i][i];
for(j=i;j<n+1;j++)
a[i][j]=a[i][j]/z;//if(a[i][j]==z)a[i][j]=0;}
for(l=i+1;l<n;l++)
{
z=a[l][i];
for(j=i+1;j<n+1;j++) a[l][j]=a[l][j]-a[i][j]*z;
a[l][i]=0;
}
for(l=i-1;l>=1;l--)
{
z=a[l][i];
for(j=i+1;j<n+1;j++) a[l][j]=a[l][j]-a[i][j]*z;
a[l][i]=0;
}
}
for (i=0;i<n;i++) {
printf("\n");
for(j=0;j<n+1;j++)
printf("%f | ",a[i][j]);
printf("\n");
getch();
}

}

що робити дальше незнаю