def cr(z):
    return (z.real(),z.imag());
verts = [cr(exp(2*pi*I*4*n/13)) for n in [0..12]];
P = polygon(verts,axes=False,color='gold');
Q = polygon(verts,axes=False,fill=False,color='tomato');
(P+Q).show()