vertici=Mfb.cubo(4)

w1=-2+1i
w2=-1+1i
w3=1+1i
w4=3+1i

f=G.lin(1.2*c(w1,w2,w3,w4)-0.25)

Gs.postscript("21-cubotesto.ps",5,2.78)
latox=c(-4,5); latoy=c(-0.2,4.8)
Gd.grafica(latox,latoy)
t=seq(0,2*pi,by=0.01)
cerchio=Gf.cerchio(t,0.1)
for (p in vertici)
polygon(cerchio+f(p),col="yellow")
r=length(vertici)
for (i in 1:(r-1)) {u=vertici[[i]]
for (j in (i+1):r) {v=vertici[[j]]
if (Mfb.hamm(u,v)==1)
Gdf(f(u),f(v),punta=0,d1=0.1,d2=0.1)}}
for (p in vertici) {z=f(p)
text(Re(z)-0.2,Im(z),
  paste(p,collapse=''),
  cex=0.4,pos=4)}
dev.off()