CUDA_full / P04_Integral / CUDA_Integral / integral

#include "integral.h"


// Подинтегральная функция: 4 / (1 + x^2).
float func(float x)
{
    return 4.0f / (1.0f + x * x);
}


/* Вычисление интеграла методом прямоугольников,
a - начало расчетной области,
n - размерность сетки (количество интервалов),
h - шаг сетки.
*/
float integral(float a, int n, float h)
{
    float result = 0.0f;

    for (int i = 0; i < n; ++i)
        result += func(a + ((float)i + 0.5f) * h);
    result *= h;

    return result;
}