%% Create surface data [X,Y] = meshgrid(-8:.5:8); R = sqrt(X.^2 + Y.^2) + eps; Z = sin(R)./R; %% Plotting surface surf(X,Y,Z,'EdgeColor','black') %% Adding your color palette and colorbar colormap(colors) colorbar()