Расчет числа pi - параллельное программирование

[VolodyaBuzin]

Добрый день! Проверьте, правильно ли я написал программу, если нет, в комментариях к коду покажите, что нужно сделать на самом деле. Как правильно распараллелить цикл прописать количество ядер? программа иногда выдаёт результат, отличный от 3.14, например, 1,4752

// Integral.cpp : Defines the entry point for the console application.
//
 
#include "stdafx.h"
#include <omp.h>
#include <time.h>
#include <iostream>
#include <cstdio>
#include <cmath>
#include <conio.h>
using namespace std;
 
double InFunction(double x) //Подынтегральная функция
{
  return (4.0/(1 + pow(x, 2)));
}
 
double CalcIntegral(double a, double b, int n, int k)
{
omp_set_num_threads(k);
  int i;
  double result, h;
 
  result = 0;
 
  h = (b - a) / n; //Шаг сетки
#pragma omp parallel for reduction(+: result)
  for(i = 0; i < n; i++)
  {
    result += InFunction(a + h * (i + 0.5)); //Вычисляем в средней точке и добавляем в сумму
  }
 
  result *= h;
 
  return result;
}
 
 
int _tmain(void)
{
    clock_t st, et;
      double integral, iv;
      int n, k;
  cout << "Enter N:\n";
  cin >> n;
  cout << "Enter number of threads:\n";
  cin >> k;
  st = clock();
  integral = CalcIntegral(0, 1, n, k);
  et = clock();
  iv = (double)(et-st)/CLOCKS_PER_SEC;
  printf("time is %f\n", iv);
  printf("The value of the integral is: %lf \n", integral);
    _getch();
    return 0;
}