[x y] = meshgrid(-3:0.1:3, -3:0.1:3);
z = 1 ./ (log(1 + 2 * x.^2 + y.^2) + 1);
surf(x, y, z);