Sra.incroci = function (MP,k,R,X)
{n=nrow(MP)
for (j in seq(1,40,2))
{P1=MP[,j]; P2=MP[,j+1]
p=runif(1,0,0.5)
for (i in 1:n) if (runif(1)<p)
{a=P1[i]; P1[i]=P2[i]; P2[i]=a}
R1=Sra.g(P1,X,k); R2=Sra.g(P2,X,k)
if ((R1<R[j])&&(R2<R[j+1]))
{MP[,j]=P1; MP[,j+1]=P2}}
MP}

Sra.muta = function (P,k,R)
{n=length(P); p=runif(1,0,0.5)
for (i in 1:n) if (runif(1)<p)
P[i]=Snc.interi(1,1,k)
P}

Sra.mutazioni = function (MP,k,R,X)
{for (j in 1:40)
{P=Sra.muta(MP[,j],k,R)
if (Sra.g(P,X,k)<R[j]) MP[,j]=P}
MP}

Sra.nuovi = function (MP,k,a,b)
{n=nrow(MP); for (j in a:b)
MP[,j]=Snc.interi(n,1,k)
MP}

Sra.rendimenti = function (MP,X,k)
apply(MP,2,Sra.g,X,k)

Sra.visualizza = function (t,dt,rend,P)
{P=paste(P,collapse=' ')
cat('\n',rend,': ',P,'   dopo ',t,
  ' generazioni\n',sep='')
a=readline('Vuoi continuare? ')
if (a=='n') 0
else {v=as(a,'numeric')
if (!is.na(v)) v else dt}}