# Il quadrato [0,1]x[0,1] e il suo trasformato con z -> z*(4+3i)
#
plot(c(-3,4),c(-1,7),type="n",xlab="", ylab="", asp=1)
abline(h=seq(-1,7,1),v=seq(-3,4,1),lty=3)
abline(h=0,v=0)
p <- seq(0,1,1/100)
for(n in p)for(m in p) points(m,n,pch=".",col="blue")
f <- function(m,n) (m*(1+0i)+n*(0+1i))*(4+3i)
for(n in p)for(m in p) points(Re(f(m,n)),Im(f(m,n)),pch=".",col="red")
#
# alternativa:
source("http://macosa.dima.unige.it/r.R")
PLANE(-4,4, -1,7)
x <- runif(5000); y <- runif(5000)
Dot(x,y, "blue")
f <- function(m,n) (m+n*1i)*(4+3i)
Dot(Re(f(x,y)),Im(f(x,y)), "orange")