untitledshuzhiwending.m

数值计算实践是我们必备的一门学科！这是我在学习之余利用MATLAB编写的程序

fprintf('数值稳定算法;\n')
y0=0.5*(1/(6*21)+1/(5*21));   %给出初值y0;
n=20;
plot(n,y0,'r*');              %画出n与y0的图示;
hold on
fprintf('y[20]=%10.6f\n',y0);   %输出y的值;
while(1)
    y1=0.2*(1/n-y0);          %给出通过y0计算出的y1的值;
    fprintf('y[%2.0f]=%10.6f',n-1,y1);
    plot(n-1,y1,'r*')
    if(n<=1)break;end
    y0=y1;
    n=n-1;
    if mod(n,2)==0,fprintf('\n'),end    %每两个值就换行;
end