From c23b87d2848526cfe478e294509780ffccdda4a0 Mon Sep 17 00:00:00 2001
From: ChrisFeatherston <chris.featherston@carleton.ca>
Date: Sat, 15 Apr 2017 16:16:09 -0400
Subject: [PATCH] Added xl histogram at every simulation to MCtraj and then
 made a gif of it

---
 MCTraj/MCtraj.m | 24 +++++++++++++++++++-----
 1 file changed, 19 insertions(+), 5 deletions(-)

diff --git a/MCTraj/MCtraj.m b/MCTraj/MCtraj.m
index 583212b..f6c5e0f 100644
--- a/MCTraj/MCtraj.m
+++ b/MCTraj/MCtraj.m
@@ -42,7 +42,7 @@
 
 
 xl = [];
-
+figure;
 for n = 1: nSims
     x(1, :) = zeros(1, nTraj);
     y(1, :) = zeros(1, nTraj);
@@ -68,10 +68,24 @@
     end
 
     xl = [xl x(end, :)];
+    histogram(xl,100)
+    title(['Simulation = ' num2str(n)])
+    axis([4e5 14e5 1 35])
+    xlabel('xl')
+    ylabel('count')
+    drawnow
+    frame = getframe(1);
+    im{n} = frame2im(frame);
 end
-
+close;
+%figure;
 if doPlot
-    imwrite(im, map, 'imagefile.gif', 'DelayTime', 0, 'LoopCount', inf);
-    figure
+    for n = 1:nSims
+    [A,map] = rgb2ind(im{n},256);
+        if n ==1
+            imwrite(A, map, 'imagefile.gif', 'gif', 'DelayTime', 0.2, 'LoopCount', inf);
+        else 
+            imwrite(A, map, 'imagefile.gif', 'gif', 'DelayTime', 0.2, 'WriteMode', 'append');
+        end
+    end
 end
-