Ткачев Ю.А. Открытые и закрытые числовые системы

64 {поменяли местами действительные числа} u2:=DiagR[k];DiagR[k]:=DiagR[j];DiagR[j]:=u2; {поменяли местами целые числа - номера проб} k:=k-1; if k>=1 then goto 10; 20: end; end;{конец BubblePas} Procedure Statisk(nf,m:integer;var X:Tcomp;var mx:real; var sigma:real;var asim:real); {расчёт основных статистик массивов компонентов: nf,m - число проб, компонентов; X - массив компонентов; mx - среднее компонентов; sigma - дисперсия компонентов; asim - коэфф. асимметрии распределения} var i,j:integer; X2,X3:TComp; sum1,sum2,sum3,m2,m3,skvo:real; begin {расчёт квадратов компонентов} for i:=1 to nf do X2[i]:=sqr(X[i]); {расчёт кубов компонентов} for i:=1 to nf do X3[i]:=X2[i]*X1[i]; {обнуление сумм компонентов,квадратов,кубов,парных прозведений и сумм для асимметрии} sum1:=0;sum2:=0;sum3:=0;asim:=0; {расчёт сумм компонентов, их квадратов} for i:=1 to nf do begin sum1:=sum1+X[i]; sum2:=sum2+X2[i]; sum3:=sum3+X3[i]; end; {расчёт среднeго, среднего квадратов, среднего кубов} mx:=sum1/nf; m2:=sum2/nf; m3:=sum3/nf; {расчёт суммы для асимметрии} for i:=1 to nf do asim:=asim+(X[i]-mx)*(X[i]- mx)*(X[i]-mx); {расчёт среднего для асимметрии}