% This is the program where you define the initial conditions, 
% initial time, final time, and call the file where
% you wrote the differential equations.
% 
% The initial conditions are x(1) = 1 and x(2) = 0.5
% 
% The integration should be performed from t = 0 to t=1
%
%

x0 = [1;0.5];


[t,x]=ode15s('nonlinear',[0 1],x0);
plot(t,x);
title('solution of sample equations');
xlabel('Time (s)');
ylabel('Concentration (gmole/l)');
legend('C1','C2');










