plot demo

[Nelson-numerical-software]

N = 100;
angle = linspace(pi/4,9*pi/4,N);
xb = cos(angle) * 0.9;
yb = sin(angle) * 0.9;
xy = cos(angle) * 0.9 + 1;
yy = sin(angle) * 0.9 - 1;
xk = cos(angle) * 0.9 + 2;
yk = sin(angle) * 0.9;
xg = cos(angle) * 0.9 + 3;
yg = sin(angle) * 0.9 - 1;
xr = cos(angle) * 0.9 + 4;
yr = sin(angle) * 0.9;
figure
hold on
plot(xb(1:3*N/4),yb(1:3*N/4),'b','LineWidth',10);
plot(xy(N/4:N),yy(N/4:N),'y','LineWidth',10)
plot(xk(1:3*N/4),yk(1:3*N/4),'k','LineWidth',10);
plot(xy(1:N/4),yy(1:N/4),'y','LineWidth',10);
plot(xb(3*N/4:end),yb(3*N/4:end),'b','LineWidth',10);
plot(xr(1:N/2),yr(1:N/2),'r','LineWidth',10);
plot(xg(1:N),yg(1:N),'g','LineWidth',10);
plot(xk(3*N/4:N),yk(3*N/4:N),'k','LineWidth',10);
plot(xr(N/2:N),yr(N/2:N),'r','LineWidth',10);
axis equal
axis off
ax = gca();
ax.XLim =  [-1.2 5.2];
ax.Color = [1 1 1];

[Nelson-numerical-software]

figure();
theta = pi / 4;
v_o = 1;
g = 9.81;

x_max = v_o^2 * sin(2*theta) / g;
y_max = 0.5 * (v_o * sin(theta))^2 / g;
x = linspace( 0, x_max, 20);
y = tan(theta) * x - 0.5 * g * x .^2 / (v_o*cos(theta))^2;
grid off
for i = 1:length(x)
	plot(x(1:i), y(1:i));
	hold on
	plot(x(i), y(i), 'o');
	axis([0 x_max 0 y_max]);
	sleep(0.1);
end