Помогите Исправить Ошибку

#include <iostream>
#include <cmath>
using namespace std;
double f1(double);
double f2(double);
double f3(double);
double B(double,double, double, double (*f)(double));
double E(double,double, double, double (*f)(double));
 
    int main(){
    double a0, a1, eps;
    double (*f_array[])(double ,double, double, double (*)(double))={B,E}; 
    double (* integ [])(double)={f1, f2, f3}; 
    int a, u;
        cout<<"Input eps "<<endl;
        cin>>eps;
        cout<<"Input a0, a1 for function "<<endl;
        cin>>a0>>a1;
        cout<<"Input the number of methods from 2, 5 "<<endl; 
        cin>>a;
        cout<<"Input the number of integral from 1 to 3 "<<endl;
        cin>>u;
        printf("%.10f",(f_array[a-1])(eps, a0, a1, integ[u-1]));
        }
 
 
double B(double eps,double a, double b, double (*f)(double))
{
    int n=4; 
    double h=(b-a)/4,k1, k0, l;
        for(int i=1; i<=n; i++)
            {
                l=a+i*h; k1+=h*f(l);
            }
                do{
                    k0=k1;
                    k1=0;
                    n*=2;
                    h=(b-a)/n;
                    for(int i=1; i<=n; i++)
                    {
                        l=a+i*h;
                        k1+=h*f(l);
                    }
            } 
                while (abs(k1-k0)>eps);
                return k1;
}
 
double E(double eps,double a, double b, double (*f)(double))
{
    int n=4;
    double h=(b-a)/4, k1, k0, l, m;
        for(int i=1; i<=n; i++)
        {
            l=a+(i-1)*h; k1+=h*f(l);
        }
            do
            {
            k0=k1;
            k1=0;
            n*=2;
            h=(b-a)/n;
                for(int i=1; i<=n; i++)
                {
                    l=a+(i-1)*h;
                    m=a+i*h; 
                    k1+=(m-l)/6*(f(l)+4*f((l+m)/2)+f(m));}
                } 
            while (abs(k1-k0)>eps);
            return k1;
}
 
double f1(double x)
{
    return (x*sqrt(1+x));
}
 
double f2(double x)
{
    return (asin(sqrt(x))/(sqrt(x*(1-x))));
}
 
double f3(double x)
{
    return (pow(2,x)/(1+pow(4,x)));
}