d7f1.frm

来自「常用数值算法集哈」· FRM 代码 · 共 88 行

VERSION 5.00
Begin VB.Form Form1 
   Caption         =   "Form1"
   ClientHeight    =   5640
   ClientLeft      =   60
   ClientTop       =   345
   ClientWidth     =   6750
   LinkTopic       =   "Form1"
   ScaleHeight     =   5640
   ScaleWidth      =   6750
   StartUpPosition =   3  'Windows Default
   Begin VB.CommandButton Command1 
      Caption         =   "Command1"
      Height          =   375
      Left            =   5160
      TabIndex        =   0
      Top             =   4920
      Width           =   1335
   End
End
Attribute VB_Name = "Form1"
Attribute VB_GlobalNameSpace = False
Attribute VB_Creatable = False
Attribute VB_PredeclaredId = True
Attribute VB_Exposed = False
Private Sub Command1_Click()
     'PROGRAM D7F1
     'Driver for routine GASDEV
     Dim TEXT(50) As String
     N = 20
     NP1 = N + 1
     NOVER2 = N / 2
     NPTS = 10000
     ISCAL = 400
     LLEN = 50
     Dim DIST(21)
     IDUM = -13
     For J = 1 To NP1
         DIST(J) = 0#
     Next J
     For I = 1 To NPTS
         J = CInt(0.25 * N * GASDEV(IDUM)) + NOVER2 + 1
         If (J >= 1) And (J <= NP1) Then DIST(J) = DIST(J) + 1
     Next I
     Print
     Print Tab(5); "Normally distributed deviate of  "; NPTS; " points"
     Print
     Print Tab(5); " x         p(x)       graph:"
     For J = 1 To NP1
         DIST(J) = DIST(J) / NPTS
         For K = 1 To 50
             TEXT(K) = ""
         Next K
         KLIM = CInt(ISCAL * DIST(J))
         If KLIM > LLEN Then KLIM = LLEN
         For K = 1 To KLIM
             TEXT(K) = "*"
         Next K
         X = J / (0.25 * N)
         Print Tab(5); Format$(X, "#.#0");
         Print Tab(15); Format$(DIST(J), "0.###0");
         Print Tab(22); "     ";
         For K = 1 To 50
             Print TEXT(K);
         Next K
         Print
    Next J
End Sub
Function GASDEV(IDUM)
    Static ISET, GSET
    If ISET = 0 Then
      Do
        Randomize
        V1 = 2# * Rnd - 1#
        V2 = 2# * Rnd - 1#
        R = V1 ^ 2 + V2 ^ 2
      Loop While R >= 1# Or R = 0
      FAC = Sqr(-2# * Log(R) / R)
      GSET = V1 * FAC
      GASDEV = V2 * FAC
      ISET = 1
    Else
      GASDEV = GSET
      ISET = 0
    End If
End Function