Математическое моделирование при ортогонализации матриц. Черный А.А. - 69 стр.

  B[5]:=SS/O[5]; SS:=0; O[6]:=0;
  for J:=1 to X do
   begin
     SS:=SS+Q[J]*Y[J];
     O[6]:=O[6]+Q[J]*Q[J];
   end;
  B[6]:=SS/O[6]; SS:=0; O[7]:=0;
  for J:=1 to X do
   begin
     SS:=SS+I[J]*Q[J]*Y[J];
     O[7]:=O[7]+(I[J]*Q[J])*(I[J]*Q[J]);
   end;
  B[7]:=SS/O[7]; SS:=0; O[8]:=0;
  for J:=1 to X do
   begin
     SS:=SS+P[J]*K[J]*Y[J];
     O[8]:=O[8]+(P[J]*K[J])*(P[J]*K[J]);
   end;
  B[8]:=SS/O[8]; SS:=0; O[9]:=0;
  for J:=1 to X do
   begin
     SS:=SS+K[J]*Q[J]*Y[J];
     O[9]:=O[9]+(K[J]*Q[J])*(K[J]*Q[J]);
   end;
  B[9]:=SS/O[9];
 end;{GB1480}
procedure GB1600;
 begin
  SS:=0; O[10]:=0;
  for J:=1 to X do
   begin
     SS:=SS+U[J]*Y[J];
     O[10]:=O[10]+U[J]*U[J];
   end;
  B[10]:=SS/O[10]; SS:=0; O[11]:=0;
  for J:=1 to X do
   begin
     SS:=SS+I[J]*U[J]*Y[J];
     O[11]:=O[11]+(I[J]*U[J])*(I[J]*U[J]);
   end;
  B[11]:=SS/O[11]; SS:=0; O[12]:=0;
  for J:=1 to X do
   begin
     SS:=SS+K[J]*U[J]*Y[J];
     O[12]:=O[12]+K[J]*U[J]*K[J]*U[J];
   end;
  B[12]:=SS/O[12];
 end;{GB1600}
procedure GB1670;
 begin
  SS:=0; O[13]:=0;
  for J:=1 to X do
   begin
     SS:=SS+V[J]*Y[J];
     O[13]:=O[13]+V[J]*V[J];
   end;
  B[13]:=SS/O[13]; SS:=0; O[14]:=0;




                                             69