ContourPlot3D[(x+y)(x+z)(yz)=xy-xz-y^2]