g = plot(x + 2*ln(abs(x/2-1)), (x, 0, 8), thickness=2, ticks=[None, 1],
gridlines='minors', axes_labels=[r'$r/m$', r'$r_*/m$'],
aspect_ratio=1, ymin=-8)
g
g.save("sch_tortoise.pdf")
g = plot(exp(x)*(x-1), (x, 0, 2), thickness=2, ticks=[None, 1],
gridlines='minors', axes_labels=[r'$x$', r'$F(x)$'])
g
g.save("sch_X2mT2.pdf")