Метод Холесского (Краута) для решения систем линейных уравнений (скрин с алгоритмом вложен))

[Bogdand]

#include <stdio.h>
#include <conio.h>
#include <stdlib.h>
main ()
{

FILE *f;

float x[4],a[4][5];
unsigned i,j,
m,
n = 4,k;
char p;

f = fopen("sist.txt","rt");

printf ("massiv \n");
for (i = 0; i < 4; i ++)
for (j = 0; j < 5; j++)
{
fscanf(f,"%f", &a[i][j]);
fscanf(f,"%c", &p);
}


for (i = 0; i < 4; i++)
{
printf("\n");
for(j = 0; j < 5; j++)
printf("%g ",a[i][j]);
}
// Для первой строки матрицы
for (j = 2; j <= n+1; j++)
a[1][j] = a[1][j]/ a[1][1];

// Для очередного l - ого столбца
m = 2;
while (m <= n)
{
for (i = m; i <= n; i++)
for (k = 1; k < m; k++)
a[i][m] = a[i][m] - a[i][k] * a[k][m];

m++;
}

// Для очередной U - ой строки
m = 2;
while (m <= n)
{
for (j = m+1; j <= n + 1; j++)
for (k = 1; k < m; k++)
a[m][j] = (a[m][j] - a[k][j]*a[m][k]) / a[m][m];

m++;
}
// Подстановка Х
for (i = n-1; i >= 1; i--)
for (k = i + 1; k <= n; k++)
x[i] = a[i][n+1] - a[i][k]*x[k];


printf("\n");
printf("Otvet \n");
for (i = n - 1; i >= 1; i--)
printf(" x = %3.2f",x[i]);

getch();


}

Выводит абы что