comment programma t.b.v. cursus WR, R 805, MLP 160762 / 2003; begin integer n; n:= read; begin integer i, j, k; real p1, p2, N; real array v1, v2, w1, w2[1: n], A[1: n, 1: n]; procedure iteratie; begin real procedure SCPR (a, b); array a, b; begin SCPR:= INPROD (j, 1, n, a[j], b[j]) end; for i:= 1 step 1 until n do begin w1[i]:= INPROD (j, 1, n, A[i, j], v1[j]); w2[i]:= INPROD (j, 1, n, A[i, j], v2[j]) end; N:= SCPR (v1, v1) + SCPR (v2, v2); p1:= (SCPR (v1, w1) + SCPR (v2, w2)) / N; p2:= (SCPR (v1, w2) - SCPR (v2, w1)) / N; N:= p1 * p1 + p2 * p2; for i:= 1 step 1 until n do begin v1[i]:= (p1 * w1[i] + p2 * w2[i]) / N; v2[i]:= (p1 * w2[i] - p2 * w1[i]) / N end end iteratie; for i:= 1 step 1 until n do begin v1[i]:= v2[i]:= 0; for j:= 1 step 1 until n do A[i, j]:= read end; v1[1]:= v2[2]:= 1; NLCR; NLCR; for k:= 1 step 1 until 20 do begin NLCR; NLCR; iteratie; FIXT (3, 0, k); FLOT (5, p1); FLOT (5, p2); NLCR; for i:= 1 step 1 until n do FLOT (5, v1[i]); NLCR; for i:= 1 step 1 until n do FLOT (5, v2[i]) end end end