function prob_2_15
% 
% Written by:
% -- 
% John L. Weatherwax                2005-05-07
% 
% email: wax@alum.mit.edu
% 
% Please send comments and especially bug reports to the
% above email address.
% 
%-----

epsilon = 1e-4;
options = odeset('Stats','on');

tic;
%[t,y]=ode15s(@ode, [ 0, 1/epsilon], epsilon, options );
ode15s(@ode, [ 0, 1/epsilon], epsilon, options );
toc

tic;
%[t,y]=ode15s(@ode, [ 0, 2/epsilon], epsilon, options );
ode15s(@ode, [ 0, 2/epsilon], epsilon, options );
toc

tic;
[t,y]=ode45(@ode, [ 0, 1/epsilon], epsilon, options );
toc

tic;
[t,y]=ode45(@ode, [ 0, 2/epsilon], epsilon, options );
toc

%===================================

function dydt = ode(t,y)
dydt = y^2*(1-y);