#include <math.h>
#include <stdio.h>
#include <iostream>
using namespace std;

int main()
{

    double mpod = 42.2, dmpod = 0.2, mp = 15.8, dmp = 0.3, ms = 58, dMp2 = 1, g = 9.8, dg = 0.05, D = 45, dD = 0.1, s = 20, ds = 0.2, d1 = 48, dd1 = 0.1, x = 500, dx = 1, n1 = 1, n2 = 5, n3 = 9;
    double a, b, y, J, m0, t11 = 3.99, t12 = 3.4, t13 = 3.01, e1, e2, e3, m1, m2, m3, J0, dJ0, t21 = 3.99, t22 = 4.00, t23 = 4.01, t24 = 3.99, t25 = 4.00;
    double  dJ0J02, Jter, Jpr, Jt32, Jt33, Jt34, Jp32, Jp33, Jp34, t32 = 4.25, t33 = 5.27, t34 = 8.52, r12, r22, r32;
    double summE, summM, summM2, summEM;
    // tabl1
    e1 = round(400 * x / (D*t11*t11)) / 100;
    e2 = round(400 * x / (D*t12*t12)) / 100;
    e3 = round(400 * x / (D*t13*t13)) / 100;
    m1 = round(100 * mpod) / 100;
    m2 = round(100 * (mpod + mp)) / 100;
    m3 = round(100 * (mpod + mp + mp)) / 100;
    summE = e1 + e2 + e3;
    summM = m1 + m2 + m3;
    summM2 = m1 * m1 + m2 * m2 + m3 * m3;
    summEM = e1 * m1 + e2 * m2 + e3 * m3;
    a = (3 * summEM - summE * summM) / (3 * summM2 - summM * summM);
    b = (summE - a * summM) / 3;

    m0 = round((-1*(b / a)) * 1000000000000000000 / ((((-1 * (b / a)) * 1000000000000000000) - 2))) / ((20 - ((-1 * (b / a)) * 1000000000000000000)));
    cout << "e1 = " << e1 << '\t' << "e2 = " << e2 << '\t' << "e3 = " << e3 << '\n';
    cout << "m1 = " << m1 << '\t' << "m2 = " << m2 << '\t' << "m3 = " << m3 << '\n';

    m0=0.21714285714286541278911564626008;

    cout << m0 << '\n';
    cout << '\n';

    cout << '\n';
    // tabl 2
    double taver = (t21 + t22 + t23 + t24 + t25) / 5;
    double dt = 2.8*sqrt(((taver - t21)*(taver - t21) + (taver - t22)*(taver - t22) + (taver - t23)*(taver - t23) + (taver - t24)*(taver - t24) + (taver - t25)*(taver - t25)) / 4) / sqrt(5);
    cout << "taver = " << taver << "   dt =   " << dt;
    double rtaver, rdt;
    cout << '\n'<<"rdt= ";
    cin  >> rdt;
    cout << '\t' << "rtaver = "; cin >> rtaver;
    cout << '\n';
    J0 = 1000 * g*D*D*1e-6*(mpod - m0)*rtaver*rtaver / (8 * x);
    dJ0J02 = sqrt((dg / g)*(dg / g) + (2 * dD / D)*(2 * dD / D) + (dmpod / (mpod - m0))*(dmpod / (mpod - m0)) + (2 * rdt / rtaver)*(2 * rdt / rtaver) + (dx / x)*(dx / x));
    dJ0 = dJ0J02 * J0;
    cout << "J0 =  " << J0 << "   dJ0 =    " << dJ0;
    cout << '\n'<< "J0 = ";
    cin >> J0;
    cout << '\t' << "dJ0 = "; cin >> dJ0;
    cout << '\n';
    // tabl 3
    r12 = round(0.1*(0.5*(d1 + s * (2 * n1 + 1)))*(0.5*(d1 + s * (2 * n1 + 1)))) / 100;
    r22 = round(0.1*(0.5*(d1 + s * (2 * n2 + 1)))*(0.5*(d1 + s * (2 * n2 + 1)))) / 100;
    r32 = round(0.1*(0.5*(d1 + s * (2 * n3 + 1)))*(0.5*(d1 + s * (2 * n3 + 1)))) / 100;
    cout << "r12 = " << r12 << '\t' << "r22 = " << r22 << '\t' << "r32 = " << r32 << '\n';
    double Jt31 = J0, Jp31 = J0;
    Jt32 = round(100 * (J0 + 4 * ms*r12*0.001)) / 100;
    Jt33 = round(100 * (J0 + 4 * ms*r22*0.001)) / 100;
    Jt34 = round(100 * (J0 + 4 * ms*r32*0.001)) / 100;
    Jp32 = round(100000 * g*D*D*0.000001*(mpod - m0)*t32*t32 / (8 * x)) / 100;
    Jp33 = round(100000 * g*D*D*0.000001*(mpod - m0)*t33*t33 / (8 * x)) / 100;
    Jp34 = round(100000 * g*D*D*0.000001*(mpod - m0)*t34*t34 / (8 * x)) / 100;
    // vivodi
    cout << "Jt32 = " << Jt32 << '\t' << "Jt33 = " << Jt33 << '\t' << "Jt34 = " << Jt34 << '\n';
    cout << "Jp32 = " << Jp32 << '\t' << "Jp33 = " << Jp33 << '\t' << "Jp34 = " << Jp34 << '\n';
    getchar();
}