Wczytywanie danych z pliku.

#include <cstdlib>
#include <iostream>
#include <fstream>
#include <vector>
#include <math.h>

using namespace std;

/*
struct punkt
{
    double x;
    double y;
};


*/

/* cout << "Wczytywanie danych x, y" << endl << endl;
    fstream plik;
    plik.open("X:\\Mikołaj Pilch Projekt\\dane.txt");
    if (plik.good())
    {
        vector<punkt> wspolrzedne;
        while (!plik.eof())
            {
                getline(plik, );
                n++;
            }
    }
    */

//plik.close();

int main( int argc, char * argv[] )
{
   
    int i;
    int n;
   
    double sum_x = 0,
    sum_y = 0,
    sum_xy = 0,
    sum_x2 = 0,
    a,
    b,
    ub,
    ua,
    s,
    pierwua,
    pierwub,
    sum_y2 = 0,
    pierws,
    sumxs2 = 0,
    sumys2 = 0,
    sredniax = 0,
    sumax = 0,
    sumay = 0,
    sredniay = 0,
    licznik = 0,
    r;
   
   
   
   
    cout << "Podaj ilosc punktow: ";
   
    cin >> n;
    cout << endl << endl;
    double x[ n ];
    double y[ n ];
   
    for( i = 0; i < n; i++ )
    {
        cout << "Podaj x: ";
        cin >> x[ i ];
        cin.ignore();
        cout << endl;
        cout << "Podaj y: ";
        cin >> y[ i ];
        cin.ignore();
        cout << endl;
    }
    cout << "Aproksymacja liniowa do postaci y=ax+b " << endl;
   
    for( i = 0; i < n; i++ )
    {
        sumax = sumax + x[ i ];
        sumay = sumay + y[ i ];
    }
   
    sredniax = sumax / n;
    sredniay = sumay / n;
   
   
    for( i = 0; i < n; i++ )
    {
        sum_x = sum_x + x[ i ];
        sum_y = sum_y + y[ i ];
        sum_xy = sum_xy + x[ i ] * y[ i ];
        sum_x2 = sum_x2 + x[ i ] * x[ i ];
        sum_y2 = sum_y2 + y[ i ] * y[ i ];
        sumxs2 = sumxs2 +( x[ i ] - sredniax ) *( x[ i ] - sredniax );
        sumys2 = sumys2 +( y[ i ] - sredniay ) *( y[ i ] - sredniay );
        licznik = licznik +( x[ i ] - sredniax ) *( y[ i ] - sredniay );
       
    }
   
    a =( n * sum_xy - sum_x * sum_y ) /( n * sum_x2 - sum_x * sum_x );
   
    b =( sum_y - a * sum_x ) / n;
   
    ua =( n *( sum_y2 - a * sum_xy - b * sum_y ) ) /(( n - 2 ) *( n * sum_x2 - sum_x * sum_x ) );
   
    pierwua = sqrt( ua );
    ub =( pierwua * pierwua * sum_x2 ) / n;
    pierwub = sqrt( ub );
   
   
   
    cout << "Rownanie lini prostej wynosi y=" << a << "x+" << b << endl;
   
    cout << "Blad u(a) wynosi: " << pierwua << endl;
    cout << "Blad u(b) wynosi: " << pierwub << endl;
   
   
   
    s = sumxs2 * sumys2;
    pierws = sqrt( s );
    r = licznik / pierws;
    cout << "Wspołczynnik korelacji Pearsona wynosi: " << r << endl;
   
   
   
   
   
    cin.get();
    return EXIT_SUCCESS;
}