disp('Higher order DEs -- recall the bungee -- are first')
disp('rewritten as a system of first order DEs -- take a')
disp('look as this example of a damped forced oscillator.')
disp('y is a column vector containing position and speed')
pause
[t,y]=ode23('forced',[0,20],[2;1])
pause
plot(t,y)
pause
plot(t,y(:,1),'ro',t,y(:,2),'b:')