q=(0:4)*2*%pi/5;
x5=sin(q);y5=cos(q);
q=(0:5)*2*%pi/6;
x6=cos(q);y6=sin(q);
ks5=[1:5,1];ks6=[1:6,1];
scf();
a=sin(2*%pi/10)*2;
x6=x6*a;y6=y6*a;
b=cos(2*%pi/10);
b=(a*cos(2*%pi/12)+b)/b;
d5x=(x5+[x5(2:5),x5(1)])/2*b;
d5y=(y5+[y5(2:5),y5(1)])/2*b;
q=2*%pi/5;
for k=1:5;
kk=modulo(8-k,5);...// mod on Matlab
x=x6*cos(kk*q)-y6*sin(kk*q)+d5x(k);...
y=x6*sin(kk*q)+y6*cos(kk*q)+d5y(k);...
plot2d(x(ks6),y(ks6));...
end;
for n=1:5;...
x=x5+d5x(n)*2;y=-y5+d5y(n)*2;...
plot2d(x(ks5),y(ks5));...
end;
for z=1:5;...
px=d5x(z)*2;py=d5y(z)*2;...
for k=1:5;...
kk=modulo(8-k,5);...
x=x6*cos(kk*q)-y6*sin(kk*q)+px-d5x(k);...
y=x6*sin(kk*q)+y6*cos(kk*q)+py-d5y(k);...
if modulo(-k+z+5,5)>1;plot2d(x(ks6),y(ks6));end;
end; end;
qx=6*d5x(2);qy=2*d5y(2);
d5x=-d5x;d5y=-d5y;
x5=-x5;y5=-y5;
for k=1:5;...
kk=modulo(8-k,5);...
x=x6*cos(kk*q)-y6*sin(kk*q)+d5x(k)+qx;...
y=x6*sin(kk*q)+y6*cos(kk*q)+d5y(k)+qy;...
plot2d(x(ks6),y(ks6));...
end;
for k=1:5;...
x=x5+d5x(k)*2+qx;...
y=-y5+d5y(k)*2+qy;...
plot2d(x(ks5),y(ks5));...
end; for z=1:5;...
px=d5x(z)*2;py=d5y(z)*2;...
for k=1:5;...
kk=modulo(8-k,5);...
x=x6*cos(kk*q)-y6*sin(kk*q)+px-d5x(k)+qx;...
y=x6*sin(kk*q)+y6*cos(kk*q)+py-d5y(k)+qy;...
if modulo(-k+z+5,5)>3;plot2d(x(ks6),y(ks6));end;...
end; end;
set(gca(),'isoview','on');
set(gca(),'axes_visible','off');
