#include<cmath>#include<cstdio>#include<iostream>usingnamespace std;#define N 20int n;longdouble a, b, e, T[N][N]={{0.0}};longdoublef(longdouble x){return1.0/(1.0+ x);}longdoublex(int i,longdouble h){return a + h * i;}intmain(){scanf("%llf%llf%llf%d",&a,&b,&e,&n);int k =0;for(; k < n; k++){longdouble h =(b - a)/pow(2, k), sum =0.0;for(int i =1; i <=pow(2, k)-1; i++) sum +=f(x(i, h));
T[k][0]=0.5* h *(f(a)+2* sum +f(b));for(int m =1; m <= k; m++)
T[k][m]=(pow(4, m)* T[k][m -1]- T[k -1][m -1])/(pow(4, m)-1);for(int m =0; m <= k; m++)printf("%.7Lf\t", T[k][m]);
cout << endl;if(k >0&&fabs(T[k][0]- T[k][k])<= e)break;}return0;}