#!/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' ...
	);