DataMuseum.dk

top - metrics - download

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

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

TextFile

opløselighed med ionstyrke
begin
  real L,s,s1,I,lnL,f;
  integer ip,in,lp,kn;
  L:=readr(<:opl.produkt:>);
OM:
  lp:=readi(<:positiv ions ladning:>);
  kn:=readi(<:negativ  -      -:>);
  if kn<0 then kn:=-kn;
  if lp<=0 or kn=0 then
  begin
    write1(out,<:illegal ladning:>);
    goto OM
  end;
  if kn mod lp=0 then
  begin
    in:=1; ip:=kn//lp 
  end
  else
  if lp mod kn=0 then
  begin
    ip:=1; in:=lp//kn
  end
  else
  begin
  ip:=kn; 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*kn*kn)/(ip+in);
  s:=0;
K:
  I:=.5*(lp**2*ip+kn**2*in)*s;
  s1:=exp(lnL+f*sqrt(I)/(1+sqrt(I)));
  if abs(s-s1)>.001*s then
  begin
    s:=s1;
    goto K
  end;
  write1(out,<:opløselighed:>,s,<: mol/l:>)
end
▶EOF◀