Das Ergebnis sieht sehr überzeugend aus,
aber der Weg dahin ist wohl eher so: (s. Kommentare)
Etwa für xo>1 :
\( |f(x) - f(x_0)| = |\frac{ x-1 }{|x^2 - 1|} - \frac{ x_0-1 }{|x_0^2 - 1|}| = |\frac{ x-1 }{|(x+1)(x - 1)|} - \frac{ x_0-1 }{|(x_0+1)(x_0 - 1)|}| \)
 \( =| \frac{1}{|x+1|} - \frac{1}{|x_0+1| } |\)  #
Für xo>1 ist ja bei hinreichend kleinem δ
nämlich δ<1  jedenfalls  x>0 und  also
|x+1| = x+1 und |xo+1|=xo+1  .
Dann setzt sich # fort mit
 \( =| \frac{1}{x+1} - \frac{1}{x_0+1} |\)
und dann Hauptnenner
\( =| \frac{x_0+1}{(x+1)(x_0+1)} - \frac{x+1}{(x+1)(x_0+1)} |\)
\( =| \frac{x_0+1 -(x+1)}{(x+1)(x_0+1)}  |\)
\( =| \frac{x_0 -x}{(x+1)(x_0+1)}  |\)
Da beide Faktoren im Nenner größer als 1 sind, ist
der Bruch also kleiner als der Zähler und du hast
\( < | x_0 -x    |\).
Mit δ=min {1 ;  ε }   [zur 1 s.o.]   hast du also jedenfalls:
\(  | x_0 -x    | < \delta  \)  ==>  \( |f(x) - f(x_0)| < \epsilon \).
Damit wäre der Fall xo>1 geklärt.
Versuche mal die anderen Fälle, vergiss auch nicht xo zwischen
-1 und +1.