Программа для нахождения определенного интеграла

#include <iostream>
#include<math.h>
#include<iomanip>
 
void metod1(double,double,double);
void metod2(double,double,double);
void metod3(double,double,double);
// x= x*atan(x), a=0, b=sqrt(3.0);
using namespace std;
int main()
{
    setlocale (LC_ALL, "ukr");
 
void (*f)(double,double,double);
double a,b,n,s;
int p;
b=sqrt(3.0);
a=0;
n=1000;
cout<<"Промiжок ("<<a<<" ; "<<b<<")\n"<<endl;
cout<<"Виберiть метод"<<endl<<" \n 1. Метод трикутникiв \n 2. Метод трапецiй \n 3. Метод Сiмпсона \n 0. Вихiд"<<endl<<" \n Зробiть свiй вибiр.\n"<<endl;
while(1){
cin>>p;
if (p==1){f=metod1; f(a,b,n);}else if(p==2){f=metod2; f(a,b,n);}else if(p==3){f=metod3; f(a,b,n);}else if(p==0) break;else cout<<" Неправильне значення"<<endl;
    }
system("pause");
return (0);
}
 
void metod1(double q,double w, double e)
{
double h,x,r;
h=(w-q)/e;
r=0;
x=q;
for (int i=0;i<e;i++)
{
r+=x*atan(x);
x+=h;
}
r*=h;
cout<<setiosflags(ios::fixed)<<setprecision(3)<<setw(6)<<r<<"\n";
}
 
void metod2(double q,double w, double e)
{
double h,x,r;
h=(w-q)/e;
r=(atan(q)+atan(w))/2;
x=q+h;
for (int i=0;i<e;i++)
{
r+=x*atan(x);
x+=h;
}
r*=h;
cout<<setiosflags(ios::fixed)<<setprecision(3)<<setw(6)<<r<<"\n";
}
 
void metod3(double q,double w, double e)
{
double h,x,r;
h=(w-q)/(2*e);
r=atan(q)-atan(w);
x=1;
for (int i=1;i<=(2*e-1);i++)
{
r+=(3+x)*(q+i*h)*atan(q+i*h);
x=(-1)*x;
}
r*=(h/3);
cout<<setiosflags(ios::fixed)<<setprecision(3)<<setw(6)<<r<<"\n";
}