Hallo Samira,
Inhomogene DGL 1. Ordnung  (Variation der Konstanten):
y '  +  f(x) · y  =  s(x)         [ bei dir:  y ' + 1/x ·y = ln(x)/x  ]
Mit der – im Folgenden hergeleiteten – untenstehenden  Formel (#)  kann man diese DGL direkt lösen.
Herleitung der Lösungsformel:
Die allgemeine Lösung der homogenen  DGL  y '  +  f(x) · y  =  0  ist
             [ bei dir: f(x) = 1/x ]
yh  =  c · e - F(x)        mit  F(x)  = ∫ f(x)  dx    (eine beliebige Stammfunktion von f)  
             [ bei dir:  F(x) = ln(x) ]
Die allgemeine Lösung y  der inhomogenen DGL erhält man aus yh , wenn man die Konstante c variabel macht, indem man sie durch c(x) ersetzt:
y  = c(x) · e - F(x)    ( jetzt muss nur noch c(x) bestimmt werden)
Wenn man jetzt y ableitet, hat man 
y' = c'(x) · e - F(x) +  c(x) · (-f(x)) · e - F(x) 
 y' und y in die inhomogene DGL eingesetzt ergibt:
 c'(x) · e - F(x) +  c(x) · (-f(x)) · e - F(x) + f(x) · c(x) · e - F(x)  =  s(x)   [ bei dir:  s(x) = ln(x)/x ]
                          (hebt sich auf)
und nach c‘(x) auflöst:    c'(x)  =  s(x) · eF(x)    
→  c(x)  =  ∫ c'(x) dx  = ∫ s(x) · eF(x) dx    
und damit
y  =  e- F(x) • ∫ s(x) · eF(x) dx   (allgemeine Lösung  der DGL)    (#)
Bei dir mit F(x) = ∫ f(x) dx = ln(x)  und    s(x) = ln(x)/x :
              y = ·e-ln(x) • ∫ ln(x)/x · eln(x) dx   =  1/x • ∫ ln(x) dx
              y =  1/x · (x · ln(x) - x + c )   =  ln(x) - 1 +  c/x   
Gruß Wolfgang