動かし方はそのまま踏襲して、下記の様に変更すれば長方形を動かす事が出来ます。
rectplotは自作の関数です。一番下に書きました。
h1 = rectplot(y1(1) , 7.5 , ' ro ');
h2 = rectplot(7.5 , y2(1) , ' bo ');
h3 = rectplot(y3(1) , -7.5 , ' ro ');
h4 = rectplot(-7.5 , y4(1) , 'bo');
h5 = rectplot(y5(1) , 7.5 , 'ro');
h6 = rectplot(7.5 , y6(1) , 'bo');
h7 = rectplot(y7(1) , -7.5 , 'ro');
h8 = rectplot(-7.5 , y8(1) , 'bo');
set(h1,'Position',[y1(n),7.5 ,10, 10])
set(h2,'Position',[7.5 ,y2(n),10, 10]);
set(h3,'Position',[y3(n),-7.5 ,10, 10]);
set(h4,'Position',[-7.5 ,y4(n),10, 10]);
set(h5,'Position',[y5(n),7.5 ,10, 10]);
set(h6,'Position',[7.5 ,y6(n),10, 10]);
set(h7,'Position',[y7(n),-7.5 ,10, 10]);
set(h8,'Position',[-7.5 ,y8(n),10, 10]);
別途下記関数を追加してください。
function h = rectplot(x,y,style)
cl = style(1); % 'ro'や'bo'の最初の一文字だけ色データとして貰う
h = rectangle('Position',[x,y,10,10],'EdgeColor',cl);
end