
% draws circle and shows random uniform distribution 
thickness=3;
for i=1:5000

%getting random points
    theta1=2*180*rand();            % problem theta is never going to be zero or 360
    %z1=(randi(10)+rand())*(-1)^(randi(2));
    z1=(randi(650)+rand())/100; % in cm units. height of the sample 
    R1=4.25+thickness*(rand())^.5;  %creates random uniform distribution
    %between the range of the thickness of the sample
    x1=R1*cosd(theta1);
    y1=R1*sind(theta1);
    plot(x1,y1)
    plot3(x1,y1,z1)
    hold on  
end 
grid on
%rotate3d on