(This article was first published on

**isomorphismes**, and kindly contributed to R-bloggers)sine

mis-calculated the angles on sine

mid-res, small cex, pch=46, sine

low-res, small cex, pch=46

z^17th power

sine in HCL rather than HSV

cosh

z+zz+zzz+zzzz+zzzzz+zzzzzz+zzzzzzz+zzzzzzzz+zzzzzzzzz+zzzzzzzzzz

playing along with Elias Wegert in R:

X <- matrix(1:100,100,100) #grid X <- X* complex(imaginary=.05) + t(X)/20 #twist & shout X <- X - complex(real=2.5,imaginary=2.5) #recentre plot(X, col=hcl(h=55*Arg(sin(X)), c=Mod(sin(X))*40 ) , pch=46, cex=6)

Found it was useful to define these few functions:

arg<- function(z) (Arg(z)+pi)/2/pi*360 #for HCL colour inputring<- function(C) C[.8 < Mod(C) & Mod(C) < 1.2] #focus on the unit circlelev<- function(x) ceiling(log(x)) - log(x)m<- function(z) lev(Mod(z))plat<- function(domain, FUN) plot( domain, col= hcl(h=arg(FUN(domain)), l=70+m(domain)), pch=46, cex=1.5, main=substitute(FUN) ) #say it directly

NB, `hcl`

's hue ∈ `[0,360]`

so `phase`

or `arg`

needs to be matched to that.

To

**leave a comment**for the author, please follow the link and comment on his blog:**isomorphismes**.R-bloggers.com offers

**daily e-mail updates**about R news and tutorials on topics such as: visualization (ggplot2, Boxplots, maps, animation), programming (RStudio, Sweave, LaTeX, SQL, Eclipse, git, hadoop, Web Scraping) statistics (regression, PCA, time series, trading) and more...