# Вторая пачка часть 190

0; FOR K := 1 TO L DO H:= H + A[I,K]*A[I,K]; (; ЕСЛИ ПАРАМЕТР H СЛИШКОМ МАЛ, ЧТОБЫ ГАРАНТИРОВАТЬ ОР¬ТОГОНАЛЬНОСТЬ, ТО ПРЕОБРАЗОВАНИЕ ХАУСХОЛДЕРА НЕ ВЫ¬ПОЛ¬НЯЕТСЯ, ПРОИСХОДИТ ПЕРЕХОД К МЕТКЕ 10 И ВЫХОД ИЗ ПРО¬ЦЕДУРЫ ;) IF H <= TOL THEN BEGIN E[I] := 0.0; E2[I]:= 0.0; GOTO 10; END; E2[I] := H; F := A[I,I-1]; IF F>=0 THEN G :=-SQRT(H) ELSE G := SQRT(H); E[I] := G; H := H - F*G; A[I,I-1] := F-G; F := 0.0; FOR J := 1 TO L DO BEGIN G := 0; (; ФОРМИРОВАНИЕ ЭЛЕМЕНТОВ ВЕКТОРА А Х U ;) FOR K := 1 TO J DO G := G + A[J,K]*A[I,K]; FOR K:= J+1 TO L DO G := G+A[K,J]*A[I,K]; G := G / H; E[J] := G; F := F+G*A[I,J]; END; (; ВЫЧИСЛЕНИЕ ПАРАМЕТРА К ;) H := F/(H+H); (; ПРЕОБРАЗОВАНИЕ МАТРИЦЫ А ;) FOR J := 1 TO L DO BEGIN F := A[I,J]; G := E[J]-H*F; E[J] := G; FOR K := 1 TO J DO A[J,K]:= A[J,K]-F*E[K]-G*A[I,K]; END; 10: H := D[I]; D[I] := A[I,I]; A[I,I] := H; END; END; {$F+} PROCEDURE PRINTR (N:INTEGER; Y:MAS); (; ПРОЦЕДУРА ДЛЯ ВЫВОДА ОДНОМЕРНОГО МАССИВА Y В ФАЙЛ ВЫВОДА;) VAR J,K1,Z:INTEGER; BEGIN J := 1; K1 := 10; WRITELN (LST,' '); IF K1>N THEN K1 := N; REPEAT FOR Z:=J TO K1-1 DO WRITE (LST,Я[Z]:12:6); WRITELN (LST,Y[K1]:12:6); J := K1+1; INC(K1,10); IF K1>N THEN K1 := N; UNTIL J>K1; END; PROCEDURE CORR (NN,MM : INTEGER; X:MAS2; VAR R : MAS1); (; ПРОЦЕДУРА РАСЧЕТА КОРРЕЛЯЦИОННОЙ ФУНКЦИИ. ДЛЯ УДОБ¬СТ¬ВА В НЕЕ ВКЛЮЧЕНЫ ВСЕ ВСПОМОГАТЕЛЬНЫЕ ПРОЦЕДУРЫ, ИС¬ПОЛЬ¬ЗУ¬ЕМЫЕ В ВЫЧИСЛЕНИЯХ. ;) FUNCTION SUMXR (NN,J,K: INTEGER; X:MAS2) : REAL; VAR SS : REAL; I : INTEGER; BEGIN SS:= 0.0; FOR I := 1 TO NN DO CASE K OF 1: SS := SS + X[I,J]; 2: SS := SS + SQR (X[I,J]); END; SUMXR := SS; END; FUNCTION SUMXYR (NN,J,K: INTEGER; X,Y:MAS2) : REAL; VAR SS : REAL; I : INTEGER; BEGIN SS:= 0.0; FOR I := 1 TO NN DO SS := SS + X[I,J]*Y[I,K]; SUMXYR := SS; END; VAR I, J, K : INTEGER; XIJ,XIK, XJK, XJ, XJ2, XI, XI2 : REAL; BEGIN FOR J := 1 TO NN DO BEGIN XJ := SUMXR (MM,J,1,X); XJ2 := SUMXR (MM,J,2,X); FOR K := 1 TO NN DO BEGIN XIJ := SUMXYR (MM,J,K,X,X); XI := SUMXR (MM,K,1,X); XI2 := SUMXR (MM,K,2,X); R[J,K] := (MM*XIJ - XI*XJ) / (SQRT (MM*XI2 - XI*XI)* SQRT (MM*XJ2 - XJ*XJ)); END; END; END; {$F+} PROCEDURE PRINTR2 (N1:INTEGER;A:MAS1;MINER:MAS3); (; ПРОЦЕДУРА ВЫВОДА ДВУМЕРНОГО МАССИВА А В ФАЙЛ ВЫ¬ВО¬ДА

Индекс
Элементарные функции Линейные уравнения Нелинейные уравнения Случайные числа