数值稳定matlab.txt

舍入误差与数值稳定性 数值稳定算法 数值不稳定算法

y0=(1/205+1/246)/2
n=40;
fprintf('y[%d]=%f\n',n,y0);
while(1)
    y1=1/(5*n)-y0/5;
    fprintf('y[%d]=%f\n',n-1,y1);
    if(n<=1)break;
    end
    y0=y1;
    n=n-1;
end