diff options
Diffstat (limited to 'useful_functions.m')
| -rw-r--r-- | useful_functions.m | 9 |
1 files changed, 6 insertions, 3 deletions
diff --git a/useful_functions.m b/useful_functions.m index e0c9e8d..bc8f709 100644 --- a/useful_functions.m +++ b/useful_functions.m @@ -236,13 +236,16 @@ endfunction function xi=susceptibility_steady_state_at_freq( atom_field_problem) % find steady state susceptibility at particular modulation frequency element % at given E_field - L0m = atom_field_problem.L0m ; - polarizability_m = atom_field_problem.polarizability_m ; - dipole_elements = atom_field_problem.dipole_elements ; + %L0m = atom_field_problem.L0m ; + %polarizability_m = atom_field_problem.polarizability_m ; + %dipole_elements = atom_field_problem.dipole_elements ; E_field = atom_field_problem.E_field ; modulation_freq = atom_field_problem.modulation_freq ; freq_index = atom_field_problem.freq_index ; + % load from file values of 'L0m', 'polarizability_m', 'dipole_elements'; + load( atom_field_problem.atom_properties_fname ); + rhoLiouville=rhoLiouville_steady_state(L0m, polarizability_m, E_field, modulation_freq); xi=susceptibility(freq_index, rhoLiouville, dipole_elements); endfunction |