# Variante di numeri_complessi_2.R, nella quale
# - 'looppiamo' sulle tre condizioni iniziali
# - mettiamo una breve pausa fra una potenza e l'altra

# definizione delle funzioni ---------------------------------------------
start.plot.im <- function(Z) {
   # Z rappresenta il max di |Re(z)| e |Im(z)| rappresentabile
   plot(NULL, xlim=c(-Z,Z), ylim=c(-Z,Z), xlab='Re( )', ylab='Im( )', asp=1)
   grid()
}

plot.im <- function(z, pch=19, col='blue', lwd=1) {
   points(z, pch=pch, col=col)
   lines( c(0, Re(z)), c(0, Im(z)), col=col, lwd=lwd)
}

# Plot -------------------------------------------------------------------
salva.su.file = !TRUE
if (salva.su.file) png("potenze_nr_complessi.png")

start.plot.im(2)

rho  <- c(1, 1.05, 0.95)
cols <- c('blue', 'red', 'green')
lwds <- c(2, 2, 3)
for (i in 1:3) {               # notare il diverso significato di 'i'
    z = rho[i] * (0+1i)^(1/4)
    for (k in 1:16) {
       plot.im(z^k, lwd=lwds[i], col=cols[i])
       if(!salva.su.file) Sys.sleep(0.1)
    }
}    

if (salva.su.file) dev.off()
#-------------------------------------------------------------------------