See ths example
clc,clear
[t,r] = meshgrid(0:10:360,0:0.2:5);
[x,y] = pol2cart(t*pi/180,r);
z = sin(r)./r;
z0 = 1; % z start coordinate
r0 = 2; % radius of a tool
for i = 1:9
x0 = -6; % x start coordinate
z0 = z0 - 0.1; % z increment
for j = 1:7
ind = ((x-x0).^2+y.^2 < r0^2) & (z>z0);
z(ind) = z0;
surf(x,y,z)
x0 = x0 + 1;
pause(0.5)
end
end
