

syms q1 q2

q=10;                          % [m3/s],outlet flow
c1=5;c2=15;c=14;               % [mol/m3], concentration       

eqn1 = q1+q2-q==0;             % mass labance
eqn2 = q1*c1+q2*c2 -q*c==0;    % component balance

eqn = [eqn1;eqn2];
var = [q1;q2];

sol = solve(eqn,var);        % find F1 and F2
q1  = sol.q1
q2  = sol.q2