Вычисление значений многочлена по схеме Горнера методом сдваивания

#include <windows.h>
#include <stdlib.h>     // включение генератора случайных чисел
#include <time.h>               // получение первого случайного числа
#include <iostream.h>           // ввод-вывод
#include <conio.h>              // для включения ожидания ввода клавиши
#include <dos.h>                // для включения sleep()
#include <math.h>
#define N 20                    // степень многочлена
double x;                       // значение переменной
unsigned long long factorial(unsigned long long n)   // для вычисления факториала
{
    if (n == 0)
        return 1;
    else
        return n * factorial (n - 1);
}
 
struct arg                      // структура для передачи параметров
{
int l, r;                       // нижний и верхний номера 2х2=матрицы
double p, q;                    // коэффициенты первой строки матрицы
};
 
double a[N+1];                      // коэффициенты многочлена
DWORD WINAPI prod(void* s)      // функция потока
{
int i, l=((arg *)s) -> l, r=((arg *)s) -> r;
Sleep(random(0.2));             // случайная задержка
if (l == r)
{
((arg *)s) -> p = x*x;      // установка
((arg *)s) ->q = a[r];      // коэффициентов
}
else
{
arg *r1 = new arg, *r2 = new arg;
HANDLE H[2];                                    // дескрипторы потоков
r1 -> l = l; r1 -> r = (l + r + 1)/2 - 1;       // установка
r2 -> l = (l + r + 1)/2; r2 -> r = r;          // границ
H[0] = CreateThread(0,0,prod, (void *)r1, 0,0); // порождение
H[1] = CreateThread(0,0,prod, (void *)r2,0,0);  // потоков
for (i = 0; i < 2; i++)                             // ожидание выполнения
WaitForSingleObject(H[i],INFINITE);             // потоков
((arg *)s) -> p = (r1 -> p)*(r2 -> p);
((arg *)s) -> q = (r1 -> p)*(r2 -> q)+(r1 -> q);
delete r1; delete r2;
}
return 1;
}
 
void main()
{
int i,g,j;
double w=0, z, v,b;
for (i=0; i<=N; i++)
{
 
a[i] = 1/factorial(2*i+1) ;     // установка случайных коэффициентов
}
arg t;                          // объявление структуры параметров
for (x = 0; x < 10; x += 1)     // перебор значений х
{
t.l = 0; t.r = N;           // границы массива матриц
prod(&t);                   // вызов подпрограммы
cout << "\n Polynom value = " << t.q*x ;    // вывод
z = 0;
 
for (i = 0; i <= N; i++)                        // цикл для проверки
{
v = 1;
for (j = 0; j <i; j++) v = v*pow(x,2);
    // вычисление xi
z += x*v*a[i];                      // вычисление значения многочлена
}
cout << "== "<< z;                  // вывод значения для проверки
getch();                                // ожидание нажатия клавиши
}
}