From e25827cbfc12e4064c6ad5596369fd7e37e1389a Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Tue, 15 Dec 2009 19:00:39 +0000 Subject: using plot output script now --- liouville.m | 25 +++++-------------------- 1 file changed, 5 insertions(+), 20 deletions(-) (limited to 'liouville.m') diff --git a/liouville.m b/liouville.m index ab583b8..c0c63db 100644 --- a/liouville.m +++ b/liouville.m @@ -64,9 +64,9 @@ for detuning_p_cntr=1:N_detun_steps+1; E_field_drive =[0, 0 , Ed, 0 , Edc, 0, 0 ]; E_field_probe =[0, Ep, 0 , Epc, 0 , 0, 0 ]; E_field_zero =[0, 0 , 0 , 0 , 0 , 0, 0 ]; - E_field.linear = E_field_zero + 0.00*E_field_probe + 0.00*E_field_drive; - E_field.right = E_field_zero + 1.00*E_field_probe + 0.00*E_field_drive; - E_field.left = E_field_zero + 0.00*E_field_probe + 1.00*E_field_drive; + E_field.linear = E_field_zero + 1.00*E_field_probe + 1.00*E_field_drive; + E_field.right = E_field_zero + 1.00*E_field_probe + 1.00*E_field_drive; + E_field.left = E_field_zero + 1.00*E_field_probe + 1.00*E_field_drive; freq_index=freq2index(wp,modulation_freq); atom_field_problem.E_field = E_field; @@ -88,24 +88,9 @@ fflush (stderr); %kappa_p=parcellfun(2, @susceptibility_steady_state_at_freq, problems_cell_array); [xi_linear, xi_left, xi_right]=parcellfun(2, @susceptibility_steady_state_at_freq, problems_cell_array); +save 'xi_vs_detuning.mat' detuning_freq xi_linear xi_left xi_right ; -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; - +output_results; elapsed_time = etime (clock (), t0) -- cgit v1.2.3