subroutine eint(a,ee,ek)
c     evaluate elliptic integrals
c-----------------------------------------------------------------------
c
      include 'ucom.h'
c
c-----------------------------------------------------------------------
c*
      b=1.0-a
      c=log(b)
      ee=1.0+b*(.4630106-0.2452740*c+b*(0.1077857-0.04125321*c))
      ek=1.38629436-.5*c+b*(0.1119697-0.1213486*c+b*(.07253230-.02887472
     1*c))
      return
      end
c++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*