figure
t=0:(1/24)*pi:2*pi;nx=1;ny=0;nz=0;
ux=0;uy=1;uz=0;
u1=ux/sqrt(ux*ux+uy*uy+uz*uz);
u2=uy/sqrt(ux*ux+uy*uy+uz*uz);
u3=uz/sqrt(ux*ux+uy*uy+uz*uz);
vx=0;vy=0;vz=1;
v1=vx/sqrt(vx*vx+vy*vy+vz*vz);
v2=vy/sqrt(vx*vx+vy*vy+vz*vz);
v3=vz/sqrt(vx*vx+vy*vy+vz*vz);
for k=1:20
  x0=50*rand(1);y0=50*rand(1);z0=50*rand(1);
    r=10*rand(1);
x=x0+r*(u1*cos(t)+v1*sin(t));
y=y0+r*(u2*cos(t)+v2*sin(t));
z=z0+r*(u3*cos(t)+v3*sin(t));
plot3(x,y,z)
fill3(x,y,z,'r');
hold on
end
grid on;
xlabel
ylabel
zlabel