equivalence (x,ix),(y,iy)
      x=78.4
      call emas3phex(ix)
      y=1.0
      do 10 j=1,10
      y=(x/y+y)/2
      print *,y
      if(abs(y*y-x).lt.0.000001*x) goto 20
  10  continue
  20  call emas3phex(iy)
      stop
      end