mirror of
https://gitlab.com/freepascal.org/fpc/source.git
synced 2025-04-05 18:07:56 +02:00
52 lines
1.0 KiB
ObjectPascal
52 lines
1.0 KiB
ObjectPascal
program eigbs1te;
|
|
|
|
uses
|
|
typ,
|
|
iom,
|
|
eig;
|
|
|
|
const
|
|
n1 = -100;
|
|
n2 = 100;
|
|
var
|
|
ex, nex, nel, p, q, n, b, term: ArbInt;
|
|
a: array[n1..n2] of ArbFloat;
|
|
lam: array[n1..n2] of ArbFloat;
|
|
begin
|
|
Write(' program results eigbs1te');
|
|
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 ex := 1 to nex do
|
|
begin
|
|
writeln('example number', ex: 2);
|
|
writeln;
|
|
Read(p, q, n, b);
|
|
nel := n * (b + 1) - (b * (b + 1)) div 2;
|
|
iomrev(input, a[p], nel);
|
|
eigbs1(a[p], n, b, lam[q], term);
|
|
writeln(' A = ');
|
|
iomwrv(output, a[p], nel, numdig);
|
|
writeln;
|
|
writeln('term=', term: 2);
|
|
if term = 1 then
|
|
begin
|
|
writeln;
|
|
writeln('lambda=');
|
|
iomwrv(output, lam[q], n, numdig);
|
|
end
|
|
else
|
|
writeln(' wrong input');
|
|
writeln;
|
|
writeln('-------------------------------------------');
|
|
end;
|
|
Close(input);
|
|
Close(output);
|
|
end.
|