g = inline('2*x-x.^2', 't', 'x');
[ts, vals] = ode23(g, [0 10], 1);
subplot(2,1,1)
plot(ts, vals)
xlabel('t value')
ylabel('solution')
title('Results using ode23')
subplot(2,1,2)
plot(ts(1:(length(ts)-1)), diff(ts),'k.')
xlabel('t value')
ylabel('stepsize')
wait

[ts, vals] = ode45(g, [0 10], 1);
subplot(2,1,1)
plot(ts, vals)
xlabel('t value')
ylabel('solution')
title('Results using ode45')
subplot(2,1,2)
plot(ts(1:(length(ts)-1)), diff(ts),'k.')
xlabel('t value')
ylabel('stepsize')
wait

[ts, vals] = ode113(g, [0 10], 1);
subplot(2,1,1)
plot(ts, vals)
xlabel('t value')
ylabel('solution')
title('Results using ode113')
subplot(2,1,2)
plot(ts(1:(length(ts)-1)), diff(ts),'k.')
xlabel('t value')
ylabel('stepsize')