% Define x
x = linspace(-2,2,1000);
% Define function
y = x + 1-2*sin(pi*x);
% Estimate roots as values less than a tolerance
tol = 0.01
ix = find(abs(y)<tol)
% Plot function
plot(x,y)
hold on
% Plot roots
plot(x(ix), y(ix), 'ro')
% Add legend
legend('F(x)', 'Roots')
