Qualitative properties of different numerical methods for the inhomogeneous geometric Brownian motion