|
|
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: 768 (0x300)
Types: TextFile
Names: »tbisec«
└─⟦00964e8f7⟧ Bits:30007478 RC8000 Dump tape fra HCØ.
└─⟦b2ec5d50f⟧
└─⟦09b4e9619⟧ »thcømat«
└─⟦this⟧
\f
message bisec
bisec=algol message.no
external boolean procedure bisec(f,x,a,b,eps);
value a,b,eps; real f,x,a,b,eps;
begin
real q,x0;
x:=a; q:=f;
bisec:=true;
if abs(q)>0.0 then
begin
x:=b;
if f*q>0 then bisec:=false
else
begin
q:=(b-a)*(if q>0 then 1 else -1)*. 5;
x0:=x:=(a+b)*.5;
eps:=if eps<0 then 0 else eps*.5;
for q:=q*.5 while abs q +abs x0> eps+abs x0 do
x:=x0:=(if f>0 then q else -q)+x0;
end
end
end bisec
; end
▶EOF◀