IC5_2

function IC5_2
K = 0.016;
ca0 = 41; %mol/L
cb0 = 29; %mol/L
cc0 = 1; %mol/L
R = @(E) (cc0 + E)./((ca0-2.*E)^2 .* (cb0 - E)) - K;
root1 = fzero(R, 1);
cc0 = 100; %mol/L
R = @(E) (cc0 + E)./((ca0-2.*E)^2 .* (cb0 - E)) - K;
root2 = fzero(R, 1);
cc0 = 1000; %mol/L
R = @(E) (cc0 + E)./((ca0-2.*E)^2 .* (cb0 - E)) - K;
root3 = fzero(R, 1);



end