d13r6a.txt

常用的数值算法的VB程序

Private Sub Command1_Click()
    'PROGRAM D13R6A
    'Driver for routine TUTEST
    NPTS = 5000
    MPTS = 1000
    EPS = 0.02
    NSHFT = 10
    VAR1 = 1#
    VAR2 = 4#
    Dim DATA1(5000), DATA2(1000)
    'Generate two Gaussian distributions of different variance
    IDUM& = -51773
    FCTR1 = Sqr(VAR1)
    For I = 1 To NPTS
        DATA1(I) = FCTR1 * GASDEV(IDUM&)
    Next I
    FCTR2 = Sqr(VAR2)
    For I = 1 To MPTS
        DATA2(I) = (NSHFT / 2#) * EPS + FCTR2 * GASDEV(IDUM&)
    Next I
    Print
    Print Tab(5); "Distribution #1 : variance =  "; Format$(VAR1, "#.00")
    Print Tab(5); "Distribution #2 : variance =  "; Format$(VAR2, "#.00")
    Print
    Print Tab(5); "Shift       T     Probability"
    For I = 1 To NSHFT + 1
        Call TUTEST(DATA1, NPTS, DATA2, MPTS, T, PROB)
        SHIFT = (I - 1) * EPS
        Print Tab(5); Format$(SHIFT, "#.#0");
        Print Tab(15); Format$(T, "#.#0");
        Print Tab(25); Format$(PROB, "#.#0")
        For J = 1 To NPTS
            DATA1(J) = DATA1(J) + EPS
      Next J
    Next I
End Sub