find : alpha4;see before (verion 1.0)
[email protected],09/10/2016
definition : log to base 10
(%i1) | log10(x):= log(x)/log(10.0); |
(%i2) | log10(100.0); |
EDTA = H4(Y) , m=4
alpha4 = (ka(1)*ka(2)*ka(3)*ka(4))/
(h^4+ka(1)*h^3+ka(1)*ka(2)*h^2+ka(1)*ka(2)*ka(3)*h+..
+ka(1)*ka(2)*ka(3)*ka(4))
ka(1) = ka_edta(1)=1.02e-2
ka(2) = ka_edta(2)=2.14e-3
ka(3) = ka_edta(3)=6.92e-7
ka(4) = ka_edta(4)=5.50e-11
p1:1*2*3*4*5 = 120
p2:1+2+3+4+5 = 15
(%i4) |
p1:'product(k, k, 1, 5)=product(k, k, 1, 5); p2:'sum(k,k,1,5)=sum(k,k,1,5); |
(%i5) |
alpha(h,i,n):='product(ka_edta[j],j,1,i)/ (h^n+sum(h^(n-k)*product(ka_edta[j],j,1,k),k,1,n)); |
(%i6) | alpha4:alpha(h,4,4); |
(%i10) |
ka_edta[1]:1.02e-2$ ka_edta[2]:2.14e-3$ ka_edta[3]:6.92e-7$ ka_edta[4]:5.50e-11$ |
(%i11) | alpha4:(alpha(h,4,4)); |
(%i12) | alpha4:alpha(h,4,4); |
pH =10,Y(m-) ,m=4,see :truebasic graph
(%i13) | alpha(10^(-10),4,4); |
(%i14) | float(4.270785644654524*10^20*product(ka_edta[j],j,1,4)), numer; |
for n:1 thru 5 do print(n);
(%i15) |
for n:1 thru 2 do print(n); |
(%i16) | conch(pH):=float(10^(-pH)); |
(%i19) |
pH:10; h:conch(pH); float(alpha(h,4,4)); |
(%i20) | ratsimp(4.270785644654524*10^20*product(ka_edta[j],j,1.0,4.0))$ |
(%i21) | float(427078564465452384256*ka_edta[1]*ka_edta[2]*ka_edta[3]*ka_edta[4]), numer; |
(%i22) |
for pH:1 thru 14 do (hh:conch(pH),alpha(hh,4,4), vv:alpha(hh,4,4), print(pH,hh,vv)); |
pH =14,Y(m-) ,m=4,see :truebasic graph
(%i23) | float(1203478442970649985024*product(ka_edta[j],j,1,4)), numer; |
pH =10,Y(m-) ,m=4,see :truebasic graph
(%i24) | float(427078564465452384256*product(ka_edta[j],j,1,4)), numer; |