#include <iostream>
#include <fstream>
#include <time.h>

using namespace std;


int main()

{
    const int n=5;

    double d, s;

    double a[n][n]= {{1,3,0,0,0},
                     {1,2,1,0,0},
                     {0,1,2,1,0},
                     {0,0,1,2,1},
                     {0,0,0,1,1}};
    double a1[n][n];
    double b[n]={1,5,3,2,2};
    double x[n];

    cout << "Macierz: " << endl;
    for (int i = 0; i < n; i++) {
        for (int j = 0; j < n; j++) {
            cout << a[i][j] << "\t";
        }
        cout << b[i];
        cout << endl;
        }
    clock_t t;
    t=clock();
    for (int k = 0; k < n; k++) {
        for (int j = k + 1; j < n; j++) {
            d = a[j][k] / a[k][k]; //mno¿enie wiersza przez liczbe
            for (int i = k; i < n; i++) {
                a[j][i] = a[j][i] - d * a[k][i]; // odejmowanie k-ego wirsza z j-ego
            }
            b[j] = b[j] - d * b[k]; // 2 czesc odejmowania
        }
    }
    for (int k=0; k<n; k++) // poszukiwanie niewiadomych
        x[k]=0;
    for (int k = n-1; k >= 0; k--){
        d = 0;
        for (int j = k ; j < n; j++){
            s = a[k][j] * x[j];
            d = d + s;
        }
            x[k] = (b[k] - d) / a[k][k];
    }
    cout << "Niewiadome: " << endl;
    for(int i = 0; i < n; i++)
        cout << "x[" << i << "]=" << x[i] << " " << endl;
    cout << endl;

    cout << "Uruchomienie funkcji 100000001 (raz z wynikami) raz:\n ";
    for (int st=0; st<10000000; st++) {

    for (int k = 0; k < n; k++) {
        for (int j = k + 1; j < n; j++) {
            d = a[j][k] / a[k][k]; //mno¿enie wiersza przez liczbe
            for (int i = k; i < n; i++) {
                a[j][i] = a[j][i] - d * a[k][i]; // odejmowanie k-ego wirsza z j-ego
            }
            b[j] = b[j] - d * b[k]; // 2 czesc odejmowania
        }
    }
    for (int k=0; k<n; k++) // poszukiwanie niewiadomych
        x[k]=0;
    for (int k = n-1; k >= 0; k--){
        d = 0;
        for (int j = k ; j < n; j++){
            s = a[k][j] * x[j];
            d = d + s;
        }
            x[k] = (b[k] - d) / a[k][k];
    }
    }
   // gauss();
    t=clock()-t;
    cout<<"Clicks: " << t <<". Czas potrzebny dla wywolania: "<<((float)t)/CLOCKS_PER_SEC;
    return 0;
}