#include <iostream>
#include <cmath>
#define g           9.81        //ускорение свободного падения
#define rho         1.2041      //плотность воздуха при норм. давлении и температуре 20 С
#define C_D         1.0         //Drag_Coefficient  wikipedia: 1.0-1.3 для человека
#define A           1.75        //Body_Surface_Area
#define mu          0.015       //коэффициент трения скольжения для коньков
#define F_d_mult    rho*C_D*A/2 //
#define F_f_mult    mu*g        //


//v0 = 10 м/с
//rot_coef = 1.1 | 3.3 | 8.5 (1-|2-|3-ой)
double lab2(double m, double v0, double a0, double R, double rot_coef) {
    double dt = 1.0 / 60.0;
    double t = 0, x = 0, L = 5;
    double v = v0;
    //прямолинейное равноускоренное движение
    while (x < L) {
        t += dt;
        double a = (a0 - F_d_mult*v0*v0 / m - F_f_mult);
        v += a*t;
        x += v*t + a*t*t / 2;
    }
    //криволинейное с радиусом R
    double v_curv = v / R;
    //доворот
    return rot_coef*v_curv;
}

int main() {
    double m = 50;
    double rot_coef = 1.1;
    double R = 1;
    double a = 1;
    for (int i = 0; i < 10; ++i) {
        double v = i;
        std::cout << "v: " << v << "; "
            << "v_f: " << lab2(m, v, a, R, rot_coef) << std::endl;
    }
    getchar();

    return 0;
}