Определённый интеграл методом Гаусса

#include<iostream>
#include<math.h>
#include <conio.h>
using namespace std;
int kol_obr;
double y(double x);//подынтыгральная функция
double perv(double x);//первообразная
double gauss (double a,double b,double eps);//метод Гаусса
void main()
{
setlocale(LC_ALL,"rus_rus.1251");
double a=0.2,b=3.3,eps,prib_gauss,prv,perv;
cout<<"Введите точность"<<endl;
cin>>eps;
prib_gauss=gauss(a,b,eps);//обращение к функции Гаусса
cout<<"Метод Гаусса "<<prib_gauss<<endl;
prv=perv(double x);//обращение к первообразной
cout<<"Первообразная "<<prv<<endl;
_getch();
}
 
double perv(double x)
{
return (1./4)*((x*x*x*x-1)*log(x*x+1))-(1./4)*(((x*x*x*x)/2)+x*x);
}
 
double gauss (double a,double b,double eps)
{
int n=1, k=1;;
double h,x,s=0,INT1,INT2=1,INT;
do
{
INT1=INT2;
k+=2;
h=(b-a)/k;//шаг интегрирования
x=a+0.5*h;
for(int j=2;j<k;j++)
s+=(x*x*x)*log(x*x+1);
INT=h*s;
}
while(abs(INT2-INT1)>eps);
return INT;
}
double y(double x)
{
kol_obr++;
return (x*x*x)*log(x*x+1);
}

INT2=h*s;

return INT2;