计算方法的MATLAB实现.docVIP

误差基础 % ex01.m clear; digit = 3; S0 = log(6)-log(5); S0 = vpa(S0,digit); S(1) = 1 - 5*S0; S0 = vpa(S(1),digit); for n = 2 : 8 S(n) = 1/n - 5*S(n-1); S(n) = vpa(S(n),digit); end S=double(S); fprintf( S0 = %11f\n,double(S0)); for n=1:8 fprintf( S%1d = %11f\n,n,S(n)); end ____________________________ % ex02.m clear; digit = 3; S(8) = (1/45 + 1/54)/2; S(8) = vpa(S(8),digit); for n = 7 : -1: 1 S(n) = 1/(5*(n+1)) - S(n+1)/5; S(n) = vpa(S(n),digit); end S0 = 1/5 - S(1)/5; S0 = vpa(S0,digit); S=double(S); fprintf( S0 = %11f\n,double(S0)); for n=1:8 fprintf( S%1d = %11f\n,n,S(