PROGRAM sig2e
*************************************************************************
* Conversion of dielectric functions to sigma function.                 *
* the standard input file must have the format x sig1 sig2              *
* The outputfile comes as: x eps.                                      *
*************************************************************************
      REAL X(100000),s1,s2
      complex sigma(100000),eps
      INTEGER I,mm    
      mm=100000 
      do 10 i=1,mm                                                  
       READ(*,*,END=11) X(i),s1,s2
       sigma(i)=cmplx(s1,s2)
10    continue   
11    mm=i-1
      do 25 i=1,mm
       eps=1+sigma(i)*(0.,1.)/(x(i)*0.5*0.0333795) 
       write(*,*) x(i),real(eps),aimag(eps)
25    continue
      END