clf
disp('Overlay the slope/vector field and the solution curves')
pause
help quiver
pause
disp('First define the grid points where to plot arrows')
pause
[tt,yy]=meshgrid(0:5:30,0:1:12)
pause
disp('Next define the slopes representing the diff equn')
disp('for slope fields always the t-components are always')
disp('equal to one; the y-components equal (dy/dt).')
pause
vt=ones(size(tt))
vy=-(yy-7)/10;
pause
disp('With basepoints and arrows defined quiver is ready:')
pause
quiver(tt,yy,vt,vy)
pause
disp('To OVERLAY the solution curves, i.e. keep the field')
pause
hold on
pause
disp('recalculate the solution curves')
pause
[t,y]=ode23('RC',[0 30],[0 2 4 6 8 10])
pause
disp('and overlay their graphs')
plot(t,y)
pause