d13r4.txt

常用的数值算法的VB程序

Private Sub Command1_Click()
    'PROGRAM D13R4
    'Driver for routine TTEST
    NPTS = 1024
    MPTS = 512
    EPS = 0.02
    NSHFT = 10
    Dim DATA1(1024), DATA2(512)
    'Generate Gaussian distributed data
    IDUM& = -5
    For I = 1 To NPTS
        DATA1(I) = GASDEV(IDUM&)
    Next I
    IDUM = -11
    For I = 1 To MPTS
        DATA2(I) = (NSHFT / 2#) * EPS + GASDEV(IDUM&)
    Next I
    Print
    Print Tab(5); "Shift      T     Probability"
    Print
    For I = 1 To NSHFT + 1
        Call TTEST(DATA1(), NPTS, DATA2(), MPTS, T, PROB)
        SHIFT = (I - 1) * EPS
        Print Tab(5); Format$(SHIFT, "#.#0"); Tab(14); Format$(T, "#.#0");
        Print Tab(24); Format$(PROB, "#.#0")
        For J = 1 To NPTS
            DATA1(J) = DATA1(J) + EPS
        Next J
    Next I
End Sub