რიცხვითის ლაბ. N1

/*
"data.txt"
----------
1 0 4
0 1 0
0 0 3

2 -7 3
----------
*/

#include <iostream>
#include <fstream>
#include <vector>
#include <cmath>

using namespace std;

const int n = 3;

int main()
{
    //A მატრიცისა და x ვექოტრის შექმნა:
    vector<vector<double>> A(n, n);
    vector<double> x(n);
    //

    //ფაილიდან წაკითხვა:
    ifstream ifs("data.txt");

    for (int i = 0; i < n; i++)
        for (int j = 0; j < n; j++)
            ifs >> A[i][j];

    for (int i = 0; i < n; i++)
        ifs >> x[i];

    ifs.close();
    //

    //Ax ვექტორის დათვლა:
    vector<double> Ax(n);

    for (int i = 0; i < n; i++)
        for (int j = 0; j < n; j++)
            Ax[i] += A[i][j] * x[j];

    //ვექტორის A ნორმა:
    double Norm(0);

    for (int i = 0; i < n; i++)
        Norm += Ax[i] * x[i];

    cout << "x veqtoris A norma = " << sqrt(Norm) << endl;
    //

    //მატრიცის პირველი ნორმა:
    Norm = 0;
    double Sum(0);
    for (int i = 0; i < n; i++)
    {
        Sum = 0;

        for (int j = 0; j < n; j++)
            Sum += fabs(A[j][i]);

        if (Sum > Norm)
            Norm = Sum;
    }

    cout << "A matricis pirveli norma = " << Norm << endl;
    //

    //მატრიცის უსასრულობა ნორმა:
    Norm = 0;
    Sum = 0;
    for (int i = 0; i < n; i++)
    {
        Sum = 0;

        for (int j = 0; j < n; j++)
            Sum += fabs(A[i][j]);

        if (Sum > Norm)
            Norm = Sum;
    }

    cout << "A matricis usasruloba norma = " << Norm << endl;
    //

    system("PAUSE");
    return (0);
}