Вот что уменя получилось:
Код:
#include "stdafx.h"
#include <fstream>
using namespace std;
int _tmain(int argc, _TCHAR* argv[])
{
int n = 4;
float a[4][4];
ifstream F("1.txt");
for (int i = 0; i < n * n; ++i)
F >> a[i / n][i % n];
F.close();
for (int i = 0; i < n; ++i){
for (int j = 0; j < n; ++j)
printf("%.f ", a[i][j]);
printf("\n");
}
int i, j;
for (int k = 0; k < n; ++k){
a[k][k] = sqrt(a[k][k]);
for (i = k + 1; i < n; ++i){
a[i][k] = a[i][k] / a[k][k];
}
printf("%d\n", k);
for (j = k + 1; j < n; ++j){
for (i = j; i < n; ++i)
a[i][j] = a[i][j] - a[i][k] * a[j][k];
// for (i = 0; i < k; ++i)
// a[j][k + 1] = a[j][k + 1] - a[i][j] * a[k + 1][i];
}
}
for (i = 0; i < n; ++i){
for (j = 0; j < n; ++j)
printf("%.f ", a[i][j]);
printf("\n");
}
return 0;
}