MATLAB——小球碰撞

clear all
syms aa bb
R=5;i=7;                 %设置小球半径以及小球数量
X=randi([10 90],1,i);       %在[10，90]间取整数作为每个小球圆心的X坐标
Y=randi([10 90],1,i);       %在[10，90]间取整数作为每个小球圆心的Y坐标
t=0:pi/100:2*pi;
x=2.*rand(1,i)+1;          %在[1，3]间取随机数作为每个小球x轴方向速度
y=2.*rand(1,i)+1;          %在[1，3]间取随机数作为每个小球y轴方向速度
for n=1:10000
    for v=1:i
        if 100-X(v)<R          %小球与右器壁相碰的速度判定
            x(v)=-abs(x(v));
        end
        if X(v)<R              %小球与左器壁相碰的速度判定
            x(v)=abs(x(v));
        end
        if Y(v)<R              %小球与下器壁相碰的速度判定
            y(v)=abs(y(v));
        end
        if 100-Y(v)<R          %小球与上器壁相碰的速度判定
            y(v)=-abs(y(v));
        end
        for b=v:i
            if sqrt((X(v)-X(b))^2+(Y(v)-Y(b))^2)<=2*R   %小球间相碰速