x = linspace(-2, 3, 501);
y = x.^2 - 1;
plot(x, y);