You can check eaxisin packagesfsmisc
x <- seq(1, 100000, len = 10)
y <- seq(1e-5, 1e-4, len = 10)
plot(x, y)

# eaxis
plot(x, y, axes = FALSE)
eaxis(side = 1)
eaxis(side = 2)

expression pretty10exp() . , :
plot(x, y, axes = FALSE)
title(pretty10exp(y[1]))
