Вычисление интеграла по формуле трапеций

/*Параметры функций : 
 @x - значение переменной x 
 @a - нижний придел интегрирования 
 @b - верхний придел интегрирования
 @n - количество кусочков, на которое разбивается интеграл 
*/
double getSimpson(long double x,double a,double b , unsigned int n ) ;
double getTrap(long double x,double a,double b , unsigned int n );

double getSimpson(long double x,double a,double b , unsigned int n )
{
    double dStep,      //Шаг
           h ;             //Величина шага
 
          int i  ;          //Переменная цикла
 
    double F[n],       //Массив, содержащий значение площади , каждого из n-отрезков
           dResult,     //Конечный результат
           dTmp;       //Промежуточный результат
 
 
    h = (b-a) / n;     //Расчитываем величену  шаг
 
    for(i=0, dStep = a; i<n; i++, dStep+=h )    //В кач. начала отр. берем нижний придел интеграла и все время перемещаемся на Шаг
        F[i] = 1 / (x*log(x) );                //Сохраняем знач. площади , каждого из кусочков
 
    dResult = 0;
 
    for(i=1; i<n; i+=2)                         //Цикл 1,3,5,7 ... n
     dResult += F[i];                           //Записываем сумму нечетных отрезков
 
    dResult *= 4;                               //Умножаем сумму нечетных отрезков на 4
    dTmp = 0;
    for(i=2; i<n-1; i+=2)                       //Цикл 2,4,6,8,10 ... n-1
     dTmp += F[i];                              //Записываем сумму четных отрезков
     dTmp*= 2 ;                                 //Умножаем сумму четных отрезков на 2
 
    dResult += dTmp + F[0] + F[n-1];            //К сумме нечетных прибовляем сумму четных и учитываем нулевой и предпосл. элементы
 
    dResult *= h/3;                             //Весь результат умножаем на Шаг/3
 
    return dResult ;                            //Возврашаем площать интеграла
}

double getTrap(long double x,double a,double b , unsigned int n )
{
    double dStep,                   //Шаг
           h ;                      //Величина шага
    int i  ;                        //Переменная цикла
    double F[n],                    //Массив, содержащий значение площади , каждого из n-отрезков
           dResult ;                //Конечный результат
 
    h= (b-a) / n;                   //Расчитываем величену  шаг
 
    for(i=0, dStep = a; i<n; i++, dStep+=h )     //В кач. начала отр. берем нижний придел интеграла и всевремя перемещаемся на Шаг
        F[i] = 1 / (x*log(x) );                 //Сохраняем знач. площади , каждого из кусочков
 
    dResult = 0;
 
    for(i=1;i<n-1;i++)
     dResult += F[i];                           // Суммируем значения отрезков, не включая 1-ый и предпоследний
    dResult += F[0]/2 + F[n-1]/2;               // Пребовляем к результату первый и последний значения отрезков , деленные на 2
    dResult *= h;                               // Умножаем рузультат на значение шага
 
    return dResult ;                            //Возврашаем площать интеграла
}