#include <stdio.h>
#include <conio.h>
#include < math.h >
using namespace std;
int P1, P2, k0, E;//îáúÿâëåíèå öåëî÷èñëåííûõ ïåðåìåííûõ
float b, I, c1, c2, L, x, Y, Mx, Qx, fi, A11, A12, A21, A22, yy, fi0, beta, B1, B2;//îáúÿâëåíèå âåùåñòâåííûõ ïåðåìåííûõ
double Cosh(float x)//ðàñ÷åò øèíóñà
{
    double Cos;
    Cos = (exp(x*beta) + exp(-x*beta)) / 2;
    return Cos;
}
double Sinh(float x)//ðàñ÷åò ÷îñèíóñà
{
    double Sin;
    Sin = (exp(x*beta) - exp(-x*beta)) / 2;
    return Sin;
}
double Y1(float x)//ïåðâàÿ ôóíêöèÿ Êðûëîâà
{
    double one;
    one = Cosh(beta * x) * cos(beta * x);
    return one;
}
double Y2(float x)//âòîðàÿ ôóíêöèÿ Êðûëîâà
{
    double two;
    two = 0.5 * (Cosh(x) * sin(beta * x) + Sinh(x) * cos(beta * x));
    return two;
}
double Y3(float x)//òðåòüÿ ôóíêöèÿ Êðûëîâà
{
    double three;
    three = 0.5 * Sinh(x) * sin(beta * x);
    return three;
}
double Y4(float x)//÷åòâåðòàÿ ôóíêöèÿ Êðûëîâà
{
    double four;
    four = 0.25 * (Cosh(x) * sin(beta * x) - Sinh(x) * cos(beta * x));
    return four;
}
int main()
{
    L = 1.3;//Çàäàíèå çíà÷åíèÿ ïåðåìåííîé L
    P1 = 3800000; //Çàäàíèå çíà÷åíèÿ ïåðåìåííîé P1
    P2 = 7400000; //Çàäàíèå çíà÷åíèÿ ïåðåìåííîé P2
    b = 0.11; //Çàäàíèå çíà÷åíèÿ ïåðåìåííîé b
    k0 = 500000000; //Çàäàíèå çíà÷åíèÿ ïåðåìåííîé K0
    E = 20000000; //Çàäàíèå çíà÷åíèÿ ïåðåìåííîé E
    I = 0.346; //Çàäàíèå çíà÷åíèÿ ïåðåìåííîé I
    c1 = 0.7; //Çàäàíèå çíà÷åíèÿ ïåðåìåííîé c1
    c2 = 0.9; //Çàäàíèå çíà÷åíèÿ ïåðåìåííîé c2;
    x = 0;
    beta = pow((k0 * b / (4 * E * I)), 0.25);
    A11 = Y1(L); //Ýëåìåíò ìàòðèöû A11
    A12 = Y2(L) / beta; //beta 'Ýëåìåíò ìàòðèöû A12
    A21 = -4 * beta * Y4(L); //Ýëåìåíò ìàòðèöû A21
    A22 = Y1(L); //Ýëåìåíò ìàòðèöû A22
    B1 = -(P1 * Y4(L - c1) + P2 * Y4(L - c2)) / E / I / pow(beta, 3); //Ýëåìåíò ìàòðèöû B11
    B2 = -(P1 * Y3(L - c1) + P2 * Y3(L - c2)) / E / I / pow(beta, 2); //Ýëåìåíò ìàòðèöû B21
    yy = (B1 * A22 - B2 * A12) / (A11 * A22 - A21 * A12); //Âû÷èñëåíèå ïàðàìåòðà y0
    fi0 = (A11 * B2 - A21 * B1) / (A11 * A22 - A21 * A12); //Âû÷èñëåíèå ïàðàìåòðà fi0
    printf("|  x   |    y    |    fi    |        Q       |        M       |\n");
    while (x <= L)
    {
        Y = yy* Y1(x) + fi0 * Y2(x) / beta; //Íà÷àëüíîå çíà÷åíèå Y
        fi = fi0 * Y1(x) - 4 * yy * beta * Y4(x); //Íà÷àëüíîå çíà÷åíèå fi
        Mx = 4 * E * I * beta * (beta * yy * Y3(x) + fi0 * Y4(x)); //Íà÷àëüíîå çíà÷åíèå M
        Qx = 4 * E * I * pow(beta, 2) * (beta * yy * Y2(x) + fi0 * Y3(x)); //Íà÷àëüíîå çíà÷åíèå Q
        if (x<c1 && x<c2)//Óñëîâèå 1
        {
            Y = Y; //Çíà÷åíèå Y ïðè âûïîëíåíèè óñëîâèÿ 1
            fi = fi; //Çíà÷åíèå fi ïðè âûïîëíåíèè óñëîâèÿ 1
            Mx = Mx; //Çíà÷åíèå M ïðè âûïîëíåíèè óñëîâèÿ 1
            Qx = Qx; //Çíà÷åíèå Q ïðè âûïîëíåíèè óñëîâèÿ 1
        }
        if (x >= c1) //Óñëîâèå 2
        {
            Y = Y + P1 * Y4(x - c1) / E / I / pow(beta, 3); //Çíà÷åíèå Y ïðè âûïîëíåíèè óñëîâèÿ 2
            fi = fi + P1 * Y3(x - c1) / E / I / pow(beta, 2); //Çíà÷åíèå fi ïðè âûïîëíåíèè óñëîâèÿ 2
            Mx = Mx - P1 * Y2(x - c1) / beta; //Çíà÷åíèå M ïðè âûïîëíåíèè óñëîâèÿ 2
            Qx = Qx - P1 * Y1(x - c1); //Çíà÷åíèå Q ïðè âûïîëíåíèè óñëîâèÿ 2
        }
        if (x >= c2) //Óñëîâèå 3
        {
            Y = Y + P2 * Y4(x - c2) / E / I / pow(beta, 3); //Çíà÷åíèå Y ïðè âûïîëíåíèè óñëîâèÿ 3
            fi = fi + P2 * Y3(x - c2) / E / I / beta / beta; //Çíà÷åíèå fi ïðè âûïîëíåíèè óñëîâèÿ 3
            Mx = Mx - P2 * Y2(x - c2) / beta; //Çíà÷åíèå M ïðè âûïîëíåíèè óñëîâèÿ 2
            Qx = Qx - P2 * Y1(x - c2); //Çíà÷åíèå Q ïðè âûïîëíåíèè óñëîâèÿ 3
        }
        printf("| %-5.2f| %-8.5f| %-9.5f| %-15.5f| %-15.5f|\n", x, Y, fi, Mx, Qx);//Âûâîä
        x = x + 0.05;
    }
    getch();
}