method2.m

来自「Matlab做的数值计算实验」· M 代码 · 共 11 行

function y = method2(n)
%利用method1 Sn' = Sn - (Sn - Sn-1)^2/(Sn-2*Sn-1＋Sn-2), n = 3,4,...
%则Ln2 ＝ Sn'
if nargin == 0
    fprintf('Not enough parameter for method2');
elseif n < 3
    n = 3;
    y = method1(n) - (method1(n) - method1(n-1))^2 / (method1(n) - 2*method1(n-1) + method1(n-2));
else
    y = method1(n) - (method1(n) - method1(n-1))^2 / (method1(n) - 2*method1(n-1) + method1(n-2));
end