method3.m

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

function y = method3(n)
%利用级数Ln2 ＝ 1/1×2 ＋ 1/2×2^2 + 1/3×2^3 + … ＝ symsum('1/(k*2^k)', 'k', 1,
%inf);
if nargin == 0
    fprintf('Not enough parameter for method3');
elseif n < 1
    n = 1;
    y = 0.5;
else
    %y = symsum(sym('1/(k*2^k)'), 'k', 1, n);
    k = 1;
    y = 0;
    while k <= n
        y = y + 1/(k*2^k);
        k = k + 1;
    end
end