summaryrefslogtreecommitdiff
path: root/liouville.m
diff options
context:
space:
mode:
authorEugeniy Mikhailov <evgmik@gmail.com>2009-12-21 22:48:07 +0000
committerEugeniy Mikhailov <evgmik@gmail.com>2009-12-21 22:48:07 +0000
commitc680b440e6963ab1f507ca0fcefdc3ce982db8af (patch)
tree935b6b2c62f5eaefa4fc407e4e343a05e3319916 /liouville.m
parent9e39b7b6b7a74414551a70ea8ccdb11a98023ef5 (diff)
downloadmulti_mode_eit-c680b440e6963ab1f507ca0fcefdc3ce982db8af.tar.gz
multi_mode_eit-c680b440e6963ab1f507ca0fcefdc3ce982db8af.zip
calculation of total absorption is added
Diffstat (limited to 'liouville.m')
-rw-r--r--liouville.m6
1 files changed, 4 insertions, 2 deletions
diff --git a/liouville.m b/liouville.m
index 29b07f9..fd991a5 100644
--- a/liouville.m
+++ b/liouville.m
@@ -87,9 +87,11 @@ fflush (stderr);
% once we define all problems the main job is done here
%kappa_p=cellfun( @susceptibility_steady_state_at_freq, problems_cell_array);
%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);
+%[xi_linear, xi_left, xi_right]=parcellfun(2, @susceptibility_steady_state_at_freq, problems_cell_array);
+total_absorption_vs_detuning=parcellfun(2, @total_field_absorption, problems_cell_array);
-save 'xi_vs_detuning.mat' detuning_freq xi_linear xi_left xi_right ;
+%save 'xi_vs_detuning.mat' detuning_freq xi_linear xi_left xi_right ;
+save 'total_absorption_vs_detuning.mat' detuning_freq total_absorption_vs_detuning;
output_results;