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

 procedure GB2000; procedure GB2020;
 procedure GB4290;
 procedure GB4340;
 procedure GB4420;
 procedure GB4690;
Implementation
function ste(Q,A:real):real;
{Функция возведения в степень}
 begin
  if A<0 then
     begin
       A:=abs(A);
       if Q=0 then Q:=1E-5;
       ste:=1/exp(A*ln(Q));
     end
    else
     begin
       if Q=0 then Q:=1E-5;
       if A=0 then ste:=1
            else ste:=exp(A*ln(Q));
     end;
 end;{ste}
 procedure GB1410;
  begin
    SS:=0; O[1]:=0;
    for J:=1 to X do
     begin
       SS:=SS+Y[J]; O[1]:=O[1]+1;
     end;
    B[1]:=SS/O[1];
    SS:=0; O[2]:=0;
    for J:=1 to X do
     begin
       SS:=SS+I[J]*Y[J]; O[2]:=O[2]+I[J]*I[J];
     end;
    B[2]:=SS/O[2];
    SS:=0; O[3]:=0;
    for J:=1 to X do
     begin
       SS:=SS+K[J]*Y[J]; O[3]:=O[3]+K[J]*K[J];
     end;
    B[3]:=SS/O[3];
  end;{GB1410}
 procedure GB1480;
  begin
    SS:=0; O[4]:=0;
    for J:=1 to X do
     begin
       SS:=SS+P[J]*Y[J];
       O[4]:=O[4]+P[J]*P[J];
     end;
    B[4]:=SS/O[4]; SS:=0; O[5]:=0;
    for J:=1 to X do
     begin
       SS:=SS+I[J]*P[J]*Y[J];
       O[5]:=O[5]+(I[J]*P[J])*(I[J]*P[J]);
     end;




                                                 68