% Decomposition of ozone. 
% Atkins and Jones, first edition p. 626


[t,u]=ode_cg1('ozone', [0, T], [0; 1; 0.5],.00001); 
figure(1)
clf
plot(t,u)

x=ozonea(u);

hold on 

plot(t,x,'k.')

axis([0 t(length(t)) 0 3])

