#include<iostream>
#include<cmath>// для сравнения
#include<Windows.h>
#include<ctime>
#include<random>

using namespace std;


int main()
{
    random_device rd;
    mt19937 gen(rd());
    double* results = new double[26];
    double summ = 0;
    srand(time(NULL));
    SetConsoleCP(1251);
    SetConsoleOutputCP(1251);


    double x;
    uniform_real_distribution<>dis(-1, 1);
    x = dis(gen);           //x = (rand() % 201 - 100.0)/ 100.0;задаём икс в диапазоне [-1;1]
    cout << " x= " << x << endl;
    double m;
    m = sqrt(x + 1);

    cout << "Результат стандартной функции равен\t" << m << endl;

    cout << "Результат моей функции равен\t";
    results[0] = 1;
    results[1] = 0.5 * x;
    const int MAX_COUNT = 25;
    for (int counter = 2; counter <= MAX_COUNT; counter++)
    {
        double n = counter;
        results[counter] = results[counter - 1] * ((2 * n) - 3) * x * (-1) / (2 * n);//даже не нужен факториал и степень
    }
    for (int counter = 0; counter <= MAX_COUNT; counter++)
        summ += results[counter];
    cout << summ << endl;
    cout << "Погрешность измерений равна: ";
    double eps;
    eps = results[25];
    cout << abs(eps) << endl;
    //result += double_fact(2 * counter - 3) * stepen(x, counter) * stepen(-1, counter-1) / double_fact(2 * counter);-на всякий случай оставлю
    delete[]results;
    system("pause");
    return 0;


}