function range = imageshow(im, range)% Show image with proper scalling% First exclude the Inf value:realim = im(find(isfinite(im)));% By default range is within 2 standard deviation from the meanif ~exist('range', 'var')    range = mean(realim(:)) + std(realim(:)) * [-2, 2];endimagesc(im, range);axis image;