% Stiffdemo demonstrates the effect of stiffness on an
%           explicit Runge-Kutta ODE solver.

close all

tspan=[0 100]; 
y0 = 0.1; 
options = odeset('Stats','on');
[t,y]=ode23('logde',tspan,y0,options);

plot(t,y),grid,axis([0 100 0.98 1.02])
xlabel('T'), ylabel('Y')
title('Solution of logistic ODE by ode23')

h = diff(t); figure, plot(h'),grid
xlabel('Step Number'), ylabel('Step Length')
title('Stepsize history in solution of logistic ODE by ode23')