function G = my_gaussian_filter(x,y,alpha)

G = zeros(x,y);
% find centers
cx = floor(x/2)+1;
cy = floor(y/2)+1;

for i=1:x
    for j=1:y
        r = sqrt((i-cx)^2 + (j-cy)^2);
        G(i,j) = exp(-alpha*r^2) ;
    end
end

end


