untitledshuzhibuwending.m

来自「数值计算实践是我们必备的一门学科！这是我在学习之余利用MATLAB编写的程序」· M 代码 · 共 13 行

fprintf('数值不稳定算法:\n')
y1=1-5*log(1.2);   %取初值y1;
n=1;
plot(n,y1,'r*')
hold on
fprintf('y[%2.0f]=%12.6f',n,y1);
for n=2:20
    y2=1/n-5*y1;
    fprintf('y[%2.0f]=%12.6f',n,y2)
    y1=y2;
    plot(n,y1,'r*')
    if mod(n,3)==0,fprintf('\n'),end
end