rootChemicals.m

function Eq = rootChemicals(x,k,A)
%Example usage: rootChemicals([x1,x2,x3,x4,x5,x6],[k1,k2,k3,k4,k5,k6,k7,k8,k9,k10],[A1,A2,A3,A4,A5,A6])

Eq(1) = A(1)-x(1)*(k(2)*k(5)*x(3)*(k(4)*k(9)+k(7)*x(4))*(k(1)+x(5))+k(1)*k(4)*k(7)*k(9)*(k(2)*(k(5)+x(4))+x(4)*x(5)))/(k(1)*k(2)*k(4)*k(5)*k(7)*k(9));

Eq(2) = A(2)-x(2)*(k(2)*k(6)*x(3)*(k(10)*k(4)+k(8)*x(4))*(k(1)+x(5))+k(1)*k(10)*k(4)*k(8)*(k(2)*(k(6)+x(4))+x(4)*x(5)))/(k(1)*k(10)*k(2)*k(4)*k(6)*k(8));

Eq(3) = A(3)-x(3)*(k(10)*k(4)*k(9)*(k(8)*(k(7)+x(1))+k(7)*x(2))+k(7)*k(8)*(k(10)*x(1)+k(9)*(k(10)+x(2)))*x(4))*(k(1)+x(5))/(k(1)*k(10)*k(4)*k(7)*k(8)*k(9));

Eq(4) = A(4)-x(4)/(k(1)*k(10)*k(2)*k(4)*k(5)*k(6)*k(9))*(k(10)*k(6)*x(1)*(k(2)*k(5)*x(3)*(k(1)+x(5))+k(1)*k(4)*k(9)*(k(2)+x(5)))+k(5)*k(9)*(k(10)*k(6)*(k(2)*x(3)*(k(1)+x(5))+k(1)*k(4)*(k(2)+x(5)))+x(2)*(k(2)*k(6)*x(3)*x(5)+k(1)*(k(2)*k(6)*x(3)+k(10)*k(4)*(k(2)+x(5))))));

Eq(5) = A(5)-x(5)*(((1+x(1)/k(5)+x(2)/k(6))*x(4))/k(2)+(k(10)*k(4)*k(9)*(k(8)*(k(7)+x(1))+k(7)*x(2))*x(3)+k(7)*k(8)*(k(10)*x(1)+k(9)*(k(10)+x(2)))*x(3)*x(4))/(k(1)*k(10)*k(4)*k(7)*k(8)*k(9))+(k(3)+x(6))/k(3));

Eq(6) = A(6)-x(6)*(1+x(5)/k(3));