DataMuseum.dk

Presents historical artifacts from the history of:

RC4000/8000/9000

This is an automatic "excavation" of a thematic subset of
artifacts from Datamuseum.dk's BitArchive.

See our Wiki for more about RC4000/8000/9000

Excavated with: AutoArchaeologist - Free & Open Source Software.


top - metrics - download

⟦3cb75c1a1⟧ TextFile

    Length: 2304 (0x900)
    Types: TextFile
    Names: »albrint1«

Derivation

└─⟦667bb35d6⟧ Bits:30007480 RC8000 Dump tape fra HCØ.
    └─⟦4334b4c0b⟧ 
        └─⟦this⟧ »albrint1« 

TextFile

beregning af ligevaegtsdata for hoejtryksforsoeg
begin
real nd,tk,tv,vmpd,vmpetv,vmpetk;
integer h,i,j,l,m,n;
array omraade(1:2);
zone ud(128,1,stderror);
day;
write(out,<:omraade   :>);
setposition(out,0,0);
readstring(in,omraade,1);
j:=1;
open(ud,4,string omraade(increase(j)),0);
read(ud,n);
 begin
 real array pd,pe,na(0:n);
 read(ud,tk,tv);
 vmpetv:=100000000;
 nd:=pe(0):=0;
 for i:=1 step 1 until n do
  begin
  read(ud,pd(i),pe(i));
  pd(i):=.05*pd(i);
  pe(i):=.05*pe(i);
vmpd:=83.143*tv/(2*pd(i))+.5*sqrt((83.143*tv/pd(i))**2+4*14.4*
  83.143*tv/pd(i));
  nd:=nd+6.15/vmpd-6.15/vmpetv;
  vmpetv:=83.143*tv/(2*pe(i))+.5*sqrt((83.143*tv/pe(i))**2+57.6*
  83.143*tv/pe(i));
  vmpetk:=83.143*tk/(2*pe(i))+.5*sqrt((83.143*tk/pe(i))**2+57.6*
83.143*tk/pe(i));
na(i):=nd-16.08/vmpetv;
  end;
 write(out,<:tabel   :>);
 setposition(out,0,0);
 readchar(in,m);
 if m=106 then
  begin
  write(out,<:foerste og sidste punkt (max. :>,n,<:)   :>);
  setposition(out,0,0);
  read(in,h,l);
  write(out,<: n     p bar       mol abs<10>:>);
  for i:=h step 1 until l do
  write(out,<<dd>,i,<:     :>,<<dd.dd>,pe(i),<:     :>,
<<d.ddd'd>,na(i),<:<10>:>);
  setposition(out,0,0);
  end;
 write(out,<:plot   :>);
 setposition(out,0,0);
 readchar(in,m);
 if m=106 then
  begin
  write(out,<:papir   :>);
  setposition(out,0,0);
setposition(in,0,0);
readchar(in,m);
  if m=106 then setplotname(<:calcm:>,0);
  write(out,<:foerste og sidste punkt (max :>,n,<:)   :>);
setposition(out,0,0);
  read(in,h,l);
  if m=106 then plotform(2,4,1) else
  plotform(0,10,10);
  plotadmini(na(h),na(l),pe(h),pe(l),0);
  for i:=h step 1 until l do
  plotpoint(na(i),pe(i),2);
  plotclose;
  end;
 end;
end
▶EOF◀