program eiggg2te; uses typ, eig, iom; const m1 = -9; m2 = 5; n1 = -10; n2 = 8; n3 = -7; n4 = 6; var i, j, l, nex, k1, k2, i1, j1, i2, j2, n, term: ArbInt; a: array[m1..m2, n1..n2] of ArbFloat; b: array[m1..m2, n3..n4] of ArbFloat; lam: array[m1..m2] of ArbFloat; begin Write(' program results eiggg2te'); case sizeof(ArbFloat) of 4: writeln('(single)'); 6: writeln('(real)'); 8: writeln('(double)'); end; Read(nex); writeln; writeln('number of examples', nex: 2); writeln; for l := 1 to nex do begin writeln('example number', l: 2); writeln; Read(i1, j1, i2, j2, n, k1, k2); for i := 1 to n do for j := 1 to i do Read(a[i1 + i - 1, j1 + j - 1]); for i := 1 to n do for j := 1 to i do Read(b[i2 + i - 1, j2 + j - 1]); eiggg2(a[i1, j1], n, n2 - n1 + 1, k1, k2, b[i2, j2], n4 - n3 + 1, lam[i1 + k1 - 1], term); writeln; writeln('A='); writeln; for i := 1 to n do iomwrv(output, a[i1 + i - 1, j1], i, numdig); writeln; writeln('B='); writeln; for i := 1 to n do iomwrv(output, b[i2 + i - 1, j2], i, numdig); writeln; writeln('term=', term: 2); writeln; if term = 1 then begin writeln('lambda', k1: 2, ' t/m', k2: 2, ' = '); iomwrv(output, lam[i1 + k1 - 1], k2 - k1 + 1, numdig); end; writeln; writeln('-------------------------------------------'); end; Close(input); Close(output); end.