Private Sub Command1_Click()
    'PROGRAM D6R7
    'Driver for routine GAMDEV
    N = 20: NPTS = 10000: ISCAL = 200: LLEN = 50
    Dim DIST(21), Text$(50)
    IDUM& = -13
    Do
        For J = 1 To 21
            DIST(J) = 0#
        Next J
        Print "  Order of Gamma distribution (n=1..20); -1 to end."
        IA = 4
        Print Tab(3); IA
        If IA <= 0 Then End
    Loop While IA > 20
    For I = 1 To NPTS
        J = Int(GAMDEV(IA, IDUM&)) + 1
        If (J >= 1) And (J <= 21) Then DIST(J) = DIST(J) + 1
    Next I
    Print "  Gamma-distribution deviate, order "; IA; "of "; NPTS; " points"
    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), "####.#0");
        Print Tab(10); Format$(DIST(J), "#####.###0");
        Print "   ";
        For M = 1 To 50
            Print Text$(M);
        Next M
    Next J
End Sub