1;


load '/tmp/xi_vs_detuning.mat' ;

figure(1);
hold off;
plot(detuning_freq, real(xi_left-xi_right), '-');  
title("differential real xi");
xlabel("two photon detuning");

figure(2);
hold off;
plot(detuning_freq, imag(xi_left-xi_right), '-');  
title("differential imag xi");
xlabel("two photon detuning");

figure(3);
hold off;
plot(detuning_freq, imag(xi_left), '-', detuning_freq, imag(xi_right), '-');  
title("imag xi");
xlabel("two photon detuning");

figure(4);
hold off;
plot(detuning_freq, real(xi_left), '-', detuning_freq, real(xi_right), '-');  
title("real xi");
xlabel("two photon detuning");

%figure(1); 
	%hold off;
	%plot(detuning_freq, imag(xi_linear), '-1;linear;');  
	%hold on;
	%plot(detuning_freq, imag(xi_left),   '-2;left;');  
	%plot(detuning_freq, imag(xi_right),  '-3;right;');  
	%title("probe absorption");
	%hold off;
%figure(2); 
	%hold off;
	%plot(detuning_freq, real(xi_linear), '-1;linear;');  
	%hold on;
	%plot(detuning_freq, real(xi_left),   '-2;left;');  
	%plot(detuning_freq, real(xi_right),  '-3;right;');  
	%title("probe dispersion");
	%hold off;