basis_transformation;

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

figure(1);
plot(B_field, real((xi_right+xi_left))/sqrt(2));
title('phase of linear x-polarization');

figure(2);
plot(B_field, imag((xi_right+xi_left))/sqrt(2));
title('absorption of linear x-polarization');

figure(3);
plot(B_field, real((1i*xi_right-1i*xi_left))/sqrt(2));
title('phase of linear y-polarization');

figure(4);
plot(B_field, imag((1i*xi_right-1i*xi_left))/sqrt(2));
title('absorption of linear y-polarization');