playing along with Elias Wegert in R: X <-…
[This article was first published on 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.
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.