Track

R=2;C=2;

vz = 10; % velocity constant
a = -32; % acceleration constant

subplot(R,C,1);

t = 0:.1:1;
z = vz*t;
% z = vz*t + 1/2*a*t.^2;

vx = 2;
x = vx*t;
%x = vx*t + 1/2*a*t.^2;

vy = 3;
y = vy*t;
% y = vy*t + 1/2*a*t.^2;

u = gradient(x);
v = gradient(y);
w = gradient(z);
scale = 0;

quiver3(x,y,z,u,v,w,scale);
view([70,18]);

subplot(R,C,2)
% z = vz*t;
z = vz*t + 1/2*a*t.^2;

vx = 2;
x = vx*t;
% x = vx*t + 1/2*a*t.^2;

vy = 3;
y = vy*t;
% y = vy*t + 1/2*a*t.^2;

u = gradient(x);
v = gradient(y);
w = gradient(z);
scale = 0;

quiver3(x,y,z,u,v,w,scale);
view([70,18]);

subplot(R,C,3)
% z = vz*t;
z = vz*t + 1/2*a*t.^2;

vx = 2;
% x = vx*t;
x = vx*t + 1/2*a*t.^2;

vy = 3;
y = vy*t;
% y = vy*t + 1/2*a*t.^2;

u = gradient(x);
v = gradient(y);
w = gradient(z);
scale = 0;

quiver3(x,y,z,u,v,w,scale);
view([70,18]);

subplot(R,C,4)
% z = vz*t;
z = vz*t + 1/2*a*t.^2;

vx = 2;
% x = vx*t;
x = vx*t + 1/2*a*t.^2;

vy = 3;
%y = vy*t;
y = vy*t + 1/2*a*t.^2;

u = gradient(x);
v = gradient(y);
w = gradient(z);
scale = 0;

quiver3(x,y,z,u,v,w,scale);
view([70,18]);