Описать метод Эйлера и обратный метод Эйлера

#include<iostream>
#include <cmath>
 
using namespace std;
 
double f(const double x, const double y){
    return y*exp(0.7*x) - exp(-1.5*x);
}
 
void P_Euler(double x, double b, double h, double y){
    double d, c, xd = 0, yd = 1, yd1 = 0, hd = h / 2;
    cout << "\nПрямий метод Ейлера \n";
    for (; x < b + h; x += h){
        y += h * f(x, y);
        for (int j = 0; j < 1; j++){
            yd1 = yd + hd * f(xd, yd);
            yd = yd1; xd += hd;
        }
        c = (yd1 - y) / (h*h - h*h / 4);
        d = c * h*h;
        cout << "x =  " << x << " |  y = " << y << " |  yd = " << yd1 << " |  c = " << d << "\n";
        cout << "--------------------------------------------------------------------------------" << "\n";
    }
}
 
 
 
void V_Euler(double x, double b, double h, double y){
    double d, c, xd = 0, yd = 1, yd1 = 0, hd = h / 2;
    cout << "\nОбратный метод Ейлера \n";
    for (; x < b + h; x += h){
        y += h * f(x, y);
        //y += h * f(x,y)/(1-h);
        for (int j = 0; j < 1; j++){
            yd1 = yd + hd * f(xd, yd);
            //yd1 = yd + hd * f( xd,yd)/(1-hd);
            yd = yd1; xd += hd;
        }
        c = (yd1 - y) / (h*h - h*h / 4);
        d = c * h*h;
        cout << "x =  " << x << " |  y = " << y << " |  yd = " << yd1 << " |  c = " << d << "\n";
        cout << "--------------------------------------------------------------------------------" << "\n";
    }
}
 
 
 
void Runge_Kutta(double x, double b, double h, double y){
    double y0 = 1, x0 = 0, k0, k1, k2, k3, c, d, hd = h / 2;
    double yd0 = 1, yd1 = 0, xd0 = 0, kd0, kd1, kd2, kd3;
    cout << "\nМетодРунге-Кутта \n";
    for (; x < b; x += h){
        k0 = h * f(x, y);
        k1 = h * f(x + h / 2, y + k0 / 2);
        k2 = h * f(x + h / 2, y + k1 / 2);
        k3 = h * f(x + h, y + k2);
 
        y += (k0 + 2 * k1 + 2 * k2 + k3) / 6;
 
        for (int j = 0; j < 2; j++){
            kd0 = hd * f(xd0, yd0);
            kd1 = hd * f(xd0 + hd / 2, yd0 + kd0 / 2);
            kd2 = hd * f(xd0 + hd / 2, yd0 + kd1 / 2);
            kd3 = hd * f(xd0 + hd, yd0 + kd2);
            yd1 = yd0 + (kd0 + 2 * kd1 + 2 * kd2 + kd3) / 6;
            xd0 += hd; yd0 = yd1;
        }
        c = abs(y - yd1) / (h*h*h*h - h*h*h*h / 16);
        d = c * h*h*h*h;
        cout << "x =  " << x << " |  y = " << y << " |  yd = " << yd1 << " |  c = " << d << "\n";
        cout << "--------------------------------------------------------------------------------" << "\n";
    }
}
 
int main(){
    setlocale(LC_ALL, "Russian");
    double a, b, h, y0;
 
    cout << " a = "; cin >> a;
    cout << " b = "; cin >> b;
    cout << " h = "; cin >> h;
    cout << " y(0)= "; cin >> y0;
 
    P_Euler(a, b, h, y0);
    V_Euler(a, b, h, y0);
    Runge_Kutta(a, b, h, y0);
 
    system("PAUSE>>void");
    return 0;
}