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

Формирование матрицы С(i) и столбца ?i (про¬¬цедура VECT). ??? 2. Решение системы С(i) = ?i методом Гаусса (про¬¬¬цедура GAUSS). ??? 3. Определение собственных значений мат¬ри¬цы А (процедура LAMBDA). ??? 4. Определение собственных векторов мат¬ри¬цы А (процедура SVECT). Если предложенную схему принять за основу, то под¬программы VECT, LАMBDA и SVECT могут вы¬гля¬деть следующим образом. PROCEDURE VECT (CONST N : INTEGER; A:MAS1; VAR C: MAS1; VAR B : MAS11); VAR I,J,K : INTEGER; X : MAS; S : REAL; BEGIN FOR I := 1 TO N DO BEGIN X[I] := 0.0; C[I,N] := X[I]; C[N,I] := 0.0; END; X[1] := 1.0; C[1,N] := X[1]; FOR I := 2 TO N+1 DO BEGIN FOR J := 1 TO N DO BEGIN S := 0.0; FOR K := 1 TO N DO S := S + A[J,K]*X[K]; IF I<>N+1 THEN C[J,N+1-I] := S ELSE B[J] := S; END; FOR J := 1 TO N DO X[J] := C[J,N+1-I]; END; END. PROCEDURA LAMBDA (CONST N : INTEGER;E:REAL; QQ : MAS11; VAR XX : MAS11); LABEL CONT, CONT1; VAR I,J,K,NN : INTEGER; Q : ARRAY [0..30] OF REAL; T,C,P,Q1,D,U,V,F,W,H,Y,Z,L,R,S,X,M : REAL; A, B, BIG : REAL; DN : BOOLEAN; BEGIN K := 1; DN := FALSE; NN:= N; X := 0.0; BIG := 1.0E-10; FOR I := 1 TO N+1 DO Q[I] := QQ[I]; REPEAT T := 1.0; C := Q[2] / Q[1]; IF NN =1 THEN BEGIN P := -C; Q1 := 0; END ELSE BEGIN IF NN = 2 THEN H := C*C/4.0 - Q[3]/Q[1]; ELSE BEGIN M := 10; C := 4.0; D := 8.0; U := 4.0; F := 1.0; W := 2.0; T := 0.0; V := 8.0; CONT: IF ABS(M-10.0)<BIG THEN BEGIN P := C; M := 0.0; Q1 := D; C := U; D := V; U := P; V := Q1; Y := C; Z := D; F := -F; END; M := M + 1.0; H := 0.0; Q1 := Q[1]; P := Q[2] - C*Q1; L := Q1; CONT1: FOR J := 3 TO NN DO BEGIN R := P; P := Q[J] - C*R - D*Q1; Q1 := R; R := L; L := Q1 - C*R - H*D; H := R; END; Q1 := Q[NN+1] - D*Q1; S := L + C*R; IF ABS(T)<=BIG THEN BEGIN X := D*R; H := R*X + S*L; IF ABS(H)<=BIG THEN GOTO CONT; END; C := C + (P*S - Q1*R)/H; D := D + (P*X+Q1*L)/H; IF ABS (C-Y+D-Z) < BIG THEN BEGIN F := -W; EXIT; END; W := -F; H := C*C/4

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