mirror of
https://gitlab.com/freepascal.org/fpc/source.git
synced 2025-04-05 14:48:18 +02:00
63 lines
1.4 KiB
ObjectPascal
63 lines
1.4 KiB
ObjectPascal
program detgpdte;
|
|
|
|
uses
|
|
typ,
|
|
iom,
|
|
det;
|
|
|
|
const
|
|
n1 = -5;
|
|
n2 = 10;
|
|
rwa = n2 - n1 + 1;
|
|
var
|
|
e, t, aantal, i, j, k, l, n, term: ArbInt;
|
|
d: ArbFloat;
|
|
a: array[n1..n2, n1..n2] of ArbFloat;
|
|
begin
|
|
iom.npos := 1000; {max. width of output to 1000, since this is piped}
|
|
Write(' program results detgpdte');
|
|
case sizeof(ArbFloat) of
|
|
4: writeln('(single)');
|
|
6: writeln('(real)');
|
|
8: writeln('(double)');
|
|
10: writeln('(Extended)');
|
|
end;
|
|
Read(aantal);
|
|
writeln;
|
|
writeln(' number of examples : ', aantal: 3);
|
|
for t := 1 to aantal do
|
|
begin
|
|
writeln;
|
|
writeln(' example nr ', t: 3);
|
|
Read(k, l, n);
|
|
for i := 1 to n do
|
|
for j := 1 to i do
|
|
Read(a[k + i - 1, l + j - 1]);
|
|
detgpd(n, rwa, a[k, l], d, e, term);
|
|
writeln;
|
|
writeln(' A =');
|
|
for i := 1 to n do
|
|
for j := 1 to i - 1 do
|
|
a[k + j - 1, l + i - 1] := a[k + i - 1, l + j - 1];
|
|
iomwrm(output, a[k, l], n, n, n2 - n1 + 1, numdig);
|
|
writeln;
|
|
writeln('term=', term: 2);
|
|
writeln;
|
|
if term = 3 then
|
|
writeln(' wrong input')
|
|
else
|
|
if term = 2 then
|
|
writeln(' matrix not pos-def')
|
|
else
|
|
begin
|
|
Write(' det =', d: numdig);
|
|
if e <> 0 then
|
|
Write(' * 8**', e: 3);
|
|
writeln;
|
|
end; {term=1}
|
|
writeln('------------------------------------------------------');
|
|
end; {t}
|
|
Close(input);
|
|
Close(output);
|
|
end.
|