double sum=0;
        for (int j = 0; j < n; j++)
        {
            sum += powf((function_values[j].y - aprox_values[j].y),2);
                //function_values[j].y = wynik prawidlowy; aprox_values[j].y = aproksymowany
        }
        sum = abs(sum);
        cout << "Blad aproksymacji = " << sum << endl;