Liviua19Ex2Krylov

#include <bits/stdc++.h>
using namespace std;

#define N 3

int main() {
    // A se seta N-ul ca si constanta la inceputul prgramului

    // Am scris matricea din exercitiu intr-o lista de liste denumita mat
    int mat[N][N] = {{1, -3, 1},
                     {2, 3, -1},
                     {1, -1, -1}};

    int y0[3] = {1, 0, 0}, y1[3] = {0, 0, 0}, y2[3] = {0, 0, 0}, y3[3] = {0, 0, 0};

    // calculam y1
    for (int i = 0; i < N; i++)
        for (int j = 0; j < N; j++)
            y1[i] += mat[i][j] * y0[j];

    // calculam y2
    for (int i = 0; i < N; i++)
        for (int j = 0; j < N; j++)
            y2[i] += mat[i][j] * y1[j];
    cout << endl;

    // calculam y3
    for (int i = 0; i < N; i++)
        for (int j = 0; j < N; j++)
            y3[i] += mat[i][j] * y2[j];
    cout << endl;

    int temp[N][N];

    // y0 ales de noi
    for (int i = 0; i < N; i++) {
        cout << y0[i] << " ";
        temp[i][0] = y0[i];
    }
    cout << endl;

    // y1 =
    for (int i = 0; i < N; i++) {
        cout << y1[i] << " ";
        temp[i][1] = y1[i];
    }
    cout << endl;
    // y2 =
    for (int i = 0; i < N; i++) {
        cout << y2[i] << " ";
        temp[i][2] = y2[i];
    }
    cout << endl;

    cout << "mat^(-1):" << endl;
    for (int i = 0; i < N; i++) {
        for (int j = 0; j < N; j++) {
            cout << temp[i][j] << " ";
        }
        cout << " |  "
             << "p" << i << "  | =" << y3[i];
        cout << endl;
    }

    return 0;
}