Численное интегрирование методом трапеций

#include "stdafx.h"
#include <iostream>
#include <iomanip>
#include <cmath>
#include <math.h>
#include <cstdlib>
using namespace std;
 
float f(float x)
{
    return 1/sqrt(pow(x,2)-3);
}
 
float fpr(float a, float b)
{
    return max( abs(3*pow(a,2)/pow(-3+pow(a,2), 5/2)-1/pow(-3+pow(a,2), 3/2)), abs(3*pow(b,2)/pow(-3+pow(b,2), 5/2)-1/pow(-3+pow(b,2), 3/2))  );
 
int _tmain(int argc, _TCHAR* argv[])
{
    setlocale(LC_ALL, "rus");
 
    float a,b, n, h, trap, t=0.000001;
    
    cout<<"Метод трапеций"<<endl<<endl;
 
    cout<<"Введите нижний предел  a="; cin>>a;
    cout<<"Введите верхний предел  b="; cin>>b;
    
 
    h=sqrt( (12*t)/((b-a)*fpr(a,b)) );
    n=(b-a)/h;
    
    float sum, sum2=0;
    for(int i=1; i<n; i++)
    {
        sum=f(a+i*h);
        sum2=sum+sum2;
    }   
 
    trap=0.5*(f(a)+f(b))*h+sum2*h;
    
    cout<<endl<<"Интеграл равен: "<<trap<<endl;
system("pause");
return 0;
}

#include "stdafx.h"
#include <iostream>
#include <iomanip>
#include <cmath>
//#include <math.h>
#include <cstdlib>
using namespace std;
 
float f (float x)
{
    return 1/sqrt(pow(x,2)-3);
}
 
float fpr(float a, float b)
{
    return max( abs(3*pow(a,2)/pow(-3+pow(a,2), 5/2)-1/pow(-3+pow(a,2), 3/2)), abs(3*pow(b,2)/pow(-3+pow(b,2), 5/2)-1/pow(-3+pow(b,2), 3/2))  );
}
 
//int _tmain(int argc, _TCHAR* argv[])
int main (void)
{
    setlocale(LC_ALL, "rus");
 
    float a,b, n, h, trap, t=0.000001;
 
    cout<<"Метод трапеций\n"<<endl;
 
    cout<<"Введите нижний предел  a="; cin>>a;
    cout<<"Введите верхний предел  b="; cin>>b;
 
 
    h=sqrt( (12*t)/((b-a)*fpr(a,b)) );
    n=(b-a)/h;
 
    float sum, sum2=0;
    for(int i=1; i<n; i++)
    {
        sum=f(a+i*h);
        sum2 += sum;
    }
 
    trap=0.5*(f(a)+f(b))*h+sum2*h;
 
    cout << "\nИнтеграл равен: " << trap << endl;
    system("pause");
    exit(0);
}