(ad * bd + ay * by)/(bd * bd + by * by);. y: = (bd * ay - ad * by)/(bd * bd + by * by); ; KOR (a: complex; var c: complex); r, f: real;: = sqrt (ad * ad + ay * ay);. d: = sqrt (r) * sqrt ((ad/r +1) /2);. y: = sqrt (r) * sqrt ((1-ad/r)/2)
END; Tij; {знаходження матриці Т з розкладання (2.3.1) за формулами (2.3.2)}
var k, i, j: integer;
BEGIN (A [1,1], T [1,1]); {знаходження 1 елемента по 1 формулою з (2.3.2)} i: = 1 to n do {знаходження елементов1 рядка і головною діагоналі по 3 формулою з (2.3.2)}
for j: = 1 to n do (i = j) and (i <> j) then
begin.d: = 0; {допоміжні змінні}
s.y: = 0;
for k: = 1 to i-1 do
begin (T [k, i], T [k, i], p); (s, p, s);; (A [i, j], s, p); ( p, T [i, j]); i
begin.d: = 0;. y: = 0; k: = 1 to i-1 do (T [k, i], T [k, i], p); (s , p, s);; (A [i, j], s, p);
DIL (p, T [i, j], T [i, j]);
endi> j then {так як матриця Т - верхня трикутна, то всі елементи під головною діагоналлю - нульові} [i, j]. d: = 0;
T [i, j]. y: = 0;;;
PROCEDURE Yi {знаходження стовпця у за формулами (2.3.6) з k, i: integer; з нижньої трикутної системи (2.3.4)}
BEGIN (b [1], T [1, 1], y [1]); {знаходимо у1}
for i: = 2 to n do {знаходимо всі інші у}. d: = 0; {допоміжні змінні}
sy: = 0; k: = 1 to i-1 do (T [k, i], y [k], p); (s, p, s);; (b [ i], s, p); (p, T [i, i], y [i]);;
END; Xi {знаходження шуканого вектор-стовпця х k, i: integer; за формулами (2.3.7) з верхньої трикутної матриці (2.3.5)} (y [n], T [n , n], x [n]); {знаходження хn} i: = n-1 downto 1 do {знаходження інших х}. d: = 0; {допоміжні змінні}
sy: = 0; k: = i +1 to n do (T [i, k], x [k], p); (s, p, s);; (y [ i], s, p); (p, T [i, i], x [i]);;
END; PROVERKA; {підставляємо отримані значення у вихідну i, j: integer; систему і вважаємо стовпець вільних членів заново. Ці значення заносимо в допоміжний вектор-i: = 1 to n do стовпець вільних членів z. Порівнюємо цей стовпець з вихідним стовпцем b. Якщо вони співпали, [i]: = 0; то система вирішена правильно. } p align="justify"> for j: = 1 to n do [i]: = z [i] + a [i, j]. d * x [j]. d; ( b , i, = , z [i]: 7: 7);;;; {очистити екран} ( vvedite razmernost matritsy n = );
readln (n); {введення розмірності матриці}
write ( vvedite pogreshnost vychislenia e =
