logtick.exp <- function(x, n, iax, laball) {
    xs <- rep(NA, times = n)
    xs[1] <- x
    for (i in 2:n) {
        xs[i] <- xs[i-1]*10
    }
    xs2 <- rep(NA, times = n)
    xs2[1] <- 2*x
    for (i in 2:n) {
        xs2[i] <- xs2[i-1]*10
    }
    xs9 <- rep(NA, times = n)
    xs9[1] <- 9*x
    for (i in 2:n) {
        xs9[i] <- xs9[i-1]*10
    }
    vec2 <- seq(xs2[1], xs9[1], xs[1])
    for (i in 2:n) {
        vec2 <- c(vec2, seq(xs2[i], xs9[i], xs[i]))
    }
    for (i in iax) {
        axis(i, at= log(xs), lab = xs)
    }
    for (i in iax) {
        if (! laball[i]) {
            axis(i, at = log(vec2), lab = NA, tcl = -0.2)
        }
        else {
            axis(i, at = log(vec2), lab = vec2, tcl = -0.2,
                 cex.axis = 0.6, mgp = c(2.1, 0.2, 0))
        }
    }
    box(bty = "l")
    return()
}
