form1.frm

计算无理数PI前一千位的VB程序

VERSION 5.00
Begin VB.Form Form1 
   Caption         =   "Form1"
   ClientHeight    =   5550
   ClientLeft      =   60
   ClientTop       =   450
   ClientWidth     =   7530
   LinkTopic       =   "Form1"
   ScaleHeight     =   5550
   ScaleWidth      =   7530
   StartUpPosition =   3  'Windows Default
End
Attribute VB_Name = "Form1"
Attribute VB_GlobalNameSpace = False
Attribute VB_Creatable = False
Attribute VB_PredeclaredId = True
Attribute VB_Exposed = False
Private Sub Form_Click()
Dim a(0 To 1005) As Long, b(0 To 1005) As Long, c(0 To 1006) As Long, sa As Long, sb As Long, sc As Long
Dim d As Long, g As Long, k As Long, q As Long, i As Long, n As Long, j As Long, m As Long
Const w = 1003
Const r = 10
b(1) = 800: c(1) = 9560
q = 1: i = 1
Do While k <= w
   If c(m) = 0 Then m = m + 1
   For j = m To w
       sc = sc * r + c(j)
       c(j) = sc \ 57121
       sc = sc - c(j) * 57121
   Next j
   If b(k) = 0 Then k = k + 1
   For j = k To w
       sa = b(j) + sa * r
       b(j) = sa \ 25
       sa = sa - b(j) * 25
       d = b(j) - c(j) + sb * r
       sb = d \ i
       a(j) = a(j) + q * sb
       sb = d - sb * i
    Next j
    i = i + 2
    q = -q
    sa = 0: sb = 0: sc = 0
    n = 0
    For j = w To 1 Step -1
      g = a(j) + n
      If (g >= 0) Then
         n = g \ r: a(j) = g - n * r
      Else
         n = (g - 9) \ r
         a(j) = g - r * n
      End If
    Next j
Loop
Print "3.";
For j = 1 To w - 1
   Print Trim(Str(a(j)));
   If (j Mod 50 = 0) Then Print
   
Next j
Print
Print m
End Sub