1;


load 'total_absorption_vs_detuning.mat' ;

zoom_factor=1e3;
figure(1);
hold off;
plot(detuning_freq, zoom_factor*total_absorption_vs_detuning, '-');  
title("total absorption");
xlabel("two photon detuning");

%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;