#include <iostream>
#include <cmath>
using namespace std;

double F(double x, double y) {
    return (cos(y))/(x+2);
}

int main() {
    double a = 0; double b = 1; double h = 0.1;
    double n = (b - a) / h;
    double X[(int)n];
    double Y1[(int)n];
    double Y2[(int)n];
    double Y3[(int)n];
    double Y4[(int)n];
    double Y[(int)n];
    //calculate
    X[0] = 0; Y[0] = 0;
    for (int i = 1; i <= n; i++) {
        X[i] = a + i * h;
        Y1[i] = h * F(X[i - 1], Y[i - 1]);
        Y2[i] = h * F(X[i - 1] + h / 2.0, Y[i - 1] + Y1[i] / 2.0);
        Y3[i] = h * F(X[i - 1] + h / 2, Y[i - 1] + Y2[i] / 2);
        Y4[i] = h * F(X[i - 1] + h, Y[i - 1] + Y3[i]);
        Y[i] = Y[i - 1] + (Y1[i] + 2 * Y2[i] + 2 * Y3[i] + Y4[i]) / 6;
    }
    //print results
    for (int i = 0; i <= n; i++) {
        cout << "X[" << i << "]=";
        cout << fixed << setprecision(5) << X[i] << endl;
    }
    cout << endl;
    for (int i = 0; i <= n; i++) {
        cout << "Y[" << i << "]=";
        cout << fixed << setprecision(5) << Y[i] << endl;
    }
    return 0;
}