PROGRAM etor
*************************************************************************
* Conversion of dielectric functions to log of reflectivity.            *
* the standard input file must have the format x eps1 eps2              *
* The outputfile comes as: x  r1 r2                                     *
*************************************************************************
      REAL X(100000),e1,e2
      complex epsil(100000),cn,cr
      INTEGER I,mm    
      mm=100000 
      do 10 i=1,mm                                                  
       READ(*,*,END=11) X(i),e1,e2
       epsil(i)=cmplx(e1,e2)
10    continue   
11    mm=i-1
      do 25 i=1,mm
       cn=csqrt((epsil(i)))
       cr=(cn-1.)/(cn+1.)
       write(*,*) x(i),real(cr),aimag(cr)
25    continue
      END