diff options
| author | Eugeniy Mikhailov <evgmik@gmail.com> | 2009-12-09 23:12:07 +0000 |
|---|---|---|
| committer | Eugeniy Mikhailov <evgmik@gmail.com> | 2009-12-09 23:12:07 +0000 |
| commit | 3c4e365adaaa98a6c967bc679ca14008c94d2581 (patch) | |
| tree | 8dd8de9ce82003ad1f51279e1015e08b004a2cb2 /useful_functions.m | |
| parent | 596fbba9a1e802f34e497f893ff6c2194cc688c3 (diff) | |
| download | multi_mode_eit-3c4e365adaaa98a6c967bc679ca14008c94d2581.tar.gz multi_mode_eit-3c4e365adaaa98a6c967bc679ca14008c94d2581.zip | |
added calculation of the adjasent sideband sucsectibilty in for wave mixing
Diffstat (limited to 'useful_functions.m')
| -rw-r--r-- | useful_functions.m | 11 |
1 files changed, 11 insertions, 0 deletions
diff --git a/useful_functions.m b/useful_functions.m index 5522afb..2ae4ced 100644 --- a/useful_functions.m +++ b/useful_functions.m @@ -141,3 +141,14 @@ function [rhoLiouville_dot, L]=constran_rho_and_match_L( rhoLiouville_dot(1)=1; endfunction +% calculate sucseptibility for the field at given frequency index +function kappa=sucseptibility(wi, rhoLiouville, dipole_elements, Nlevels, Nfreq) + + rho=rhoOfFreq(rhoLiouville, wi, Nlevels, Nfreq); + kappa=0; + for i=1:Nlevels + for j=1:Nlevels + kappa+=dipole_elements(j,i)*rho(i,j); + endfor + endfor +endfunction |