|
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: 1536 (0x600) Types: TextFile Names: »algprsep«
└─⟦621cfb9a2⟧ Bits:30002817 RC8000 Dump tape fra HCØ. Detaljer om "HC8000" projekt. └─⟦0364f57e3⟧ └─⟦7e928b248⟧ »algbib« └─⟦this⟧
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◀