one3.m

来自「数值分析方法,主要里面涵盖了计算方法方面matlab程序！」· M 代码 · 共 31 行

clear
clc
%有效算法
fprintf('有效算法结果:\n')
y0=1.0/4.0*(log(5)-log(4));
fprintf('y[%d]=%f\n',0,y0);
n=1;
while(1)
    y1=1.0/(4.0*n)-1.0/4.0*y0;
    fprintf('y[%d]=%f\n',n,y1);
    if(n>=10) break;
    end
    y0=y1;
    n=n+1;
end

%无效算法
fprintf('无效算法结果:\n')
y0=3.0/55.0;
n=10;
fprintf('y[%d]=%f\n',n,y0);
while(1)
    y1=1.0/n-4.0*y0;
    fprintf('y[%d]=%f\n',n-1,y1)
    if(n<=1) break;
    end
    y0=y1;
    n=n-1;
end