diedai_jaitken.m

来自「非线性方程组的求解方法：使用matlab编写」· M 代码 · 共 18 行

    function [x]=DieDai_aitken(x0,epsilon)
    k=1;x(1)=f1(x0);
    while 0<1
    xx(k,1)=f1(x(k));xx(k,2)=f1(xx(k,1));
    x(k+1)=xx(k,2)-(xx(k,2)-xx(k,1))^2/(xx(k,2)-2*xx(k,1)+x(k));
    if abs(x(k+1)-x(k))<epsilon
        return
    end
    x(k)=x(k+1);
    k=k+1;
    end
    end
    %  [x]=DieDai_aitken(1,0.000000001)
    % function f=f1(x)
    % f=x^3+2*x^2+10*x-20;
    % end
    % x =1.369255425256832   1.368808110910821   1.368808107821373
    % 1.368808107821373   1.368808107821373