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

  For j:= 1 to M do Y1[i.j]:=Y[i,j]+H*(P0[i,j]+2*P1[i,j]+2*P2[i,j]+P3[i,j])/6;
  X:=X+H;
end;
{Процедура общего вида для расчета пластины методом начальных параметров}
Procedure PLAST;
Var
Y:Matr;
Z :Amatr;
R,C1,C2,MR,MT,W,TETA:real;
I,J.K,K1,IK,NP,KP,IP:integer;
         S:array[1..N] of real;
RZ,HZ,DZ,QRZ:array[1..50] of real;
Begin
Clrscr;
K:=0; K1:=0;
NUS(Y);
For i:=1 to NU do K1:=K1+KR[i];
K1:=K1+NU;
For NS:=1 to NU do
  Begin
      R:=RG[NS];
      H:=(RG[NS+1]-RG[NS])/NN[NS];
      NP:=NN[NS] div KR [NS];
      KP:=KR[NS]+1;
      For ip:=1 to KP do
        Begin
            K:=K+1;
            RZ[K]:=R;
            HZ[K]:=TOL(NS,R);
            DZ[K]:=DD(R);
            QRZ[K]:=Q(NS,R);
            For i:=1 to N do
               Begin
                       For j:=1 to M do Z[K,I,J]:=Y[I,J];
                       End;
               IF(IP<>KP) Then
                                For IK:=1 to NP do
                                       RKM (R,Y)
                                       Else break;
                       End;
               End;
            GUS(Y,C1,C2,Z);
            Writeln (`C1=`,C1,`C2=`,C2);
            Writeln;
            Writeln(` R(м)`,` H(м)`,`Q(H/м)`,`W(м)`,
               `       Mr(H*м/м)`,` Mt(H*м/м)`,` Teta(рад)`);
            Writeln;
            For i:=1 to K1 do
               Begin
                       For j:=1 to N do S[j]:=C1*Z[i,j,1]+Z[i,j,2];
                       W:=S[1]+C2;
                       MR:=S[3]/RZ[i];



                                     263