contourplot[z=√(x*(y+1/y)),x/y>=0]