faraday/tasks_pp.m


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
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));