По заданным действительным числам вычислить интегралы

#include <iostream>
#include <iomanip>
#include <fstream>
# include <cmath>
# include <time.h>
# define M_PI 3.14159265358979323846264338327950
# define N 2000 // максимально возможное число точек в которых проводим вычисления 
# define r 7 // количество узлов в формуле гаусса
using namespace std;
 
ofstream ff("res.txt");
 
int n=20; // число точек в которых проводим вычисления 
double usl_x[N];
double X[N][r];// хранит значения узлов разбиения по х
double x[N];
double koef_gauss[r]={0.129484966, 0.279705391, 0.381830051, 0.417959184, 0.381830051, 0.279705391, 0.129484966};
double korni_mhogoch_legandra[r]={-0.949107912, -0.741531185, -0.405845151, 0, 0.405845151, 0.741531185, 0.949107912};
 
 
double f(double x0, int v)// v- номер интегралла
{
if(v==1) return pow(atan(x0),2.);
    else return sin(exp(10*x0));
}
 
double integr(int v, double a ,double b)
{
    double h;
    int n1=300;
    h=abs((b-a))/n1;
    x[0]=a;
    for(int i=1;i<=n1;i++) x[i]=x[i-1]+h;
    
    for(int i=0;i<n1;i++)
        for(int j=0;j<r;j++) X[i][j]=(x[i+1]+x[i])/2.+((x[i+1]-x[i])/2.)*korni_mhogoch_legandra[j];
 
    double sum=0,sum1=0;
    for(int i=0;i<n1;i++)
    {
        sum1=0;
        for(int j=0;j<r;j++)
        {
            sum1+=koef_gauss[j]*f(X[i][j],v);
        }
        sum+=(h/2.)*sum1;
    }
    
    return sum;
}
 
void setka(double a, double b)
{
    double h;
h=(double)(b-a)/n;
for (int i=0;i<=n;i++)    usl_x[i]=a+i*h;
}
 
void main()
{
    
    setlocale(LC_ALL, "Russian");// подключение русской раскладки для вывода всех сообщений на русском языке
    cout<<"Введите а, b,n\n";
    double a,b,n;
    double res;
    cin>>a>>b>>n;
    setka(a,b);
    res=integr(1,a,b);
    setka(0,n);
    res+=integr(2,0,n);
    cout<<"Результат\n"<<res<<endl;
    system("pause");
}