d13r8.txt

常用的数值算法的VB程序

Private Sub Command1_Click()
    'PROGRAM D13R8
    'Driver for routine CHSONE
    NPTS = 2000
    NBINS = 10
    Dim BINS(10), EBINS(10)
    IDUM& = -15
    For J = 1 To NBINS
        BINS(J) = 0#
    Next J
    For I = 1 To NPTS
        X = EXPDEV(IDUM&)
        IBIN = Int(X * NBINS / 3#) + 1
        If IBIN <= NBINS Then BINS(IBIN) = BINS(IBIN) + 1#
    Next I
    For I = 1 To NBINS
        EBINS(I) = 3# * NPTS / NBINS * Exp(-3# * (I - 0.5) / NBINS)
    Next I
    Call CHSONE(BINS(), EBINS(), NBINS, -1, DF, CHSQ, PROB)
    Print
    Print Tab(5); "  Expected      Observed"
    For I = 1 To NBINS
        Print Tab(9); Format$(EBINS(I), "###.#0");
        Print Tab(22); Format$(BINS(I), "###.#0")
    Next I
    Print
    Print Tab(5); "Chi-squared: "; Format$(CHSQ, ".####0E+00")
    Print Tab(5); "Probability: "; Format$(PROB, ".####0E+00")
End Sub