用fortran或者matlab计算1+1/4+1/9……+1/10000000000,使其误差小于10^-6.发现溢出导致没有结果怎么办？
2个回答
展开全部
用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
直接做计算，数据太多了，很难保证误差那么小。
不过可以试下
ans = 0.d0
do i = 1, 100000000000
ans = 1.d0/dble(i)
enddo
已赞过
已踩过<
评论
收起
你对这个回答的评价是？
推荐律师服务：
若未解决您的问题，请您详细描述您的问题，通过百度律临进行免费专业咨询
