function plota_celulas(K,tam,CD,LD,NC,cor)

LD_c = LD ./ NC; %tamanho de uma celula
MD = CD - 0.5*LD; %canto inferior do dominio
nk = size(K,1);

for k = 1:nk
    posk = K(k,:);
    MD_c = MD + (posk .* LD_c); %canto inferior da celula 
    MD_r = MD_c + (1 - tam)/2*LD_c; %canto inferior do retangulo
    LD_r = tam*LD_c; %tamanho do retangulo
    hold on
    rectangle('Position',[MD_r(1),MD_r(2),LD_r(1),LD_r(2)],...
        'FaceColor',cor, 'EdgeColor','black')
end



