Как написать программу для метода Рунге-Кутты 4 порядка с автоматическим выбором шага ?

#include <iostream>
#include <math.h>
#include <conio.h>
 
using namespace std;
 
float dydx(float x, float y){
    return y*y - 2 / (x*x);
}
 
// Находит значение y для заданного x, 
//используя размер шага h и 
//начальное значение y0 при x0.
float rungeKutta(float x0, float y0, float x, float h){
    // Подсчитайте количество итераций с 
    //использованием размера шага h
    int n = (int)((x - x0) / h);
 
    float k1, k2, k3, k4, k5;
    float y = y0;
    for (int i = 1; i <= n; i++)
    {
        // формулы Рунге Кутты, для нахождения 
        //следующего значения y
        k1 = h*dydx(x0, y);
        k2 = h*dydx(x0 + 0.5*h, y + 0.5*k1);
        k3 = h*dydx(x0 + 0.5*h, y + 0.5*k2);
        k4 = h*dydx(x0 + h, y + k3);
        cout << "k1: " << k1 << endl;
        cout << "k2: " << k2 << endl;
        cout << "k3: " << k3 << endl;
        cout << "k4: " << k4 << endl;
        // следующий у
        y = y + (1.0 / 6.0)*(k1 + 2 * k2 + 2 * k3 + k4);
        cout << "y(" << x0+h << ")=" << y << endl;
        // Update next value of x
        x0 = x0 + h;
    }
 
    return y;
}
 
int main(){
    float x0, y0, x, h;
    cout << "Input x0: "; cin >> x0;
    cout << "Input x: "; cin >> x;
    cout << "Input y(x0)=y0: "; cin >> y0;
    cout << "Input h: "; cin >> h;
    cout << "method Runge-Kutta: " << rungeKutta(x0, y0, x, h) << endl;
 
    system("pause");
    return 0;
}