program dobbeltall
      double precision a,b
c Double representation
2     continue      
      write(*,*) 'Double ='
      read (*,*) a
      call dint(a,b)
      write(*,*) b
      goto 2
      end
c
      subroutine dint(a,b)
      dimension a(2),b(2)
      call rint(a,b)
      return
      end
c      
      subroutine rint(ia,ib)
      dimension in(2,32),ia(2),ib(2)
      do j=1,2
      do i=1,32
      in(j,i)=iand(1,ishft(ia(j),1-i))
      enddo
      enddo
      write(*,1) ((in(j,i),i=32,1,-1),j=1,2)
1     format(1x,64(i1))
      ib(1)=ia(1)
      ib(2)=ia(2)
      return
      end