source("fdim.R")

#fractal dimension of a square
side=400
a=matrix(0,nrow =side,ncol=side)
a[50:200,50:200]=1
image(a,col=c("white","black"))
fdim(a,1)

#fractal dimension of a line
side=400
a=matrix(0,nrow =side,ncol=side)
a[50:200,100]=1
image(a,col=c("white","black"))
fdim(a,1)

#fractal dimensions of random lattices
finout=NULL
is.fractal=NULL
side=200
lowdens = 0.05
highdens=0.95
stp=0.05

rge1=seq(lowdens,highdens, by=stp)
for (limit in rge1) { 
  a=matrix(as.numeric(runif(side^2)<limit),nrow =side)
  xx=fdim(a,0) #replace zero with one to see fractal dimension plots
  #image(a,col=c("white","black"),main=limit) #delete comment to see the lattices
  #Sys.sleep(1)  #delete comment to keep the figure for one second
  finout=c(finout,as.numeric(xx[1]))
  is.fractal=c(is.fractal,as.numeric(xx[2]))
  
}
clrs = c("red","green")
plot(rge1,finout,xlab="Density",ylab="Dimension",pch=19,col=clrs[is.fractal+1])