Matlab and Octave

z = 4*sin(2*y+sin(2*y+x))+cos(5/2*(x)) ;
[x,y] = meshgrid(-1:0.1:1);
z = 4*sin(2*y+sin(2*y+x))+cos(5/2*(x));
surf(x,y,z)