Теория, планы, алгоритм и компьютерные программы при пяти уровнях факторов. Черный А.А. - 61 стр.

  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;
  G3:=S5+I0*R7+M0*N9;
  G4:=T9*N0-T8;
  G5:=Z9-T9*S0-G4*P0;
  G6:=T9*R0-T7+G4*U0;
  G7:=G4*Z0+T9*W0-G3;
  G0:=(Z6*G7+Z8*G6)/(Z6*G5-Z7*G6);
  H0:=(G0*Z7+Z8)/Z6;
  K0:=G0*P0+H0*U0+Z0;
  L0:=-(W0+G0*S0+H0*R0+K0*N0);
 end;{GB3860}
procedure GB4290;
 begin
  for J:=1 to X do
   begin
     I[J]:=ste(F[J],J1)+V1;
     K[J]:=ste(F[J],O1)+U1*ste(F[J],J1)+Q1;
     L[J]:=ste(F[J],P1)+I1*ste(F[J],O1)+M1*ste(F[J],J1)+F1;
     M[J]:=ste(F[J],T1)+G1*ste(F[J],P1)+H1*ste(F[J],O1)
     +K1*ste(F[J],J1)+L1;
   end;
 end;{GB4290}