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

[ts2, vals2] = ode45(g, [0 10], 1);
subplot(2,1,1)
plot(ts2, vals2)
xlabel('t value')
ylabel('solution')
title('Results using ode45')
subplot(2,1,2)
plot(ts2(1:(length(ts2)-1)), diff(ts2),'k*', 'markersize', 2)
xlabel('t value')
ylabel('stepsize')
pause