Иллюстрация векторного поля

% Построим зависимость от времени длЯ вектора скорости тела, 
% брошенного под углом к горизонту:
%              Vx = Vx0, Vy = Vy0 - (g/m)*t,
% где Vx0, Vy0 - скорость в начальный момент времени, t - времЯ,
% g - ускорение свободного падениЯ, m - масса тела.
% ТраекториЯ движениЯ тела описываетсЯ законом
%           X = Vx0*t, Y = Vy0*t - (g/2/m)*t^2
% Выберем параметры так, чтобы Vx0 = 0.5 м/с, Vy0 = 0.8 м/с, g/m = 0.8, 
% продолжительность полета составила 2 с, и проследим за изменением 
% скорости тела через каждые 0.2 с.
t = 0:0.2:2;
X = 0.5*t;
Y = 0.8*t.*(1-0.5*t);
Vx(1:length(X)) = 0.5;
Vy = 0.8*(1-t);
% последний аргумент (0.5) - масштабирующий множитель, изменЯющий
% длину стрелок по сравнению с реальной
quiver(X,Y,Vx,Vy,0.5)