用二分法求函数fx=x^2-x-1在区间(1,2)内一个零点

clearclcf=@(x)x。^2-x-1;xab=[1;2];pab=sign(f(xab));while(diff(xab)>1e-6)x=mean(xab);p=sign(f(x));if p,break;endxab(p==pab)=x;enddisp(x)结果：1。6180