iterate1.m

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

%递推计算的稳定性
function y = iterate1(n, a)
%method1: I(n) = -a*I(n-1) + 1/n, n = 1,2, ..., 10
%I(0) = log( (a+1)/a ), a = 15 | a = 0.05
if nargin == 2 ,%只有输入参数为2个时才计算
    if n == 0,
        y = log((a+1)/a);
    end
    if n > 0,
        y = -a * iterate1(n - 1, a) + 1/n;
    end
end