d6r5.txt

VB常用数值算法集 内含有解线性代数方程组

Private Sub Command1_Click()
    'PROGRAM D6R5
    'Driver for routine EXPDEV
    NPTS = 10000: EE = 2.718281828
    Dim TRIG(21), X(21)
    For I = 1 To 21
        TRIG(I) = (I - 1) / 20#
        X(I) = 0#
    Next I
    IDUM& = -1
    For I = 1 To NPTS
        Y = EXPDEV(IDUM&)
        For J = 2 To 21
            If (Y < TRIG(J)) And (Y > TRIG(J - 1)) Then
                X(J) = X(J) + 1#
            End If
        Next J
    Next I
    TOTAL = 0#
    For I = 2 To 21
        TOTAL = TOTAL + X(I)
    Next I
    Print "Exponential distribution with "; NPTS; " points:"
    Print "      interval     observed    expected"
    For I = 2 To 21
        X(I) = X(I) / TOTAL
        EXPECT = Exp(-(TRIG(I - 1) + TRIG(I)) / 2#)
        EXPECT = EXPECT * 0.05 * EE / (EE - 1#)
        Print Tab(5); Format$(TRIG(I - 1), ".#0"); Tab(12); Format$(TRIG(I), ".#0");
        Print Tab(21); Format$(X(I), ".###0"); Tab(33); Format$(EXPECT, ".###0")
    Next I
End Sub