1 files changed, 19 insertions, 0 deletions
diff --git a/faraday/tasks_pp.m b/faraday/tasks_pp.m
new file mode 100644
index 0000000..bc81b07
--- /dev/null
+++ b/faraday/tasks_pp.m
@@ -0,0 +1,19 @@
+
+fnames=glob('results/*.mat');
+
+Nsteps=length(fnames)
+B_field=zeros(1,Nsteps);
+xi_linear=zeros(1,Nsteps);
+xi_left=zeros(1,Nsteps);
+xi_right=zeros(1,Nsteps);
+
+% read the information from resulting files
+for i=1:Nsteps
+	d=load(fnames{i});
+	B_field(i)=d.B_field;
+	xi_linear(i)=d.xi_linear;
+	xi_left(i)=d.xi_left;
+	xi_right(i)=d.xi_right;
+endfor
+
+plot(B_field, real(xi_right+xi_left));