OpenMP, тройной интеграл

#include <iostream>
#include <Windows.h>
#include <omp.h>
 
#define PARALLEL
 
double inline Function(double value)
{
    return(value * value * value * value);
}
 
double FindInt(double low, double high, double step)
{
    double Value = 1;
 
    int Count = (high - low) / step + 1;
 
    for (int i = 0; i < 3; i++)
    {
        double Temp1 = 0;
        double Temp2 = 1;
        double temp = low;
#ifdef PARALLEL
#pragma omp parallel for reduction(+:Temp1, temp) reduction(*:Temp2)
#endif
        for (int i = 0; i < Count; i++)
        {
            Temp2 = 1;
 
            Temp2 *= temp;
            Temp2 *= temp;
            Temp2 *= temp;
            Temp2 *= temp;
            Temp2 *= step;
 
            Temp1 += Temp2;
 
            temp += step;
        }
        Value *= Temp1;
    }
 
    return(Value);
}
 
void main()
{
    int count = 0;
    double step;
    std::cout << "Input Step: ";
    std::cin >> step;
 
    DWORD Start = GetTickCount();
    std::cout << FindInt(0,1,step);
    std::cout << std::endl << "Time: " << (GetTickCount() - Start);
 
    std::cin >> step;
}