plot z(x,y) =(x^2+y^2+x*y+x+1)