x=0.2:0.001:0.5;
y=1./x;
[X,Y,Z]=cylinder(y,30);
mesh(X,Y,Z)