plot[ [y+x,y-x,1] . [[2,1,1],[1,2,1],[1,1,2]] . [y+x,y-x,1] /sqrt[2] ]