Как сделать окружности разный цветов

clear all    %   очистка переменных.
r=1;  % неизменная переменная
 
for x0=-1:1/12:1  % цикл для центра окружности
 i=1;                        % счетчик второго цикла
 for phy=0:pi/24:2*pi;       % цикл для переменной phi
   X1(i)=x0+r*cos(phy);
   Y1(i)=r*sin(phy);
   i=i+1;
 end
 grid on
 hold on
 plot(X1,Y1);
end

clear all % очистка переменных.
r=1; % неизменная переменная
color = 0;
for x0=-1:1/12:1 % цикл для центра окружности
    i=1; % счетчик второго цикла
    for phy=0:pi/24:2*pi; % цикл для переменной phi
        X1(i)=x0+r*cos(phy);
        Y1(i)=r*sin(phy);
        i=i+1;
    end
    grid on
    hold on    
    plot(X1,Y1,'Color',[1-1/25*color,0,1-1/25*color]);
    color = color+1;
end

clear all % очистка переменных.
r=1; % неизменная переменная
str = ['y' , 'm', 'c', 'r', 'g', 'b', 'k']; % перечень цветов
color =  randi(8,1,25);
str(randi(8,1,1));
for x0=-1:1/12:1 % цикл для центра окружности
    i=1; % счетчик второго цикла
    for phy=0:pi/24:2*pi; % цикл для переменной phi
        X1(i)=x0+r*cos(phy);
        Y1(i)=r*sin(phy);
        i=i+1;
    end
    grid on
    hold on    
    plot(X1,Y1,'Color',str(randi(7,1,1)));   % случайный цвет из перечня
end

clear all % очистка переменных.
r=1; % неизменная переменная
% str = ['y' , 'm', 'c', 'r', 'g', 'b', 'k'];
str = {[0 0 1], [1 0 1], [0 0 0], [1 0 0], [0 1 0], [0.5 1 0.5]}; % Цвета сами подбирайте. Это обычная RGB-политра
for x0=-1:1/12:1 % цикл для центра окружности
    i=1; % счетчик второго цикла
    for phy=0:pi/24:2*pi; % цикл для переменной phi
        X1(i)=x0+r*cos(phy);
        Y1(i)=r*sin(phy);
        i=i+1;
    end
    grid on
    hold on
%     plot(X1,Y1,'Color',str(randi(size(str,2),1)));
    plot(X1,Y1,'Color',str{randi(size(str,2),1)});   
end

clear all % очистка переменных.
r = 1;
x0=-1:1/12:1; % центры окружности
fi=0:360; % переменная phi
[X0 FI] = meshgrid(x0,fi); % создаем матрицы
X1 = X0+r*cosd(FI);
Y1 = r*sind(FI);
plot(X1,Y1)
grid on
axis equal