%% Plot animation of FPI main code
%Dejen Ketema, March, 2019
function p=fixed_point_plot_2019(g,p0,xmin,xmax)
xx=linspace(xmin,xmax,1000);
plot(xx,g(xx),xx,xx,'linewidth',2)
axis equal,grid on,axis([xmin,xmax,xmin,xmax])

fprintf(' n         p          |p-p0|   \n');
fprintf('-------------------------------\n');
for n=1:20
  p=g(p0);
  fprintf('%2d  %12.8f  %12.8f\n',n,p,abs(p-p0));
  line([p0,p0],[p0,p],'color','r','linewidth',2); pause(1);
  line([p0,p] ,[p,p], 'color','r','linewidth',2); pause(1);
  p0=p;
end