function z = u(x,y)
       
      z = -x.^3  + 3*x.*y.^2 +2;