Параллельное программирование: нахождение интеграла

// intagral_omp.cpp : Defines the entry point for the console application.
//
 
#include "stdafx.h"
#include <omp.h>
#include <iostream>
#include <time.h>
#include <math.h>
#include <conio.h>
using namespace std;
#define n 2000
 
 
double f(double x)
{
    return (x*x)*4*pow(sin(x),5);
}
    
int _tmain(int argc, _TCHAR* argv[])
{
    omp_set_dynamic(0);
    omp_set_num_threads(4);
    double a=-5.0;
    double b=10.0;
    double h;
    double* x=new double[n];
    double result=0;
    int thr=0;
 
 
    x[0]=a;
    x[n]=b;
    h=fabs(a-b)/n;
    
    cout<<"Integral calculator"<<endl;
    cout<<"\n\na = "<<a<<endl;
    cout<<"\nb = "<<b<<endl;
    cout<<"\nh = "<<h<<endl;
    cout<<"\n\n Number of divisions = "<<n<<endl;
    for(int i=0;i<=n;i++)
    {
        x[i]=a+i*h;
        //cout<<x[i]<<endl;
    }
    
#pragma omp parallel for reduction(+: result) 
    for(int i=1;i<=n;i++)
    {
        result+=h*f((x[i-1]+x[i])/2);
        if(omp_get_num_threads()>=thr)
        {
        thr=omp_get_num_threads();
        }
 
    }
    cout<<"\n\nMax number of used threads: "<<thr<<endl;
    cout<<result;
    //system("pause");
    _getch();
 
 
    return 0;
}