Метод Рунге-Кутта

    float func= 0, F = 0, P = 0, U = 0, t = 0, e = 0, lambda = 0, gamma = 0; // Переменные для функции
    float Runge, K1 = 0, K2 = 0, K3 = 0, K4 = 0; //Переменные для Рунге-Кутта 4го порядка
 
    cout << "Insert p: ";
    cin >> P;
    cout << "\n";
    cout << "Insert U: ";
    cin >> U;
    cout << "\n";
    cout << "Insert F: ";
    cin >> F;
    cout << "\n";
    cout << "Insert e: ";
    cin >> e;
    cout << "\n";
    cout << "Insert gamma: ";
    cin >> gamma;
    cout << "\n";
    cout << "Insert lmbda: ";
    cin >> lambda;
    cout << "\n";
    cout << "Insert t: ";
    cin >> t;
    cout << "\n";
 
    func = (F*(sqrt(P/U))) * (sin(gamma)*cos(lambda) + ((e*pow(cos(gamma), 2) + 2*cos(gamma) + e)/(1 + cos(gamma)))*sin(lambda)); // Функция
 
    K1 = func + t*func;
    cout << "Answer for K1: " << K1 << endl; //Элементы из формулы Рунге-Кутта, их вычисление
    K2 = func + 0.5*K1;
    cout << "Answer for K2: " << K2 << endl;
    K3 = func + 0.5*K2;
    cout << "Answer for K3: " << K3 << endl;
    K4 = func + K3;
    cout << "Answer for K4: " << K4 << endl;
 
    Runge = func + (K1 + K2 + K3 + K4)/6; // Вычисление Рунге-Кутта
    cout << "Answer for Runge-Kutt(4th): ";
    cout << Runge << endl;
 
    system("pause");
    return 0;
}

...........................................................
cout << "Insert t: ";
    cin >> t;
    for(int i = 0; i < 10; i++)
    {       
        t = t + 70;
        cout << "\n";
 
        func = (F*(sqrt(P/U))) * (sin(gamma)*cos(lambda) + ((e*pow(cos(gamma), 2) + 2*cos(gamma) + e)/(1 + cos(gamma)))*sin(lambda)); // Функция
 
        K1 = func + t*func;
        cout << "Answer for K1: " << K1 << endl; //Элементы из формулы Рунге-Кутта, их вычисление
        K2 = func + 0.5*K1;
        cout << "Answer for K2: " << K2 << endl;
        K3 = func + 0.5*K2;
        cout << "Answer for K3: " << K3 << endl;
        K4 = func + K3;
        cout << "Answer for K4: " << K4 << endl;
 
        Runge = func + (K1 + K2 + K3 + K4)/6; // Вычисление Рунге-Кутта
        cout << "Answer for Runge-Kutt(4th): ";
        cout << Runge << endl;
    }
.........................................................

float step = 0, t;
    float x1;
    float x;
 
    cout << "Insert Integrator's step: ";
    cin >> t;
    cout << endl;
    cout << "Insert x: ";
    cin >> x;
    cout << endl;
 
    for(int i = 0; i < 10; i++)
    {       
        step = step + t;
        cout << "\n";
        float f(float x);
        {
            x1 = 2*x + 4;
        }
 
        K1 = step*f(x1);
        K2 = step*f(x1+K1/2);
        K3 = step*f(x1+K2/2);
        K4 = step*f(x1 + K3);
 
        Runge = f(x) + (1/6)*(K1 + 2*K2 + 2*K3 + K4); // Вычисление Рунге-Кутта
        cout << "Answer for Runge-Kutt(4th): ";
        cout << Runge << endl;
    }
    system("pause");
    return 0;