gexian.m

来自「Matlab 编辑的 高斯 二分法 割线法 牛顿法」· M 代码 · 共 28 行

 %割线法
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;