matlab误差传播和算法稳定性

算法描述:

    

方案二：

递推公式结果：

y(1)=0.212647           

y(2)=0.071838           

y(3)=0.065374           

y(4)=0.046157           

y(5)=0.038461           

y(6)=0.032051           

y(7)=0.027701           

y(8)=0.024325           

y(9)=0.021697           

y(10)=0.019576           

y(11)=0.017833           

y(12)=0.016375           

y(13)=0.015137           

y(14)=0.014073           

y(15)=0.013148           

y(16)=0.012338           

y(17)=0.011621           

y(18)=0.010984           

y(19)=0.010412           

y(20)=0.009897 

积分结果显示：

积分代码：

syms x
y1=x/(1+4*x);
y2=int(y1,0,1);
y2
double(y2)

递推公式：

y0 = 1/4 - log(5)/16;
y1 = 0;
n= 1;
while n<=20
    y1 = 1/4*(1/n-y0);
    fprintf('y(%d)=%-20f',n,y1);
    y0 = y1;
    n = n+1;
    if mod(n,1) == 0
        fprintf('\n');
    end
end