#include <iostream>
#define _USE_MATH_DEFINES
#define _USE_GENERIC_MATH1
#include <math.h>
using namespace std;

int N = 100;
double I_anslit = (sqrt(2)/4) * (2*atan(sqrt(2)/2) + 2*atan(sqrt(2)));
double a = 0;
double b = 1.0;

double f(double x) {
return (1 + x*x) / (1 + x*x*x*x);
}

double F() {
double I = 0, h;
double x, alpha;
h = (b - a) / N;
alpha = a + (h / 2);
for (int i = 1; i <= N; i++) {
x = a + i * h;
I += f(alpha + (i - 1) * h);
}
I *= h;
return I;
}

int main() {
cout << "I_anslit= " << I_anslit << endl;
cout << "q = " << F() << endl;
cin.get();
return 0;
}

double F() {
 double I = 0, h;
 double x, alpha;
 h = (b - a) / N;
 alpha = a + (h / 2);
 for (int i = 1; i <= N; i++) {
  x = a + i * h;
  I += f(alpha + (i - 1) * h);
 }
 I *= h;
 return I;
}