[x,y]= meshgrid(-4:0.02:4, -4:0.02:4);
z = ( y.^4-18* y.^2 + 6*x.^2+ 2*x.^3  )  ;
grid on
mesh(x,y,z)

