#include <cstdio>
#include <iostream>
#include <cmath>
#include <algorithm>
using namespace std;
double F(double x){
return x * x + x;
}
double simpson(double a,double b){
double c = a + (b - a) / 2.0;
return (F(a) + 4 * F(c) + F(b)) * (b - a) / 6.0;
}
double asr(double a,double b,double eps,double A){
double c = a + (b - a) / 2.0;
double L = simpson(a,c), R = simpson(c,b);
if(fabs(L + R - A) <= 15 * eps) return L + R + (L + R - A) / 15.0;
return asr(a,c,eps/2.0,L) + asr(c,b,eps/2.0,R);
}
double asr(double a,double b,double eps){
return asr(a,b,eps,simpson(a,b));
}
int main(){
cout<<asr(0,1,0.000001)<<endl;
cout<<5.0/6<<endl;
return 0;
}