%割線法
function y=gexian(x0,x1)    

    %初始值
    x0=2;
    x1=4;
    %根據公式求x0和x1
    x2=x1-fun(x1)*(x1-x0)/(fun(x1)-fun(x0));         
    n=1;
    %迭代運算
    while (abs(x1-x0)>=1.0e-4) & (n<=100000000)  
        %將x1賦給x0
        x0=x1;     
        %將x2賦給x1
        x1=x2;      
        %重新進行計算x0和x1
        x2=x1-fun(x1)*(x1-x0)/(fun(x1)-fun(x0));     
        n=n+1;
    end
    
    %輸出結果
    y=x2
    n
   
    %定義函數
function y=fun(x)
    y=x^3+2*x^2-x-3;