西安交通大學-計算方法上機報告

1、計算方法上機報告姓 名： 名： 學 號： 號：班 級： 級：上課班級： 上課班級： 計算方法上機報告2（3）Matlab 源程序 源程序clear; %清除工作空間變量clc; %清除命令窗口命令m=inpu

2、t('請輸入有效數字的位數m='); %輸入有效數字的位數s=0; for n=0:50t=(1/16^n)*(4/(8*n+1)-2/(8*n+4)-1/(8*n+5)-1/(8*n+6));if t<=10^(-m) %判斷

3、通項與精度的關系break;endend;fprintf('需要將n值加到n=%d\n',n-1); %需要將n值加到的數值for i=n-1:-1:0t=(1/16^i)*(4/(8*i+1)-2/(8*i+4)-1/(8*i+5)-1/(8*i+6));s=s+t; %求和運算ends

4、=vpa(s,m) %控制s的精度 （4）結果與分析 ）結果與分析當保留 11 位有效數字時，需要將 n 值加到 n=7， s =3.1415926536；當保留 30 位有效數字時，需要將 n 值加到 n=22, s =3.14159265358979323846264338328。通過上面的實驗結果可以看出，通過從后往前計算，這種算法很好的保