sc=7;
phi=pi/7;
co = @(z) cos(z*4*pi/7+phi)*sc;
si = @(z) sin(z*4*pi/7+phi)*sc;
[sx,sy,sz] = sphere(6);
col=[
1 0.89 0 % a
0.5 0.5 0.5 % b
0.5 0.5 0.5 % c
0.79 0.37 0.58 % d
0.44 0.6 0 % e
0.5 0.5 0.5 % f
1 0.6 0 % g
1 0.89 0 % a
];
z=0:0.01:7;
x=co(z);
y=si(z);
n=length(z);
figure(1)
clf
for i=1:n-1
plot3(x(i:i+1),y(i:i+1),z(i:i+1),...
'LineWidth',2, 'Color',[0.5,0.5,0.6])
hold on
end
for z=0:7
p=surf(sx+co(z),sy+si(z),sz+z);
colormap(col(z+1,:))
alpha(p, 0.4)
end