  function P4()
% Checks out AllRoots
alpha = 33;
z = AllRoots(alpha);
x = linspace(-1,1,300);
y = sin(alpha*x)-x;
plot(x,y,'m',x,0*x,'k',z,0*z,'.r','Markersize',20)
shg
clc
fprintf('alpha = %10.6f\n\nThe roots of f(x) = sin(alpha*x)-x...\n\n',alpha)
fprintf('%20.15f\n',z)
   
