d13r9.txt

常用的数值算法的VB程序

Private Sub Command1_Click()
    'PROGRAM D13R9
    'Driver for routine CHSTWO
    NPTS = 2000
    NBINS = 10
    Dim BINS1(10), BINS2(10)
    IDUM& = -17
    For J = 1 To NBINS
        BINS1(J) = 0#
        BINS2(J) = 0#
    Next J
    For I = 1 To NPTS
        X = EXPDEV(IDUM&)
        IBIN = Int(X * NBINS / 3#) + 1
        If IBIN <= NBINS Then BINS1(IBIN) = BINS1(IBIN) + 1#
        X = EXPDEV(IDUM&)
        IBIN = Int(X * NBINS / 3#) + 1
        If IBIN <= NBINS Then BINS2(IBIN) = BINS2(IBIN) + 1#
    Next I
    Call CHSTWO(BINS1(), BINS2(), NBINS, -1, DF, CHSQ, PROB)
    Print
    Print Tab(5); "Dataset 1     Dataset 2"
    For I = 1 To NBINS
        Print Tab(5); Format$(BINS1(I), "###.#0");
        Print Tab(19); Format$(BINS2(I), "###.#0")
    Next I
    Print
    Print Tab(5); "Chi-squared:  "; Tab(19); Format(CHSQ, ".###0E+00")
    Print Tab(5); "Probability: "; Tab(19); Format(PROB, ".###0E+00")
End Sub