DataMuseum.dk

top - metrics - download

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

└─⟦621cfb9a2⟧ Bits:30002817 RC8000 Dump tape fra HCØ.  Detaljer om "HC8000" projekt.
    └─⟦0364f57e3⟧ 
        └─⟦this⟧ »iontxt« 

TextFile

opløselighed med ionstyrke
begin
  real L,s,s1,I,LNL,f,I0;
  integer ip,in,LP,LN;
  L:=readr(<:opl.produkt:>);
OM:
  LP:=readi(<:positiv ions ladning:>);
  LN:=readi(<:negativ  -      -:>);
I0:=readr(<:evt. fremmed ionstyrke:>);
  if LN<0 then LN:=-LN;
  if LP<=0 or LN=0 then
  begin
    write1(out,<:illegal ladning:>);
    goto OM
  end;
  if LN mod LP = 0 then
  begin
    in:=1; ip:=LN//LP; 
  end
  else
  if LP mod LN=0 then
  begin
    ip:=1; in:=LP//LN
  end 
  else
  begin
  ip:=LN;
  in:=LP; ip:=LN;
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*LN*LN)/(ip+in);
  s:=0;
K:
  I:=.5*(LP**2*ip+LN**2*in)*s+I0;
  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◀