# xaxt = "n" tar bort default labels för x-axeln
plot(N, p_n(N), type='l', col="red", log="x", xaxt = "n", yaxt = "n", xlab="n", ylab="t(n)", xlim=c(5*10^7, 2.8*10^11), ylim=c(0, 6.5*10^-5))
xlabels = c(10^8, 10^9, 10^10, 10^11)
ylabels = c(0, 2*10^-5, 4*10^-5, 6*10^-5)
axis(1, at=xlabels, labels=c(expression(10^8), expression(10^9), expression(10^10), expression(10^11)))
axis(2, at=ylabels, labels=c(expression(0), expression(2 %*% 10^{-5}), expression(4 %*% 10^{-5}), expression(6 %*% 10^{-5})))