# Variante di numeri_complessi_1.R, nella quale definiamo
# una funzione per disegnare puntino e segmento
# e una per creare il plot iniziale

# 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)

z = (0+1i)^(1/4); for (k in 1:16) plot.im(z^k, lwd=2, col='blue')
z = 1.05*(0+1i)^(1/4); for (k in 1:16) plot.im(z^k, lwd=2, col='red' )
z = 0.95*(0+1i)^(1/4); for (k in 1:16) plot.im(z^k, lwd=3, col='green' )

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