DataMuseum.dk

top - metrics - download

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

└─⟦00964e8f7⟧ Bits:30007478 RC8000 Dump tape fra HCØ.
    └─⟦b2ec5d50f⟧ 
        └─⟦this⟧ »hnret« 

TextFile

opløselighed med ionstyrke
begin
  real L,s,s1,I,lnegL,f;
  integer ip,in,lp,lneg;
  L:=readr(<:opl.produkt:>);
OM:
  lp:=readi(<:positiv ions ladning:>);
  lneg:=readi(<:negativ  -      -:>);
  if lneg<0 then lneg:=-lneg;
  if lp<=0 or lneg=0 then
  begin
    write1(out,<:illegal ladning:>);
    goto OM
  end;
  if lneg mod lp=0 then
  begin
    in:=1; ip:=lneg//lp 
  end
  else
  if lp mod lneg=0 then
  begin
    ip:=1; in:=lp//lneg
  end
  else
  begin
  ip:=lneg; in:=lp;
end;
  write1(out,<:sammensætning:  P:>,ip,<:  N:>,in,"nl",1);
  lnegL:=ln(L/ip**ip/in**in)/(ip+in);
  f:=.509*ln10*(ip*lp*lp+in*lneg*lneg)/(ip+in);
  s:=0;
LO:
  I:=.5*(lp**2*ip+lneg**2*in)*s;
  s1:=exp(lnegL+f*sqrt(I)/(1+sqrt(I)));
  if abs(s-s1)>.001*s then
  begin
    s:=s1;
    goto LO
  end;
  write1(out,<:opløselighed:>,s,<: mol/l:>)
end
▶EOF◀