DataMuseum.dk

top - metrics - download

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

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

TextFile

;ali time 5 0
mode list.yes
message tfcalc on.rydiv
clear fcalc
fcalc=set 70
permanent fcalc.17
fcalc=hcalg list.yes
1980-03-23
f-value calculation 2-channel case
begin
integer c,char,i,j;
real fpure,a,b,f0,r1,r2,r12,z1,z2,ang1,ang2;
if -,readifp(<:j:>,j) then j:=1;
read(in,ang1,ang2);
comment for S series in B I etc ang1=2 and ang2=2/3,
        for D series in B I etc ang1=4 and ang2=10/3,
See Sobelman p. 312 and p. 318;
write(out,nl,2,<:angular factors :>,ang1,ang2);
write(out,nl,2,<:f(pure):>,sp,2,<:f(calc)+:>,sp,1,<:f(calc)-:>,sp,3,
      <:r1:>,sp,6,<:r2:>,sp,6,<:z1:>,sp,6,<:z2:>);
for c:=readchar(in,char) while char<>25 and char<>101 do
begin
  repeatchar(in);
  read(in,fpure,r1,r2);
  r12:=r1*r1;
  f0:=fpure/r12;

  write(out,nl,1);
  for i:=1 step 1 until j  do
  begin
    read(in,z1,z2);
    a:=z1*r1;
    b:=ang2*z2/ang1*r2;
    write(out,nl,1,<< -d.dddd>,fpure,f0*(a+b)**2,f0*(a-b)**2,r1,r2,z1,z2);
  end;
  repeatchar(in);
  for c:=readchar(in,char) while c>6 and char<>25 do;
  repeatchar(in);
end;
end;
head 1
lookup b1sfval
if ok.no
contract from.rydiv b1sfval
fcalc b1sfval
head 1
lookup al1dfval
if ok.no
contract from.rydiv al1dfval
fcalc al1dfval j.1
mode list.no
▶EOF◀