function x2=GX(x0,x1)
x2=x1-ff(x1)*(x1-x0)/(ff(x1)-ff(x0));   %第一次割線選代計算
n=1;
while (abs(x1-x0)>=1e-4)&(n<=100000)
    x0=x1;
    x1=x2;
    x2=x1-ff(x1)*(x1-x0)/(ff(x1)-ff(x0));    %割線迭代求解
    n=n+1;
end