% Kinky plots
[t, x] = ode45('test', [0 10*pi], 3); 
plot(t, x, '+-')
fprintf('Press <Enter> to continue\n')
pause
[t, x] = ode45('test', [0 10*pi], 3, odeset('refine',10));
l = (length(t) - 1) / 10
tp = t(10*(0:l) + 1);
xp = x(10*(0:l) + 1);
plot(t,x, tp,xp,'o', tp,xp,'--')
