diff --git a/analyze.m b/analyze.m
new file mode 100644
index 0000000..cb3c6c0
--- /dev/null
+++ b/analyze.m
@@ -0,0 +1,17 @@
+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)
+