Руководство к решению задач по механике материалов и конструкций. Егодуров Г.С - 274 стр.

               YY1:=Y2[1];
               End
       If (X1>DZ1) And (X1<=DZ2) THEN
               Begin
               D[1]:=Y2[2];
               D[2]:=(-PP*Y2[1]+R0*X1-CC1*(X1-DZ1)*YY1)/F2;
               YY2:=Y2[1];
               End;
       If (X1>DZ2) and (X1<=1) THEN
               Begin
               D[1]:=Y2[2];
               D[2]:=(-PP*Y2[1]+RO*X1-CC1*(X1-DZ1)*YY1-CC2*(X1-DZ2)*YY2)/F3;
               End
End; {FCT}
       {Процедура общего виды для интегрирования уравнения методом Рунге-Кутта}
Procedure RK(Var X0: real; Var Y0: Wect);
       Var
       I:integer;
       Y1,P0,P1,P2,P3:Wect;
Begin {RK}
       FCT(X0,Y0,P0);
       For I:=1 to 2 Y1[I]:=Y0[I]+H*P0[I]/2;
       FCT(X0+H/2,Y1,P1);
       For I:=1 to 2 do Y1[I]:=Y0[I]+H*P1[I]/2;
       FCT(X0+H/2, Y1, P2);
       For I:=1 to 2 do Y1[I]:=Y0[I]+H*P2[I];
       FCT(X0+H,Y1,P3);
       For I:=1 to 2 do 0[1]:=Y0[I]+H*(P0[I]+2*P1[I]+2*P2[I]+P3[I]/6;
       X0:=X0+H;
End; {RK}
                                           {Основная программа}
       BEGIN
       Clrscr;
       Writeln (`Введите DL1, DL2,DL в метрах`);
               Read (DL1,DL2,DL);
       Writeln(`Введите жесткость упругих опор C1,C2 в H/м`); Readln(C1,C2);
       Writeln(`Введите изгибную жесткость EI1,EI2,EI3 в H*м*м`);
       Readln(EI1,EI2,EI3);
       Writeln(`Введите начальное значение силы P в H`);
       Readln(P);
       Writeln(`Введите шаг приращения нагрузки DP вd H`);
       Readln(DP);
       Writeln(`Введите шаг интегрирования НН в метрах`);
       Readln(HH);
       CC1:=C1*DL*DL*DL/EI1; CC2:=DL*DL*DL/EI1; DZ1:=DL1/DL;
       DZ2:=DL2/DL; F2:=EI2/EI1; F3:=EI3/EI1; PP:=P*DL*DL/EI1;
       DD:=DP*DL*DL/EI1; H:=HH/DL; M:1;
2:     X:=0; Y[1]:=0; Y[2]:=1; R0:=0;
       Writeln X<1 do RK(X,Y);
       A1:=Y[1];
       A2:=Y[2]; X:=0; Y[1]=0; Y[2]:=0; R0:=0;
       While X<1 do RK(X,Y);



                                         274