%     A.-H. DISPERSION CURVES
%     THIS CALCULATES THE SQUARES OF THE REFRACTIVE
%     INDEX WITH THE PLUS/MINUS SIGNS.  NO ATTEMPT IS
%     MADE TO DETERMINE ORDINARY, EXTRAORDINARY RAYS
%     OUTPUT in FILE ah.out
 
% Magnetic field parameters     
YL=0.1;
YT=0.2;
Z=0.001;

% Start, Stop, and Step in X scan
XMIN=0;
XMAX=2;
DELX=0.01;

i1=1;
for x=XMIN:DELX:XMAX
  X=x;
  if (X==1.0) X=0.9999; end;
  BR=0.5*YT^2/(1-X-j*Z); 
  DENP=1-j*Z-BR+sqrt(BR^2+YL^2); 
  DENM=1-j*Z-BR-sqrt(BR^2+YL^2); 
  N2P=1-X/DENP;
  N2M=1-X/DENM; 
  XX(i1)=X;
  N2PP(i1)=N2P;
  N2MM(i1)=N2M;
  i1=i1+1;
end
  
  Q=[XX; real(N2PP); imag(N2PP); real(N2MM); imag(N2MM)]';

save ah.out Q -ascii