DataMuseum.dk

top - metrics - download

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

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

TextFile

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