equivalence (x,ix),(y,iy),(z,iz),(v,iv)
      read *,x
      print *,'x =',x,' hex value ='
      call emas3phex(ix)
      v=0.000001*x
      write(*,1)
      print *,'0.000001*x =',v,' hex value ='
      call emas3phex(iv)
      y=1.0
      do 10 j=1,10
      y=(x/y+y)/2
      write(*,1)
      print *,'y =',y,'hex value ='
      call emas3phex(iy)
      z=abs(y*y-x)
      write(*,1)
      print *,'abs(y*y-x) =',z,'hex value ='
      call emas3phex(iz)
      if(z.lt.0.000001*x) goto 20
  10  continue
  20   stop
  1   format(1h ,)
      end