#include <iostream>
#include <fstream>

using namespace std;
//float Aitken(float wsspol[1][0], float wsspol[0][1],float wsspol[2][0], float wsspol[0][2], float wsspol[3][0],
//              float wsspol[0][3], float wsspol[4][0], float wsspol[0][4], float argument) //zadeklarowanie funkcji i zmiennych
// (xi,yi) - wspolrzedne pierwszego punktu, (xj,yj) - wspolrzedne drugiego punktu, (xk,yk) - wspolrzedne trzeciego punktu
float Aitken(float tablica[5][5], float argument)
{
    /*float w01,w02,w012,w03,w013,w0123; //zadeklarowanie zmiennych wielomianów
    w01=(yi*(xj-argument)-(yj*(xi-argument)))/(xj-xi); // wielomian pierwszego stopnia Wij
    w02=(yi*(xk-argument)-(yk*(xi-argument)))/(xk-xi); // wielomian pierwszego stopnia Wik
     w03=(yi*(xl-argument)-(yl*(xi-argument)))/(xl-xi); // wielomian pierwszego stopnia Wil
    w012=(w01*(xk-argument)-(w02*(xj-argument)))/(xk-xj); // wielomian drugiego stopnia Wijk
    w013=(w01*(xl-argument)-(w03*(xj-argument)))/(xl-xj); // wielomian drugiego stopnia Wijl
    w0123=(w012*(xl-argument)-(w013*(xk-argument)))/(xl-xk); // wielomian trzeciego stopnia Wijkl*/
    float w[4][4];
    for(int i=1;i<=3;i++)
    {
        w[0][i]=(tablica[0][i]*(tablica[i+1][0]-argument)-(tablica[0][i+1]*(tablica[i][0]-argument)))/(tablica[i+1]-tablica[i]);
    }


 return w[0][1];

}

int main()
{
   float wspol[5][5],argument; // zadeklarowanie zmiennych do których przypiszemy wartosci z pliku
   ifstream plik; // otwarcie strumienia umozliwiajacego obsluge wejscia plikowego
   plik.open("dane.txt"); // otworzenie pliku
   if(plik.good()==false) //sprawdzenie czy plik istnieje
   {
       cout<<" plik nie istnieje "<<endl;
   }

   else{

plik>>argument>>wspol[1][0]>>wspol[0][1]>>wspol[2][0]>>wspol[0][2]>>wspol[3][0]>>wspol[0][3]>>wspol[4][0]>>wspol[0][4];  // wprowadzenie danych z pliku do zadeklarowanych wczesniej zmiennych
cout<<"wartosc wielomianu interpolacyjnego trzeciego stopnia przechodzacego przez punkty ";
cout<<'('<<wspol[1][0]<<','<<wspol[0][2]<<')'<<'('<<wspol[2][0]<<','<<wspol[0][2]<<')'<<'('<<wspol[3][0]<<','<<wspol[0][3]<<')'<<'('<<wspol[4][0]<<','<<wspol[0][4]<<')'<<" dla argumentu x = "<<argument<<" wynosi ";

 cout<<Aitken(wspol[5][5],argument);}

    return 0;
}