#include <iostream>
#include <iomanip>
#include <cmath>

using namespace std;

void process(double Xn, double Xk, double dX, double eps) {
    while (Xn < Xk) {
        int n{};
        double res = 0, prev_res = 0, sum = 0;
        res = pow(Xn, 2 * n + 1) / (2 * n + 1);
        while (fabs(res - prev_res) > eps) {
            sum += res;
            prev_res = res;
            n++;
            res = pow(Xn, 2 * n + 1) / (2 * n + 1);
        }
        cout << Xn << " " << setprecision(4) << sum << " " << (n + 1) << " " << setprecision(4)
             << atanh(Xn) << endl;

        Xn += dX;
    }
}

int main() {
    double Xn, Xk, dX, eps;
    cin >> Xn >> Xk >> dX >> eps;

    process(Xn, Xk, dX, eps);

    return 0;
}