
% tests the g(sigma_delta) function

function img = footest(image)

image = imresize(image,[size(image,1)*2 size(image,2)*2]);
image = blurGauss(image,1.0);

k = 2.^(18/3);
sigma = 1.6;
sigmadelta = sigma*sqrt(k.^2-1);

eric1 = blurGauss(image,sigma);
eric2 = blurGauss(eric1,sigmadelta);

eric3 = blurGauss(image,sigma*k);

imshow(abs(eric2-eric3))
img = abs(eric2-eric3);