Private Sub Command1_Click()
    'PROGRAM D6R9
    'Driver for routine BNLDEV
    N = 20: NPTS = 10000: ISCAL = 200: LLEN = 50: NN = 100
    Dim DIST(21), Text$(50)
    IDUM& = -133
    For J = 1 To 21
        DIST(J) = 0#
    Next J
    Print "  Mean of binomial distribution (0..20); -1 to end."
    XM = 5
    Print Tab(3); XM
    If XM < 0 Then End
    PP = XM / NN
    For I = 1 To NPTS
        J = Int(BNLDEV(PP, NN, IDUM&))
        If (J >= 0) And (J <= 20) Then DIST(J + 1) = DIST(J + 1) + 1
    Next I
    Print "   x      p(x)         graph:"
    For J = 1 To 20
        DIST(J) = DIST(J) / NPTS
        For K = 1 To 50
            Text$(K) = " "
        Next K
        KLIM = Int(ISCAL * DIST(J))
        If KLIM > LLEN Then KLIM = LLEN
        For K = 1 To KLIM
            Text$(K) = "*"
        Next K
        Print Tab(3); Format$(CSng(J - 1), "####.#0");
        Print Tab(10); Format$(DIST(J), "#####.###0");
        Print "   ";
        For M = 1 To 50
            Print Text$(M);
        Next M
    Next J
End Sub