% LMA400, 2017, Föreläsning 2, Exempel 1
f = @(x) (x.^2 - 1)./(x.^3 - x.^2 + x - 1);
x = 0:0.1:2;
plot(x, f(x), 'b', 'linewidth', 2)
grid on
shg