d6r6.txt

来自「VB常用数值算法集 内含有解线性代数方程组」· 文本 代码 · 共 43 行

Private Sub Command1_Click()
     'PROGRAM D6R6
     '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