[t,y]=meshgrid(0.1:.1:2); dt = ones(20); dy = (-y)./(t+1); quiver(t,y,dt,dy,1); hold on; for i=0:.4:2 initial = i; time = [0.1,2]; [T,Y] = ode45(@ode,time,initial); plot(T,Y,'LineWidth',2,'Color',[1 0 0]); end