Найти ошибку в программе. Интеграл по формуле Симпсона,

// Lab14Y.cpp: определяет точку входа для консольного приложения.
//
 
#include "stdafx.h"
#include <iostream>
#include "math.h"
#include "conio.h"
#pragma warning(disable : 4996)
#define M_PI 3.1415926535897932384626433832795
 
using namespace std;
 
 
 
double f(double x)
{
    return sqrt(pow(x, 2) + 0.16) / x;
}
 
 
 
 
double c(double(*func)(double), double a = 1, double b = 2)
{
    double s = 0;
    for (double i = a; i<b; i += 0.01)
    {
 
        s += ((*func)(i)+(*func)(i + 0.01))*0.5*0.01;
    }
 
    return s;
}
 
 
 
double Simps(int n) // formula Simpsona
{
    double m;
    double x, x1, x2, h, sum = 0;
    m = n >> 1;
    h = (double)1 / n;
    for (int i = 1; i <= m; i++)
    {
        x = 2 * i*h;
        x1 = (2 * i + 1)*h;
        x2 = (2 * i + 2)*h;
        sum += (f(x) + 4 * f(x1) + f(x2));
    }
    return h*sum / 3;
}
 
 
 
void main()
{
    
    double a = 0;
    double b = M_PI/3;
    int g = 2.2;
    cout << c(f, a, b) << endl;
    cout << c(f, a) << endl;
    cout << c(f) << endl;
    cout << Simps(g) << endl;
 
    
}

// Lab14Y.cpp: определяет точку входа для консольного приложения.
//
 
#include "stdafx.h"
#include <cmath>
#include <iostream>
#include <conio.h>
#pragma hdrstop
#define pi 3.1415926535897932384626433832795
 
using namespace std;
 
//Підінтегральна функція
double func(double x)
{
    double r = sin(x);
    
    return r;
}
 
 
double simpson(double(*Fx)(double), double a, double b, int n);
double newton(double(*Fx)(double), double a, double b);
 
int main(int argc, char* argv[])
{
    double I;
    double S;
    double a, b;
    int n;
 
    cout << "\n Input integral atributes:\n" << endl;
    a = 0;
    b = pi / 3;
 
    cout << "\t-> Enter n for Method Simpson = ";
    cin >> n;
 
    // вичислення інтеграла
 
    I = simpson(&func, a, b, n);
    S = newton(&func, a, b);
    // Вивід результату
    cout << endl << "Method Newton \n\n I= " << S;
    cout << endl << "\nMethod Simpson \n\n I= " << I;
    cin.get();
}
 
 
// Метод Сімпсона
double simpson(double(*Fx)(double), double a, double b, int n)
{
 
    double h;
    h = (b - a) / n;
 
    double I, I2 = 0, I4 = 0;
    I4 = Fx(a + h);
    for (int k = 2; k < n; k += 2)
    {
        I4 += Fx(a + (k + 1)*h);
        I2 += Fx(a + k*h);
    }
    I = Fx(a) + Fx(b) + 4 * I4 + 2 * I2;
    I *= h / 3;
 
    return I;
}
 
// Метод Ньютона Лейбніца
double newton(double(*Fx)(double), double a, double b)
{
    double s = 0;
    for (double i = a; i<b; i += 0.000001)
    {
 
        s += ((*Fx)(i)+(*Fx)(i + 0.000001))*0.5*0.000001;
    }
 
    return s;
}