DataMuseum.dk

top - metrics - download

    Length: 1536 (0x600)
    Types: TextFile
    Names: »tlobachev«

└─⟦00964e8f7⟧ Bits:30007478 RC8000 Dump tape fra HCØ.
    └─⟦b2ec5d50f⟧ 
        └─⟦09b4e9619⟧ »thcømat« 
            └─⟦this⟧ 

TextFile

;lobachevtxt  rene moss 15-10-73/1710-73

lobachevski=algol
external
real procedure lobachevski(x);
value x; real x;
if x=0.0 then lobachevski:=0
else
begin
  boolean sign; real l;
  if x<0 then
  begin
    sign:=false; x:=-x-x
  end
  else
  begin
    sign:=true; x:=x+x
  end;
  if x>pi then system(9,0,<:<10>***call:>);
  if x<pi/2 then
  begin
    l:=x*x;
    l:=(((((((((((((
        0.147995167768'-16 *l+0.168953112864'-15)*l
       +0.195097545200'-14)*l+0.228324378368'-13)*l 
       +0.271489636400'-12)*l+0.329060828336'-11)*l
       +0.408350548176'-10)*l+0.522014514256' -9)*l
       +0.693549955712' -8)*l+0.970769199936' -7)*l
       +0.146398258376' -5)*l+0.248015873008' -4)*l
       +0.520833333328' -3)*l+0.208333333344' -1)*l*x
  end
  else
  begin
    x:=pi-x;
    l:=x*x;
    l:=(((((((
       -0.318631822176'-13 *l-0.169365068544'-11)*l
       -0.948943499392'-10)*l-0.574110817152' -8)*l
       -0.393675988896' -6)*l-0.347222222224' -4)*l
       -0.694444444448' -2)*l
       +(ln(x/2)-1)/2)*x+pi/2*ln(2)
  end;
  lobachevski:=(if sign then l else -l)
end lobachevski
; end
▶EOF◀