function gaussian(sigma_x, sigma_y)
x = -4:0.1:4;
y = -4:0.1:4;
[x, y] = meshgrid(x, y);
z = exp(-(x.^2 / 2 * sigma_x^2 + y.^2 / 2 * sigma_y^2)) / (2 * pi * sigma_x * sigma_y);
mesh(x, y, z)
end



function gaussian(sigma_x, sigma_y)
x = -4:0.1:4;
y = -4:0.1:4;
[x, y] = meshgrid(x, y);
z = exp(-(x.^2 / 2 * sigma_x^2 + y.^2 / 2 * sigma_y^2)) / (2 * pi * sigma_x * sigma_y);
mesh(x, y, z)
end