用fortran或者matlab计算1+1/4+1/9……+1/10000000000,使其误差小于10^-6.发现溢出导致没有结果怎么办？

用fortran或者matlab计算1+1/4+1/9……+1/10000000000,使其误差小于10^-6.发现溢出导致没有结果怎么办？

用MATLAB没有溢出的问题啊？
>> a=sum(1./(1:1e5).^2)
a =
          1.64492406689824
>> b=symsum(1/k^2,1,1e5);
>> double(b)
ans =
          1.64492406689823
>> double(b)-a
ans =
    -1.59872115546023e-014以上分别用两种方法求级数和，得到的结果相差很小（1e-14量级）。其中b为符号解，理论上可以得到任意精度的结果。例如，如果希望保留30位有效数字，可以用
>> vpa(b,30) 
ans = 
1.64492406689822626980574850332

不要用数组， 直接做计算，数据太多了，很难保证误差那么小。 不过可以试下 ans = 0.d0 do i = 1, 100000000000 ans = 1.d0/dble(i) enddo 再看看别人怎么说的。