DataMuseum.dk

top - metrics - download

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

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

TextFile

opløselighed med ionstyrke
begin
  real L,s,s1,I,cnL,f;
  integer ip,in,cp,cn;
  L:=readr(<:opl.produkt:>);
OM:
  cp:=readi(<:positiv ions ladning:>);
  cn:=readi(<:negativ  -      -:>);
  if cn<0 then cn:=-cn;
  if cp<=0 or cn=0 then
  begin
    write1(out,<:illegal ladning:>);
    goto OM
  end;
  if cn mod cp=0 then
  begin
    in:=1; ip:=cn//cp 
  end
  else
  if cp mod cn=0 then
  begin
    ip:=1; in:=cp//cn
  end
  else
  begin
  ip:=cn; in:=cp;
  end;
  write1(out,<:sammensætning:  P:>,ip,<:  N:>,in,"nl",1);
  cnL:=ln(L/ip**ip/in**in)/(ip+in);
  f:=.509*ln10*(ip*cp*cp+in*cn*cn)/(ip+in);
  s:=(-0);
LO:
  I:=.5*(cp**2*ip+cn**2*in)*s;
  s1:=exp(cnL+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◀