Add a separable Laplacian of Gaussian filter

My favorite filter.  This implementation takes about 0.030 seconds to
filter a 512 x 512 Uint16 image with a 21 x 21 kernel.  The standard
filter, `imfilter(img, imlog(2.5))`, takes 8 seconds for the same image
and kernel sizes.