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