My problem can be explained with the following two plots:
plot(factorial(x),x,0,9)
and
plot(factorial(x),x,0,10)
When plotting up to 9, the plot gives me exactly what i expect, a
linear-linear plot.
When plotting up to 10 though, sage silently chose to make the
y-axis log10. This should not be the standard behavior as it leads
to confusion.