clear all;
close all;

v = [];
velocity0;
vr = v;

v = [];
simulation0;
vs = v;

clear v;

plot(vs(:,1),vs(:,4),vr(:,1),vr(:,4));

figure(2)