Теория и практика эффективного математического моделирования. Черный А.А. - 360 стр.

    for J:=1 to X do
      begin
        SS:=SS+Q[J]*V[J]*Y[J];
        O[23]:=O[23]+(Q[J]*V[J])*(Q[J]*V[J]);
      end;
    B[23]:=SS/O[23]; SS:=0; O[24]:=0;
    for J:=1 to X do
      begin
        SS:=SS+I[J]*Q[J]*V[J]*Y[J];
        O[24]:=O[24]+(I[J]*Q[J]*V[J])*(I[J]*Q[J]*V[J]);
      end;
    B[24]:=SS/O[24]; SS:=0; O[25]:=0;
    for J:=1 to X do
      begin
        SS:=SS+P[J]*K[J]*V[J]*Y[J];
        O[25]:=O[25]+(P[J]*K[J]*V[J])*(P[J]*K[J]*V[J]);
      end;
    B[25]:=SS/O[25]; SS:=0; O[26]:=0;
    for J:=1 to X do
      begin
        SS:=SS+U[J]*K[J]*Q[J]*Y[J];
        O[26]:=O[26]+(U[J]*K[J]*Q[J])*(U[J]*K[J]*Q[J]);
      end;
    B[26]:=SS/O[26]; SS:=0; O[27]:=0;
    for J:=1 to X do
      begin
        SS:=SS+K[J]*Q[J]*V[J]*Y[J];
        O[27]:=O[27]+(K[J]*Q[J]*V[J])*(K[J]*Q[J]*V[J]);
      end;
    B[27]:=SS/O[27];
  end;{GB2040}
procedure GB3660(A,B,E,N,R:real; var V0,U0,Q0:real);
  var L2,N0,R0,N3,N4,N5:real;
  begin
    N0:=(ste(A,N)+ste(B,N)+ste(E,N))/3;
    R0:=(ste(A,R)+ste(B,R)+ste(E,R))/3;
    L2:=2*N;
    N3:=(ste(A,L2)+ste(B,L2)+ste(E,L2))/3;
    N4:=N+R;
    N5:=(ste(A,N4)+ste(B,N4)+ste(E,N4))/3;
    V0:=-N0;
    U0:=(N0*R0-N5)/(N3-N0*N0);
    Q0:=-(R0+U0*N0);
  end;{GB3660}
procedure GB3710(A,B,C,D,N,R,S:real; var V0,U0,Q0,I0,M0,F0:real);
  var
    N0,R0,S0,L2,N3,K2,R3,N4,N5,N6,N7,R4,R5,P0,Z1,Z2,Z3:real;
  begin
    N0:=(ste(A,N)+ste(B,N)+ste(C,N)+ste(D,N))/4;
    R0:=(ste(A,R)+ste(B,R)+ste(C,R)+ste(D,R))/4;
    S0:=(ste(A,S)+ste(B,S)+ste(C,S)+ste(D,S))/4;
    L2:=2*N;
    N3:=(ste(A,L2)+ste(B,L2)+ste(C,L2)+ste(D,L2))/4;
    K2:=2*R;
    R3:=(ste(A,K2)+ste(B,K2)+ste(C,K2)+ste(D,K2))/4;
    N4:=N+R;
    N5:=(ste(A,N4)+ste(B,N4)+ste(C,N4)+ste(D,N4))/4;
    N6:=N+S;




                                  360