data_dir='results/';
output_dir='results/';
detuning_freq=0;

gmg=.7; % gyro magnetic ration for ground level
zeeman_splitting=+0.1;
Nsteps=100;
B_fields=linspace(-zeeman_splitting/gmg, zeeman_splitting/gmg, Nsteps);

%[psr_rad]=psr_vs_detuning(Ep, psi_el, B_field, theta, phi) 

% phi is angle between linear polarization and axis x
%phi=pi/4;
phi=0;
% theta is angle between lab z axis (light propagation direction) and magnetic field axis (z')
theta=0;
% psi_el is the ellipticity parameter (phase difference between left and right polarization)
psi_el=-30/180*pi*e-8;

Ep=sqrt(0.1);

%[psr_rad_smEp_pos_el]=psr_vs_detuning(detuning_freq, Ep, psi_el, B_field, theta, phi) ;
[psr_rad_smEp_pos_el]=faraday_vs_B(detuning_freq, Ep, psi_el, B_fields, theta, phi) ;

plot(B_fields, psr_rad_smEp_pos_el)