log_tick_marks <- function(min,max)
{
    nsplit <- abs(round(log10(max-min)))
    i <- 0
    nurange <- c()
    while(i<=nsplit) {
        nurange <- c(nurange,sapply(1:10,function(x) x*(10^i)))
        i <- i+1;
    }
    nurange
}

axis(2,log_tick_marks(1,100000),FALSE,TRUE,NA,NA)