clear all; [x,y]=meshgrid([-10:0.1:10],[-10:0.1:10]); z = exp(1-0.05*(x-1).^2-0.2*(y+2).^2)+sin(x).^2.*sin(y).^2.*exp(0.05*(1-x.^2-y.^2)).*y.^2; h=surfc(x,y,z,'EdgeColor','none'); zlim([-8,8]); for i = 2:length(h) newz = get(h(i),'zdata') - 8; set(h(i),'zdata',newz); end;