圆周率之高精度运算

zopey

本例是入门级的（10年前就有的代码） ，运算要点：
1、傅里叶级数，展开式有很多 选个收敛快的
2、数组 存储结果 ，低位到高位
3、模拟 手动四则运算的过程，  step by step
4、计算速度一般

圆周率之高精度运算.rar (11.67 KB, 下载次数: 56)

2019-6-20 21:28 上传

点击文件名下载附件

2015年 Chudnovsky(丘德诺夫斯基)方法 是计算PI非常快速并且精度很高的方法， 计算三千两百万位 6分钟。

数学家太牛了。

zopey

1. Sub getpi()
   
2. Rows("2:65535").ClearContents
   
3. 
   
4. Dim k0, nums, max As Long, f(), g, i As Long, j As Long
   
5. nums = [d1]
   
6. k0 = [b1]
   
7. max = 18 * nums
   
8. ReDim f(0 To max)
   
9. 
   
10. g = k0 / 5
    
11. For i = 0 To max
    
12.     f(i) = g
    
13. Next
    
14. 
    
15. 'Stop
    
16. '''''''''''''''
    
17. Dim nRow, t, d As Long, k As Long, arr()
    
18. ReDim arr(0 To Int(nums / 20), 1 To 20)
    
19. 
    
20. For j = max To 1 Step -18
    
21. '采用普遍的傅立叶级数展开方法 pi=2+1/3*(2+2/5*(2+3/7*(2+4/9*(2+5/11  .....))))))
    
22. t = 0
    
23. For i = j To 1 Step -1
    
24.     t = t + f(i) * k0
    
25.     d = 2 * i + 1
    
26.     f(i) = t - Int(t / d) * d
    
27.     t = Int(t / d) * i
    
28. Next i
    
29. k = k + 1
    
30. If k = 21 Then k = 1: nRow = nRow + 1
    
31. 
    
32. arr(nRow, k) = Int(g + t / k0) Mod k0
    
33. g = t Mod k0
    
34. Next j
    
35. 
    
36. [a2].Resize(Int(nums / 20) + 1, 20) = arr
    
37. [a2].Resize(Int(nums / 20) + 1, 20).NumberFormatLocal = Right(k0, Len(k0) - 1)
    
38. End Sub

复制代码