# playing along with Elias Wegert in R: X <-…

[This article was first published on

Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.

**isomorphismes**, and kindly contributed to R-bloggers]. (You can report issue about the content on this page here)Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.

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 input ring <- function(C) C[.8 < Mod(C) & Mod(C) < 1.2] #focus on the unit circle lev <- 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 their blog:**isomorphismes**.R-bloggers.com offers

**daily e-mail updates**about R news and tutorials about learning R and many other topics. Click here if you're looking to post or find an R/data-science job.

Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.