Построить Тор с анимацией шариков на его поверхности

clear;
m=100;
n=40;
a=0;
b=2*pi;
c=0; d=2*pi;
R1=2;
R2=2;
R3=2;
R=0.5+max(max(R1,R2),R3);
for i=1:n
    u=a+(b-a)*(i-1)/(n-1);
    for j=1:n
        v=a+(d-c)*(j-1)/(n-1);
        x1(i,j)=R1*sin(u)*cos(v);
        y1(i,j)=R2*sin(u)*sin(v);
        z1(i,j)=x1(i,j)^2-y1(i,j)^2;
        
    end
end
mesh(x1,y1,z1);
hold on;
for i=1:m
    fi=2*pi*rand(1);
    psi=pi*rand(1);
    x(i)=cos(fi)*sin(psi);
    y(i)=sin(fi)*sin(psi);
    z(i)=x(i)^2-y(i)^2;
    ball(i)=plot3(x,y,z,'ro');
end
pause
 
for t=1:200
    for i=1:m
        fi=2*pi*rand(1);
        psi=pi*rand(1);
    x(i)=2*cos(fi)*sin(psi);
    y(i)=2*sin(fi)*sin(psi);
    z(i)=x(i)^2-y(i)^2;
    set(ball(i),'XData',x,'YData',y,'ZData',z);
end
pause(0.1);
end

clear;
R=5;
r=3;
n=50;
m=20;
q=8;
h1=2;
h2=5;
for i=1:h1
    for j=1:h2
fi = 2*pi*rand(1);
psi = 2*pi*rand(1);
X1(i,j) = (R + r.*cos(fi)).*cos(psi);
Y1(i,j) = (R + r.*cos(fi)).*sin(psi);
Z1(i,j) = r.*sin(fi);
end
end
balls=plot3(X1,Y1,Z1,'blacko'); 
 %-------
pause;
for t=1:20 %vremya krycheniya
  for i=1:h1
    for j=1:h2
fi = 2*pi*rand(1);
psi = 2*pi*rand(1);
X1(i,j) = (R + r.*cos(fi)).*cos(psi);
Y1(i,j) = (R + r.*cos(fi)).*sin(psi);
Z1(i,j) = r.*sin(fi);
end
end
%TOR
for i=1:q
    for j=1:q
fi = 2*pi*rand(1);
psi = 2*pi*rand(1);
X2(i,j) = (R + r.*cos(psi)).*cos(fi);
Y2(i,j) = (R + r.*cos(psi)).*sin(fi);
Z2(i,j) = r.*sin(psi);
end
end
%VSE
mesh(X2,Y2,Z2);
grid on; hold on;
plot3(X1,Y1,Z1,'blacko');  hold off;
pause(0.1)
end