/*Do file for calculating chemical equilibrium concentrations*/
keep Ao Bo Co Do Keq
gen a=1-Keq
gen b=Co+Do+Keq*Ao+Keq*Bo
gen c=Co*Do-Keq*Ao*Bo
gen x=(-b+sqrt(b^2-4*a*c))/(2*a)
gen Aeq=Ao-x
gen Beq=Bo-x
gen Ceq=Co+x
gen Deq=Do+x
gen x_Ao=x/Ao