Вычисление интеграла

[DeathReject]

Всем привет. Нужна ваша помощь
Есть код программы для решения интеграла с помощью процедуры параллель 1/1+x*x
Нужно переделать код для решения интеграла f(x) = 7/(x3 + 1), количество процессов 8
с++ консольное приложение

#include "stdafx.h"
#include  <stdio.h>
#include  <omp.h>
int   main(int   argc ,   char   **argv)
 {
int   me,   size ;
               unsigned   intervals ;
               int   i,   j;
               double  h,   Ik,   x;
               printf ( " Enter   intervals : " ) ;
               scanf ( "%u" , &intervals ) ;
               h = 1.0/ intervals ;
               Ik = 0;
#pragma omp   parallel
{
#pragma omp for private (x)   reduction (+:Ik)
 for (j =0;   j  < intervals;   j++) 
               {
                  printf ( " thread N%d j=%d\n" ,   omp_get_thread_num () , j);
	x = (j + 0.5)*h;
                  Ik += 4/(1 + x  * x) ;
                  }
}
                  Ik  *= h;
printf ( " value = %.10 lf \n", Ik);
return 0;
}

[Ezhik Kihze]

Вы не видете где в коде написано 1/1+x*x? Или вставить туда свой ленитесь?