t=0:0.1:2*pi+0.1;
x=cos(t);
y=sin(t);
kk=length(t);
pp=0:0.005:1;
N=length(pp)
figure
hold on;
for k=1:1:N
       plot3(x,y,pp(k).*ones(1,kk))
end
grid
colormap(cool)
view([-30 -6 50])

cylinder(1,30)