function z = v(x,y)

       z = (x-1).^2 + y.^2 + 2;