f := x-> (x-1)/(x^2 - 1) ;
a := 1;

> f(a);

Error, (in f) division by zero

> plot( f(x), x=a-1..a+1);

> Limit(f(x), x=a):%=value(%);