%solve ODE 
[V,x] = ode23s('prob624fun',[0 200],[6;14;0;4;6;0]); 
plot(V,x(:,1),'o',V,x(:,2),'d',V,x(:,3),'s'); 
ylabel('Molar Flow Rates'); 
xlabel('Volume'); 
legend('Fa','Fb','Fc'); 