R lang test

mandelbrot.R

xmin <- -2
xmax <- 2
ymin <- -2
ymax <- 2
size <- 1000 #点の数はsize*size
n <- 100 #発散判定の試行回数
X <- matrix(rep(seq(xmin, xmax, length=size), times=size), nrow=size)
Y <- matrix(rep(seq(xmin, xmax, length=size), times=size), nrow=size, byrow=TRUE)
C <- X + 1i * Y
Z <- matrix(rep(numeric(size), times=size), nrow = size)
Res <- matrix(rep(numeric(size), times=size), nrow = size)

for(i in 1:n) {
  Z = Z^2 + C
  mask = abs(Z) > 2
  Z[mask] <- C[mask] <- 0 + 0i
  Res[mask] <- i
  print(i)
}

image(Res)