dsolve[x^2y''-xy'+y==ln(x)]