find : alpha4;see before (verion 1.0)
peter.vlasschaert@gmail.com,09/10/2016

definition : log to base 10

log10(x):= log(x)/log(10.0);

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

p1:'product(k, k, 1, 5)=product(k, k, 1, 5);
p2:'sum(k,k,1,5)=sum(k,k,1,5);

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));

alpha4:alpha(h,4,4);

ka_edta[1]:1.02e-2$
ka_edta[2]:2.14e-3$
ka_edta[3]:6.92e-7$
ka_edta[4]:5.50e-11$

alpha4:(alpha(h,4,4));

alpha4:alpha(h,4,4);

 pH =10,Y(m-) ,m=4,see :truebasic graph

alpha(10^(-10),4,4);

float(4.270785644654524*10^20*product(ka_edta[j],j,1,4)), numer;

for n:1 thru 5 do print(n);

for n:1 thru 2
 do print(n);

conch(pH):=float(10^(-pH));

pH:10;
h:conch(pH);
float(alpha(h,4,4));

ratsimp(4.270785644654524*10^20*product(ka_edta[j],j,1.0,4.0))$

float(427078564465452384256*ka_edta[1]*ka_edta[2]*ka_edta[3]*ka_edta[4]), numer;

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

float(1203478442970649985024*product(ka_edta[j],j,1,4)), numer;

 pH =10,Y(m-) ,m=4,see :truebasic graph

float(427078564465452384256*product(ka_edta[j],j,1,4)), numer;