Математическое моделирование. Черный А.А. - 231 стр.

procedure GB2020;
 begin
  SS:=0; O[5]:=0;
  for J:=1 to X do
    begin
     SS:=SS+M[J]*Y[J];
     O[5]:=O[5]+M[J]*M[J];
    end;
  B[5]:=SS/O[5];
 end;{GB2020}
procedure GB3860(A,B,C,D,E,N,R,S,W:real;
var V0,U0,Q0,I0,M0,F0,G0,H0,K0,L0:real);
 var
  N0,R0,S0,W0,L2,N3,K2,R3,M2,S3,N4,N5:real;
  N6,N7,N8,N9,R4,R5,R6,R7,S4,S5,P0:real;
  Z1,Z2,Z3,Z4,Z5,Z6,Z7,Z0,Z8,Z9,T7:real;
  T8,T9,G3,G4,G5,G6,G7:real;
 begin
  N0:=(ste(A,N)+ste(B,N)+ste(C,N)+ste(D,N)+ste(E,N))/5;
  R0:=(ste(A,R)+ste(B,R)+ste(C,R)+ste(D,R)+ste(E,R))/5;
  S0:=(ste(A,S)+ste(B,S)+ste(C,S)+ste(D,S)+ste(E,S))/5;
  W0:=(ste(A,W)+ste(B,W)+ste(C,W)+ste(D,W)+ste(E,W))/5;
  L2:=2*N;
  N3:=(ste(A,L2)+ste(B,L2)+ste(C,L2)+ste(D,L2)+ste(E,L2))/5;
  K2:=2*R;
  R3:=(ste(A,K2)+ste(B,K2)+ste(C,K2)+ste(D,K2)+ste(E,K2))/5;
  M2:=2*S;
  S3:=(ste(A,M2)+ste(B,M2)+ste(C,M2)+ste(D,M2)+ste(E,M2))/5;
  N4:=N+R;
  N5:=(ste(A,N4)+ste(B,N4)+ste(C,N4)+ste(D,N4)+ste(E,N4))/5;
  N6:=N+S;
  N7:=(ste(A,N6)+ste(B,N6)+ste(C,N6)+ste(D,N6)+ste(E,N6))/5;
  N8:=N+W;
  N9:=(ste(A,N8)+ste(B,N8)+ste(C,N8)+ste(D,N8)+ste(E,N8))/5;
  R4:=R+S;
  R5:=(ste(A,R4)+ste(B,R4)+ste(C,R4)+ste(D,R4)+ste(E,R4))/5;
  R6:=R+W;
  R7:=(ste(A,R6)+ste(B,R6)+ste(C,R6)+ste(D,R6)+ste(E,R6))/5;
  S4:=S+W;
  S5:=(ste(A,S4)+ste(B,S4)+ste(C,S4)+ste(D,S4)+ste(E,S4))/5;
  V0:=-N0;
  U0:=(N0*R0-N5)/(N3-N0*N0);
  Q0:=-(R0+U0*N0);
  P0:=(N0*S0-N7)/(N3-N0*N0);
  Z1:=R0*S0-R5+P0*(N0*R0-N5);
  Z2:=U0*(N0*S0-N7)+U0*P0*(N0*N0-N3);
  Z3:=R3-R0*R0+2*U0*(N5-N0*R0);
  I0:=(Z1+Z2)/(Z3+(N3-N0*N0)*U0*U0);
  M0:=I0*U0+P0;
  F0:=-(S0+I0*R0+M0*N0);
  Z4:=R0+U0*N0;
  Z5:=Z4*N0-N5-U0*N3;
  Z6:=R3+U0*N5-Z4*R0-Z5*U0;
  Z7:=Z4*S0+Z5*P0-R5-U0*N7;
  Z0:=(N0*W0-N9)/(N3-N0*N0);
  Z8:=Z5*Z0+Z4*W0-R7-U0*N9;
  Z9:=S3+I0*R5+M0*N7;
  T7:=R5+I0*R3+M0*N5;
  T8:=N7+I0*N5+M0*N3;
  T9:=S0+I0*R0+M0*N0;



                                            231