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

⟦2008513e9⟧ TextFile

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

Derivation

└─⟦621cfb9a2⟧ Bits:30002817 RC8000 Dump tape fra HCØ.  Detaljer om "HC8000" projekt.
    └─⟦0364f57e3⟧ 
        └─⟦7e928b248⟧ »algbib« 
            └─⟦this⟧ 

TextFile


clear prsep
prsep=set 50
prsep=algol index.no 
\f


PRSEP   PAL  MAJ 1975
begin
comment program til beregning af PR- og QQ-separation

inputregler:
l)
2);
integer n;
real A,B,C,ro,kappa,T,deltaQQ,c,cf1,cf2,sepB;
c:=2.997925'10;cf1:=16.8575; cf2:=27.9908'40;


read(in,A,B,C,T); read(in,n);
A:=case n of(A, A*'6/c, cf1/A, cf2/A);
B:=case n of(B, B*'6/c, cf1/B, cf2/B);
C:=case n of(C, C*'6/c, cf1/C, cf2/C);
kappa:=(2*B-A-C)/(A-C); ro:=(A-C)/B;
deltaQQ:=ro*(0.93 - 0.63*abs(kappa)**1.3)*sqrt((A+B+C)*T);
write(out,<:delta(QQ) =  :>,<<dd.dd>,deltaQQ,<:  1/cm:>);
sepB:=10*sqrt(C*T/9);
if ro>1/3 and ro<=1/2 then deltaQQ:= (1/3)*ro*sepB*sqrt(2);
if ro>1/2 then deltaQQ:= (2/3)*(ro**2)*sepB*sqrt(2);
write(out,<:<10>:>,<:delta(QQ) =  :>,<<dd.dd>,deltaQQ,<:  1/cm:>);
end
▶EOF◀