method1.m

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

function y = method1(n)
%利用级数Ln2 ＝ 1 － 1/2 + 1/3 - 1/4 + ... = symsum((-1)^(k-1)/k, 1, inf)
%计算Ln2
if nargin == 0
    fprintf('Not enough parameter for method\n');
elseif n < 1
    n = 1;
else
    %y = symsum(sym('(-1)^(k-1) / k'), 'k', 1, n);
    k =1;
    y = 0;
    while k <= n
        y = y + (-1)^(k-1) / k;
        k = k + 1;
    end
end