|
|
DataMuseum.dkPresents historical artifacts from the history of: RC4000/8000/9000 |
This is an automatic "excavation" of a thematic subset of
See our Wiki for more about RC4000/8000/9000 Excavated with: AutoArchaeologist - Free & Open Source Software. |
top - metrics - download
Length: 1536 (0x600)
Types: TextFile
Names: »js3«
└─⟦667bb35d6⟧ Bits:30007480 RC8000 Dump tape fra HCØ.
└─⟦4334b4c0b⟧
└─⟦this⟧ »js3«
opløselighed med ionstyrke
begin
real L,s,s1,I,LNL,f;
integer ip,in,lp,LN;
L:=readr(<:opl.produkt:>);
OM:
lp:=readi(<:positiv ions ladning:>);
LN:=readi(<:negativ - -:>);
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;
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;
LL:
I:=.5*(lp**2*ip+LN**2*in)*s;
s1:=exp(LNL+f*sqrt(I));
if abs(s-s1)>.001*s then
begin
s:=s1;
goto LL
end;
write1(out,<:opløselighed:>,s,<: mol/l:>)
end
write1(out,<:opløselighed:>,s,<: mol/l:>)
end
▶EOF◀