a = imread('PSWebImageBitMap0001.bmp',  'bmp');
b = double(a(:,:,2));

f = inline('fLap(x)');
c = nlfilter(b,[3 3],f);

figure(1)
clf
imagesc(c)
colormap(gray)
title('Laplacian gradient mask')
