z=peaks(50);
surfl(z)
shading interp
colormap copper