#include <stdio.h>
#include <math.h>

double f(double x){
 //   return cos(x);
    return x*x*x*x;
}
double riemann(int n,double a,double b){
    double h=(b-a)/n;
    double r=0;
    for(int k=0;k<n;k++){
        r+=f(a+h*k);
    }
    return r*h;
}
double richard(int n,double a,double b){
    return 2*riemann(2*n,a,b)-riemann(n,a,b);
}
double richard2(int n,double a,double b){
    return (4*richard(2*n,a,b)-richard(n,a,b))/3;
}
int main(){
    int n=1;
    for(int i=0;i<20;i++){
        printf("%d %g %g %g\n",n,
            fabs(riemann(n,0,1)-1.0/5),
            fabs(richard(n,0,1)-1.0/5),
            fabs(richard2(n,0,1)-1.0/5));
        n*=2;
    }
    return 0;
}