[x,y]=meshgrid(-7.5:7.5,-15:15);z=1./(1+sqrt(1+x.^2+(y-7.5).^2))+1./(1+sqrt(1+x.^2+(y+7.5).^2));mesh(x,y,z);
