[C++] Sprawdzenie poprawności użycia OpenMP przy całkowaniu

#include <iomanip>
#include <iostream>
#include <cstdlib>
#include <cmath>
#include <omp.h>

using namespace std;

//*******************************
//** Tutaj definiujemy funkcję **
//*******************************

double f( double x, int a, int b, int c )
{
    return( a * pow( x, 2 ) + b * x + c );
}
//********************
//** Program główny **
//********************

int main()
{
    const int N = 1000; //liczba punktów/trapezów podziałowych
    double xp, xk, s, dx;
    int a, b, c, i;
   
   
    cout << setprecision( 3 ) // 3 cyfry po przecinku
    << fixed; // format stałoprzecinkowy
   
    cout << "Podaj wartosc a \n";
    cin >> a;
   
    cout << "Podaj wartosc b \n";
    cin >> b;
   
    cout << "Podaj wartosc c \n";
    cin >> c;
   
    cout << "f(x) = a * x * x + b * x + c x" << endl
    << "Podaj poczatek przedzialu calkowania" << endl
    << "xp = ";
    cin >> xp;
    cout << "Podaj koniec przedzialu calkowania" << endl
    << "xk = ";
    cin >> xk;
   
    #pragma omp parallel firstprivate(xk, xp) private (dx, i, a, b, c) shared (s)
    {
       
        s = 0;
        dx =( xk - xp ) / N;
       
        #pragma omp for reduction (+:s)
        for( i = 1; i < N; i++ ) {
           
            s += f( xp + i * dx, a, b, c );
            s *= dx;
            cout << "Wartosc calki wynosi : " << setw( 8 ) << s
            << endl << endl;
           
        }
    }
    system( "pause" );
    return 0;
}