%calcul
options=odeset('abstol', 1.0e-4, 'reltol', 1.0e-4);
[t y]=ode45('tricl', [0:0.001:0.7], [1 0 0 ], options);
q = y(:,1);
r = y(:,2);
s = y(:,3);
t = 1-y(:,1)-y(:,2)-y(:,3);
p = r+2*s+3*t;
[m i] = max(r);
%prsentation
delete(gcf);
figure('name','trichloration du benzne', 'number', 'off', 'Menu', 'none');
subplot('Position',[0.1 0.4 0.8 0.55]);
plot(p, q, 'y-',p, r, 'r:',p, s, 'b-.', p, t, 'g--')
legend('benzne', 'C_6H_5Cl', 'C_6H_5CL_2', 'C_6H_5CL_3');
subplot('Position',[0.1 0.01 0.8 0.35]);
axis off
st{1} = sprintf('Le maximum de C_6H_5Cl vaut %6.4f et s''obtient pour',m);
st{2} = sprintf('une concentration de C_6H_6 gale %6.4f', q(i));
st{3} = sprintf('La quantit de chlorine par mole de C_6H_6 pour l''obtenir,')
st{4} = sprintf('est de %6.4f', p(i));
text(0, 0.55, st, 'fontsize', 14)
