Автоматический выбор шага и метод Рунге-Кутта - C++ - Ответ 6356982

#include <iostream>
using namespace std;
float fun(float x, float y)
{
    return x*x + y*y;
}
 
void run_kut(float a, float b, float y0, int n)
{
    float x, y, x1, y1, h, k1, k2, k3, k4;
    x = a;
    y = y0;
    h = (b - a) / n;
    for (int i = 0; i<n; i++)
    {
        k1 = h*fun(x, y);
        x1 = x + h / 2.;
        y1 = y + k1 / 2.;
        k2 = h*fun(x1, y1);
        y1 = y + k2 / 2.;
        k3 = h*fun(x1, y1);
        x += h;
        y1 = y + k3;
        k4 = h*fun(x, y1);
        y += (k1 + 2. * k2 + 2. * k3 + k4) / 6.;
        cout << "x[" << i + 1 << "]=" << x << "\t" << "y[" << i + 1 << "]=" << y << "\n";
    }
}
void main()
{
    setlocale(LC_ALL, "Russian");
    int n;
    float a, b, y0,e;
    system("Cls");
    cout << "Ведите начало a="; cin >> a;
    cout << "Ведите конец b="; cin >> b;
    cout << "Ведите начальное положение y0="; cin >> y0;
    cout << "Ведите точность e="; cin >> e;
    n = (b-a)/pow(e,1./4);//тут конечно лучше приводить к наибольшему целому(в такойм варианте получается приведение к наименьшему целому)
    cout<<"Для заданной точность, количесво разбиений = " << n<<endl;
    cout << "\nМетод Рунге-Кутты:\n";
    run_kut(a, b, y0, n);
    cout << "\n";
    system("pause");
}