%main driver for exhaust_compo
%Gasoline CHmOn;
m=1.87;n=0;
Lambda=[0.6:0.1:1.4];
[xco2,xh2o,xco,xh2,xo2,xn2]=exhaust_compo(m,n,Lambda);

h=plot(Lambda,[xco2,xh2o,xco,xh2,xo2],'-');
xlabel('\lambda');
ylabel('Exhaust Mole fractions');
legend('CO2','H2O','CO','H2','O2');