%% Falling object model

%f = @(t,y) 9.8 - y/5;
%f = @(t,y) -(y-1)*t;
f = @(t,y) y-t;

X =-1:.5:15;
%X=-1:1:15;
%Y = 20*(0:.1:1.5);
Y = 2*(-1:.1:1);


% Plot direction field
dirfield(f,X,Y)
xlabel('t');
ylabel('y');
title('dy/dt=y-t');