龙贝格算法.frm

完整的数值模拟软件

VERSION 5.00
Begin VB.Form Form19 
   BackColor       =   &H00FFFFFF&
   Caption         =   "龙贝格算法"
   ClientHeight    =   9210
   ClientLeft      =   60
   ClientTop       =   345
   ClientWidth     =   13365
   LinkTopic       =   "Form19"
   ScaleHeight     =   9210
   ScaleWidth      =   13365
   StartUpPosition =   3  '窗口缺省
   Begin VB.CommandButton Command1 
      Caption         =   "龙贝格算法"
      BeginProperty Font 
         Name            =   "宋体"
         Size            =   14.25
         Charset         =   134
         Weight          =   700
         Underline       =   0   'False
         Italic          =   0   'False
         Strikethrough   =   0   'False
      EndProperty
      Height          =   615
      Left            =   5280
      TabIndex        =   0
      Top             =   7320
      Width           =   2535
   End
End
Attribute VB_Name = "Form19"
Attribute VB_GlobalNameSpace = False
Attribute VB_Creatable = False
Attribute VB_PredeclaredId = True
Attribute VB_Exposed = False
Private Sub Command1_Click()
Dim a, b, g As Integer
'a = 0
'b = 1

a = InputBox("输入积分下限a的值", , 100)
b = InputBox("输入积分上限b的值", , 100)
Cls
Print
EPS = 10 ^ (-10)
'g = InputBox("输入g得值", , 20)
h = b - a
Dim T(500000), s(500000), c(500000), r(500000)
T(1) = h * (4 / (1 + a * a) + 4 / (1 + b * b)) / 2
'T(1)=h*
k = 1
e = T(1)
Print "----------------------------------------------------------------------------------------------------------------------------------"

Print "分半次数", "区间等分数", "     T", , "       S", , "        C", , "      R"

Print "----------------------------------------------------------------------------------------------------------------------------------"

Print " 0", " 1", T(1),
While e > EPS
   Print
   n = 2 ^ (k - 1)
   ss = 0
   h = (b - a) / n
   i = 1
   While i <= n
      X = a + (2 * i - 1) * h / 2
      ss = ss + 4 / (1 + X * X)
      i = i + 1
   Wend
   T(2 ^ k) = T(2 ^ (k - 1)) / 2 + ss * h / 2
   s(2 ^ (k - 1)) = T(2 ^ k) * 4 / 3 - T(2 ^ (k - 1)) / 3
   If k = 1 Then
      Print k, 2 ^ k, T(2 ^ k), , s(2 ^ (k - 1)),
   Else
      Print k, 2 ^ k, T(2 ^ k), s(2 ^ (k - 1)),
   End If
   If k >= 2 Then
      c(2 ^ (k - 2)) = (16 * s(2 ^ (k - 1)) - s(2 ^ (k - 2))) / 15
      Print c(2 ^ (k - 2)),
   End If
   If k >= 3 Then
      r(2 ^ (k - 3)) = c(2 ^ (k - 2)) * 64 / 63 - c(2 ^ (k - 3)) / 63
      Print r(2 ^ (k - 3)),
   End If
   
   e = Abs(T(2 ^ (k - 1)) - T(2 ^ k))
   r1 = T(2 ^ k)

   k = k + 1
Wend
Print
Print
Print "                 r=", r1
End Sub

Private Sub Label2_Click()

End Sub

Private Sub Command2_Click()

End Sub

Private Sub Form_Load()

End Sub