j=3
      i=7
      k=iand(j,i)
      print *,k
      k=ibits(i,0,2)
      print *,k
      k=ishft(j,4)
      print *,k
      stop
      end