# 09-figura-2

n=31
X=0:(n-1)
f = function (x)
{(x^3+23*x^2+7*x+24)%%n}

rota = function (x)
{if (x==Din.omega(x,f)) u=c(0,0.6)
else u=c(0,1)
P.quale(x,27,c(0,0.4),18,c(10,0.4),
  8,c(0,0.4),10,c(210,1),22,c(210,1),
  7,c(-90,0.8),1,c(130,0.7),29,c(140,0.4),
  23,c(130,0.2),26,c(-140,1.2),
  9,c(-30,0.7),13,c(0,0.5),28,c(220,0.6),u)}

scambi=c(16,17,10,22,30,11,8,18,19,20,9,13)

L=Din.componenti(X,f)
X=c(L,recursive=T)

Gs.postscript("09-figura-2.ps",10,10)
lato=c(-8,8)
par(bg="#ccff99")
Gd.grafica(lato,lato,cornice=T)
Din.grafo(X,f,r=7,
  rotazioni=rota,scambi=scambi)
dev.off()