matlab 稳定性,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