added easy field decomposition in left,right and linear polarizations

1 files changed, 17 insertions, 4 deletions

diff --git a/liouville.m b/liouville.m
index 7d1bd08..ab583b8 100644
--- a/liouville.m
+++ b/liouville.m
@@ -89,10 +89,23 @@ fflush (stderr);
 [xi_linear, xi_left, xi_right]=parcellfun(2, @susceptibility_steady_state_at_freq, problems_cell_array);
 
 
-figure(1); plot(detuning_freq, real(xi_left)); title("probe dispersion");
-figure(2); plot(detuning_freq, imag(xi_left)); title("probe absorption");
-%figure(3); plot(detuning_freq, real(kappa_m)); title("off resonant sideband dispersion");
-%figure(4); plot(detuning_freq, imag(kappa_m)); title("off resonant absorption");
+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;
+
 
 elapsed_time = etime (clock (), t0)