数值不稳定matlab.txt

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

y0=log(6.0/5.0);
fprintf('y[%d]=%f\n',0.y0)
n=1;
while(1)
y1=1.0/n-5*y0;
fprintf('y[%d]=%f\n',n,y1)
if(n>=40)break;
end
y0=y1;
n=n+1;
end