plot[log(1-2x-2y+2x^2+2y^2)+2arctan((x-y)/(1-x-y))=0,x>0,y>0]