% Define max and min t-values and spacing
tmin = 0;
tmax = 1.5;
nt = 15;
dt = (tmax-tmin)/nt;

% Define max and min y-values and spacing
ymin=0;
ymax=1.5;
ny = 15;
dy = (ymax-ymin)/ny;

% Set up (t,y) grid points
t=[tmin:dt:tmax];
y=[ymin:dy:ymax];
[T,Y]=meshgrid(t,y);

% Define vector field (1,F) at meshpoints
% Gives direction field for y' = f(t,y) with F(T,Y) = f(t,y) 
One = 0.*T + 1;
F = T - Y.^2;

% Plot direction field
quiver(t,y,One,F)