Just some functions

library(lattice)

tau <- 2 * pi
n <- 100
x <- y <- seq(-2 * tau, 2 * tau, by=tau/32)
m <- matrix(rep(NA, n * n), nrow=n)

k <- c(1/4, 1/2, 1, 2, 4, 8, 16, 32)
for(k in k) {
    for(i in 1:n) {
        for(j in 1:n) {
            m[i,j] <- cos(k*x[i])^2 * -tan(y[j])
        }
    }
    
#     png(paste('images/test', k, '.png', sep=''), width=1612, height=806)
    g <- wireframe(m, shade = TRUE,
              main='I really dont know what is it',
              xlab='Adorable function 1',
              ylab='Neat function 2',
              zlab='Tada!')
    plot(g)
#     dev.off()
}