|
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: »fkursus5202«
└─⟦00964e8f7⟧ Bits:30007478 RC8000 Dump tape fra HCØ.
└─⟦b2ec5d50f⟧
└─⟦this⟧ »fkursus5202«
Bestemmelse af pH for n syre/base par ved bisection
begin
integer n,i;
real h,ha,hb,Kw,Kw2;
MORE:
n:=readi(<:antal syre/base par:>);
begin
array K,C(1:n);
real procedure sum(K,C,n,h);
value n,h;
real h; integer n; array K,C;
begin
real s;
integer i;
s:=0;
for i:=1 step 1 until n do
s:=K(i)*C(i)/(h+K(i)) +s;
sum:=s;
end sum;
repeat
write1(out,<:K1.....Kn,C1....Cn:>,n);
outendcur('nl');
i:=read(in,K,C);
until i=2*n;
write1(out,"sp",2,<:K(i):>,"sp",6,<:C(i):>,"nl",1,"sp",2,"-",4,"sp",6,
"-",4,"nl",1);
for i:=1 step 1 until n do
write1(out,<<d.dd'-dd>,K(i),"sp",2,<<d.dd'dd>,C(i),"nl",1);
Kw:='-14;Kw2:=sqrt(Kw);
ha:=0;
for i:=1 step 1 until n do
ha:=ha:=ha+C(i);
hb:=Kw2-ha;
ha:=Kw2+ha;
bisec(sum(K,C,n,h)-h,h,ha,hb,'-4);
write1(out,<:pH =:>,<<dd.ddd>,-ln(h)/ln(10));
if readb(<: MORE:>) then goto MORE;
end
end
▶EOF◀