помогите с функциями в С++ СodeBlocks

[Valeriyaa]

нужно вычислить интеграл методом Симпсона и трапеции
подскажите что делать дальше и в чём ошибки?
файл main.cpp:

#include <iostream>
#include "integrall.h"
inline double ff(double x){return ((x+1)*(x+1)*(x+1);}
 
using namespace std;
 
int main()
{cout<<integrall(exp,1,0)<<endl;
cout<<integrall(ff,-3,5)<<endl;
 
    return 0;
}

файл integrall.cpp:

#include <iostream>
#include "integrall.h"
 
 
using namespace std;
 
double integrall (p_Func,double a,double b,doubl eps=(E-7));
//int main()
{double xi,dx=b-a;
double R,I=(pf(a)+pf(b))*dx*0.5;}
 
{
    R=0;
    for (xi=a+dx*0.5;xi<b;xi+=dx)
        R+=dx*pf(xi);
    I=(I+R)*0.5;
    dx/=2;
}
while (fabs(I-R)<2*eps)
 
    return I;
}

файл integrall.h:

#ifndef INTEGRALL_H_INCLUDED
#define INTEGRALL_H_INCLUDED
#define double(*pf)(double)
double integrall (double(*pf)(double),double a,double b,double eps=(E-7));
 
 
 
#endif // INTEGRALL_H_INCLUDED