Вычислить значение 3х интегралов с заданной точностью С++

[ksly]

Даны 3 интеграла.Необходимо написать программу,вычисляющую их значение с заданной точностью.Головная функция, обеспечивающая ввод с клавиатуры точности вычислений (в пределах 10^-1... 10^-7). Вывести количество разбиений, на котором достигается заданная точность. Проблемы с подчеркнутым

#include "stdafx.h"
#include <iostream>
#include <cmath>

using namespace std;
double F1(double x);
double F2(double x);
typedef double(*pointFunc)(double);


double M_Simpson(pointFunc F, double a, double b, unsigned N)
{
	double Res = 0.0;
	double h = (b - a) / N;
	for (double x = a + h; x<b; x += 2 * h)
		Res += 2 * (*F)(x);
	for (double x = a + 2 * h; x<b; x += 2 * h)
		Res += 4 * (*F)(x);
	Res = h*((*F)(a)+(*F)(b)+Res) / 3;
	return Res;
}

double F1(double x)
{
	return 1 / (x + 1)*(sqrt(pow(x, 2) + 1));
	
}

double F2(double x)
{
	return sqrt(pow(2, x) - 1);
}

double F3(double x)
{
	double e = 2.71;
	return sqrt(pow(e, x) - 1);
}

int main()
{
	cout <<"I1="<<M_Simpson(&F1)<< endl;
	cout <<"I2="<<M_Simpson(&F2)<< endl;
	cout <<"I3="<<M_Simpson(&F3)<< endl;
	system("pause");
	return 0;
}