Решение матриц Методом Гаусса

#include "stdafx.h"
#include <iostream>
#include <stdio.h>
 
using namespace std;
 
void print(double *a);
void release(double *a);
 
int _tmain(int argc, _TCHAR* argv[])
{
    double a[3][4] = { { 0.42, -1.13, 7.05, 6.15 },
    { 1.14, -2.15, 5.11, -4.16 },
    { -0.71, 0.81, -0.02, -0.17 } };
    print(&a[0][0]);
 
    for (int i = 0; i < 3; i++)
    {
        for (int j = i + 1; j < 3; j++)
        {
            double r = -(a[j][i] / a[i][i]);
            for (int k = i; k < 4; k++)
            {
                a[j][k] += r * a[i][k];
            }
        }
        print(&a[0][0]);
    }
 
    double res[4];
    for (int i = 0; i < 3; i++)
        res[i] = 0;
 
    for (int i = 2; i >= 0; i--)
    {
        double r = 0.0;
        for (int j = i + 1; j <= 2; j++)
        {
            r += a[i][j] * res[j];
        }
        res[i] = (a[i][3] - r) / a[i][i];
    } 
 
    for (int i = 0; i < 3; i++)
    {
        cout << "x" << i + 1 << " = ";
        printf("%.4f\n", res[i]);
    }
 
    getchar();
    return 0;
}
 
 
void print(double *a)
{
    cout.precision(4);
    for (int i = 0; i < 3; i++)
    {
        for (int j = 0; j < 4; j++)
        {
            cout << *(a + i * 4 + j);
            (j < 2) ? cout << " + " : (j == 2) ? cout << " = " : cout << " ";
        }
        cout << endl;
    }
    cout << endl;
}