Проверить значение первообразной интеграла Римана

# include <iostream>
# include <conio.h>
# include <math.h>
using namespace std;
double m, c, d, d1;
double f1(double);
double f2(double);
double F1(double);
double F2(double);
double integral (double a, double b, int n, double(*f)(double))
{ int k;
double s, h, xk; s = 0; h = (b - a) / n; xk = a + h;
for (k = 1; k <= n - 1; k++)
{
    s = s + f(xk); xk += h;
}
s = h*(f(a) + f(b)) / 2 + s;
return s;
}
double f1(double x)
{
    return ((exp (3 * x) + c) / (exp( x) + d)) + (c - d)*x;
}
double f2(double x)
{ 
    return sqrt(c - (0.16 / x*x)) + (c*c - d1*d1 + m)*x;
}
double F1(double x)
{
    return (exp(2*x)/2)-(exp(x)+x);
}
double F2(double x)
{ 
    return sqrt(c - 0.16)- 0.4*(acos( 0.4/x ));
}
void main()
{
    cin >> c>> d1 >> d >> m;
    double i1, i2,I1,I2;
    i1 = integral(0, 2, 100, f1);
    i2 = integral(1, 2, 100, f2);
******* // здесь нужно проверить правильность ответа по формуле интегрaла F(b)- F(a) для 2 функций F1 и F2 НЕ знаю как сделать( 
    cout << "i1=" << i1 << " " << "\n" << "i2=" << i2 << "\n"<< "I1=" << I1 << " " << "\n" << "I2=" << I2 << "\n";
 
    _getch();
}

cout<<i1<<" "<<F1(2)-F1(0)<<endl;
cout<<i2<<" "<<F2(2)-F2(1)<<endl;