Вычисление интеграла

#include <iostream>
 
 
 
double Stepen(double a, int b) 
{
    int tmp = a;
    for (int i = 0; i < b; i++) 
    {
        a *= tmp;
    }
    return a;
}
 
int main()
{
    setlocale(LC_ALL, "Russian");
    
    double res;
    double xUp, xDown;
    std::cout << "Введите нижний аргумент интеграла: "; std::cin >> xDown; std::cout << "Ввеидте верхний аргумент интеграла: "; std::cin >> xUp;
 
    res = (1 / Stepen(xUp, 4) + 10 + Stepen(xUp, 2) + 9) - (1 / Stepen(xDown, 4) + 10 + Stepen(xDown, 2) + 9);
 
    std::cout << res;
    
    system("pause");
    return 0;
}

#include <iostream>
#include <cmath>
 
float func(float x)
{
    return 1. / (std::pow(x, 4) + 10 * (x * x) + 9);
}
 
int main()
{
    float a = 0; // lower bound
    float b = 1; // upper bound
    float d = 0.01; // step
    int n = static_cast<int>((b - a) / d);
    double result = 0;
    for (int i = 1; i <= n; ++i)
        result += d * func(a + d * (i - 0.5));
 
    std::cout << "result: "<< result << std::endl;
}

res = (1/( Stepen(xUp, 4) + 10 + Stepen(xUp, 2) + 9)) - (1/(Stepen(xDown, 4) + 10 + Stepen(xDown, 2) + 9));