XProfan

Título de la ventana upper$("Binomialverteilungsdichte: P(n Versuche, k veces Erfolg,"+\
" en Einzelerfolgs-Wahrscheinlichkeit p)")
'(CL) CopyLeft 2018-04 by P.Pájaro carpintero, Wien. OHNE GEWÄHR!
Ventana de Estilo 24:Ventana %maxx/4,10 - %maxx/2,%maxy-60

Proc Binom_nCr :parámetros N&,k&

    caso k&=0:volver 1.0:caso N&<=0:volver 0.0:caso (2*k&)>N&:k&=N&-k&
    var P!=1:whileloop k&:P!=P!*(N&-k&+&Loop)/&Loop:endwhile:volver P!

ENDPROC

Proc Fac :parámetros N&:casenot N&:volver 1.0

    if N&>100000:Imprimir " *** FAC() TOO BIG! *** ";:beep:volver 0.0:endif

        if N&<171:var p!=1:whileloop N&:p!=p!*&Loop:endwhile:volver p!

            más :var s!=0:whileloop N&:s!=s!-lg(&Loop):endwhile:set("decimals",17)
            volver str$( round(10^(int(s!)-s!),13))+"e+"+str$(int(-s!) ):endif

        ENDPROC

        Proc BinVtlg :parámetros N&,k&,p!

            caso (n&<=0) or (k&>n&) or (p!<=0) (p!>=1):volver (k&=0)
            caso (p!>=1):volver (k&>0)
            volver Binom_nCr(N&,k&)*p!^k&*(1-p!)^(n&-k&)

        ENDPROC

        Declarar p!,n&,k&,su!
        Nochmal:
        CLS:font 2
        imprimir "\n Nötige Vorinformation:"
        imprimir "\n           Erfolgswahrscheinlickeit uno Einzelversuchs:":font 0
        imprimir "\n     (z.B. Münze p = 50 %; ´6´= 1 Würfelseite de 6 Möglichen: p = 16.6667 %)"
        font 2:imprimir "\n p_Einzel [%] ?: ";
        input p!:caso (p!<0) or (p!>100):goto "Nochmal"

        REPEAT

            Imprimir "\n    Gesamtzahl Versuche [0=Restart]?: ";:input n&:caso n&=0:BREAK
            imprimir   "      Zufällige Anzahl a Erfolgen ?: ";:input k&
            imprimir   " Wahrscheinlichkeit para genau esta Anzahl: "; 100*BinVtlg(N&,k&,p!/100);" [%]"

            if (N&>1) and (N&<1000):imprimir   "             ... para esta Anzahl y mehr: ";:su!=0

                Whileloop k&,N&:su!=su!+100*BinVtlg(N&,&Loop,p!/100):endwhile:caso su!>100:su!=100

                    imprimir su!,"[%]"

                endif

            UNTIL 0

            goto "Nochmal"