function imshow(im)
image(uint8(cat(3,im,im,im)));
axis equal
axis off