1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
|
#!/usr/bin/octave -qf
% reads parameters file (argv1) and output results of the calculations into output file (argv(2)
arg_list = argv ();
if (nargin < 2)
error("need at least two parameters");
end
params_file=arg_list{1};
output_fname=arg_list{2};
load(params_file);
[xi_linear, xi_left, xi_right, E_field_pos_freq, light_positive_freq]=susceptibility_problem(detuning_freq, Ep, psi_el, B_field, theta, phi);
save(output_fname, 'detuning_freq', 'Ep', 'psi_el', 'B_field', 'theta', 'phi' ...
, 'xi_linear', 'xi_left', 'xi_right' ...
, 'E_field_pos_freq', 'light_positive_freq' ...
);
|
