#include <iostream>

using namespace std;

double f(double x, double y) {
    return x+x*y+y+1;
}

void prosty(double x, double y, double h, int k) {
    cout << "Prosty Euler" << endl << endl;
    for (int i = 1; i <= k; i++) {
        y = y + h*f(x, y);
        x = x + h;
        cout << "x" << i << ": " << x << endl << "y" << i << ": " << y << endl;
    }
}

void ulepszony(double x, double y, double h, int k) {
    cout << "Ulepszony Euler" << endl << endl;
    for (int i = 1; i <= k; i++) {
        y = y + h/2*(f(x, y) + f(x+h, y+h*f(x,y)));
        x = x + h;
        cout << "x" << i << ": " << x << endl << "y" << i << ": " << y << endl;
    }
}



int main() {
    double h = 0.25;
    double x = -1;
    double y = 1;
    int k = (0 - x)/h;

    cout << "Liczba krokow: " << k << endl << endl;

    prosty(x, y, h, k);
    cout << endl << "------------------" << endl << endl;
    ulepszony(x, y, h, k);

    return 0;
}