|
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: »tdigamma«
└─⟦621cfb9a2⟧ Bits:30002817 RC8000 Dump tape fra HCØ. Detaljer om "HC8000" projekt.
└─⟦0364f57e3⟧
└─⟦e4d872f9f⟧ »cproc«
└─⟦00964e8f7⟧ Bits:30007478 RC8000 Dump tape fra HCØ.
└─⟦b2ec5d50f⟧
└─⟦e4d872f9f⟧ »cproc«
└─⟦this⟧
(mode list.yes
lookup tdigamma
clear user digamma
digamma=set 2
global digamma
mode 15.yes
digamma=algol
testdigamma=set 160
testdigamma=algol
testdigamma
mode 10.no 15.no list.no)
1979-1-4
external
real procedure digamma(x);
value x; real x;
if x<0 then digamma:=maxreal else
begin
integer i; real t,xm,term;
comment Abramowitz and Stegun p. 258-259;
xm:=x-1;
t:=-euler;
for i:=1,i+1 while abs term> '-9 do
begin
term:=xm/i/(i+xm);
t:=t+term;
end;
digamma:=t;
end digamma;
end;
begin
real x;
if lookupentry(<:digammalist:>)=0 then
begin
for x:=.1 step .1 until 4,4 do
write(out,nl,1,x,digamma(x));
end
else
begin
setplotname(<:tek603:>,0);
plotform(0,10,14);
plotadmini(0,4,-5,1.8,0);
plotgraph(x,digamma(x),.1,4,'-2);
end
end
▶EOF◀