Интеграл функции методом трапеций, прямоугольника и Симпсона

#include <iostream>
#include <conio.h>
#include <iomanip>
 
 
 
 
using namespace std;
float Pr(float f);
float Tr(double f);
float Si(double F);
 
 
int b,e=1000,n,k;
    double a,f=1.0,tmp1,tmp2,h,x=1.0,P=0,s=0,F;
 
void main()
{
    tmp1=1+pow(3*x,8.0);
    tmp2=pow(tmp1,0.5);
    f=(pow(x,15.0)*tmp2);
        cout<<"vvedit grani [a,b]"<<endl;
    cin>>a; cout<<endl;
    cin>>b; cout<<endl;
    cout<<"Vvedit kilkist krokiv";
    cin>>k;
        
        
    a=Tr(a);
    cout<<a;
    
 
 
            _getch();
 
}
 
    
float Pr(int a,int b, unsigned k,double h, double(f)){
    
     h=(b-a)/k;
    f=(pow(x,15.0)*a+tmp2*b);
        for(x=a;x<=b;x+=h)
            f+=pow(x,15.0)*tmp2;
        f*=h/e;
        return float(f);
}
 
 
 
 
 
float Tr(int a,int b, unsigned k,double h, double(f)){
 
    
    
 
         h=(b-a)/k;
        
        f=(pow(x,15.0)*a+tmp2*b)/2;
        for(x=a+h;x<=b;x+=h)
            f+=pow(x,15.0)*tmp2;
        f*=h/e;
        return float(f);
}
 
 
 
 
 
 
float Si(double a, double b, unsigned k,double (*f)(double x), double (F))
    {
        
        
        double h=(b-a)/k;
        double x0=a;
        double x1=a+h;
        
        
        for(unsigned i=0;i<=k-1;i++){
            s=s+f(x0)+4*f(x0+h/2)+f(x1);
            x0+=h;
            x1+=h;
        }
        F=(h/6)*s;
        return float(F);
    }

float Tr( double (*f),int a=0,int b=1, unsigned k=1,double h=0.1);

float Tr( double (*f),int a,int b, unsigned k,double h){

double func(double x)
{
    double r=pow(x,14)*(57*pow(x,8)+15)*sqrt(3*pow(x,8)+1);
    return r;
}
double SimpsonMethod(double W, double Z, double (*r)(double x)) {
    int k=1000;
    double h =(W-Z)/k;
    double sum = 0;
 
    double x0 = W;
    double x1 = W + h;
 
    for (int i=0; i<=k-1; i++) {
        sum += r(x0) + 4*r(x0 + h/2) + r(x1);
 
        x0 += h;
        x1 += h;
    }
 
    return (h/6)*sum;
}