|
DataMuseum.dkPresents historical artifacts from the history of: RC4000/8000/9000 |
This is an automatic "excavation" of a thematic subset of
See our Wiki for more about RC4000/8000/9000 Excavated with: AutoArchaeologist - Free & Open Source Software. |
top - metrics - download
Length: 1536 (0x600)
Types: TextFile
Names: »chebycheftx«
└─⟦621cfb9a2⟧ Bits:30002817 RC8000 Dump tape fra HCØ. Detaljer om "HC8000" projekt.
└─⟦0364f57e3⟧
└─⟦787c125fb⟧ »adjprocfile«
└─⟦00964e8f7⟧ Bits:30007478 RC8000 Dump tape fra HCØ.
└─⟦b2ec5d50f⟧
└─⟦787c125fb⟧ »adjprocfile«
└─⟦this⟧
; proj proc * page 1 31 08 77, 11.05;
; chebychef_c
; ***********
if listing.yes
char 10 12 10
chebychef_c = set 1
chebychef_c = algol
external real procedure chebychef_c
___________________________________
_ (mode, t, g, arg, funct, coef);
value mode, g;
integer mode, t, g;
real arg, funct, coef;
begin
integer i, nd;
real sum, n, c_arg;
array f_a(0:g);
nd := 4*(g + 1);
n := nd/2;
_
comment argument computation;
_____________________________
if mode extract 1 = 1 then
for t := 0 step 1 until g do
arg :=cos((2*t + 1)/n*pi);
_
comment function and coeff;
___________________________
if (mode shift (-1)) extract 1 = 1 then
begin
sum := 0;
for t := 0 step 1 until g do
begin
f_a(t) := funct;
sum := f_a(t) + sum;
end;
t := 0; <*jensen dev*>
chebychef_c :=
coef := 2*sum/n;
comment * page ;
for t := 1 step 1 until g do
begin
sum := 0;
c_arg := pi/n;
for i := 0 step 1 until g do
sum := f_a(i)*cos((((2*i + 1)*t) mod nd)*c_arg)
_ + sum;
coef := 4*sum/n;
end t-loop;
end mode 2
else
chebychef_c := 0;
end chebychef_c;
end
if warning.yes
(mode 0.yes
message chebychef_c not ok
lookup chebychef_c)
▶EOF◀