f=@(x)x.^2-x-1;
 a = 1;
 b = 2;
 tol=0.001;
 fa=f(a);
 fb=f(b);
  p=False_postion_2019_dej(f,fa,fb,a,b,tol);