d=1;
sigma[x_,y_]=-1/(2 Pi (x^2 + y^2 + d^2)^(3/2))
Plot3D[sigma[x,y],{x,-4,4},{y,-4,4},PlotPoints->100, PlotRange->All]