1;

load 'xi_vs_detuning.mat' ;

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;