XProfan

Titre de la fenêtre "Hypergeometrische Verteilung"
Fenêtre Style 24:set("decimals",15):font 2
Fenêtre %maxx/10,0 - %maxx/2,%maxy-40

Proc BinCoeff :parameters N&,k&:var h!=1

    'Pascalsches Dreieck: Binärkoeffizient = ´N sur k ´

    whileloop k&,1,-1:h!=h!*(N&-&Boucle+1)/&Boucle

        endwhile:return h!

    ENDPROC

    Proc HypGFormel :parameters N&,P!,s&,x&'P(x) = (p x)*(q s-x)/(N s)

        Cas p!<0:p!=abs(p!/100):Cas (p!>=0) and (p!<1):P!=int(p!*N&)
        Retour BinCoeff(P!,x&)/BinCoeff(N&,s&) * BinCoeff(N&-P!,s&-x&)

    ENDPROC

    Proc Hyp0 :parameters N&,P!,s&'Startwert x=0 qui Rekursionsformel

        'P!=Stückzahl | si p=0..1 Wahrscheinlichkeitsfaktor, si <0: Prozentanteil
        'probabilité, dans einem Sample s aus N Elementen avec P besonderen
        'aucun qui besonderen Elemente pour trouver
        cas p!<0:p!=abs(p!/100):Cas (p!>=0) and (p!<1):P!=int(p!*N&)
        var H!=1
        var Q&=N&-P!

        WhileLoop s&,1,-1

            H!=H!*(Q&-&Boucle+1)/(N&-&Boucle+1)

        Endwhile

        Retour H!

    ENDPROC

    Proc HypGNext :parameters N&,P!,s&,x&,f_lastx!'REKURSIONSFORMEL

        Cas p!<0:p!=abs(p!/100):Cas (p!>=0) and (p!<1):P!=int(p!*N&)
        cas (x&<0) or (x&>s&):return 0
        var Q&=N&-P!
        return f_lastx!*(P!-x&)/(x&+1)*(s&-x&)/(Q&-(s&-x&)+1)

    ENDPROC

    'Main
    Déclarer f!,Addf!,w$,N&,P!,s&,x&
    Lup:
    CLS
    imprimer "\n Anzahl Grundgesamtheit: ";:input w$:cas w$>»:N&=val(w$)

    si N&<1:N&=100:P!=5:s&=20:x&=0:CLS:imprimer "\n TESTWERTE 100,5,20":goto "skip":endif

        imprimer " Anzahl Besondere Elem.: ";:input w$:cas w$>»:P!=val(w$)
        imprimer "           Sample-Taille: ";:input w$:cas w$>»:s&=val(w$)
        imprimer "    Besondere im Sample: ";:input w$:cas w$>»:x&=val(w$)
        Skip:
        imprimer "\n                      Verteilungsdichte"
        Imprimer " HyperGVtlg-Rekursionsformel:      Check avec BinCoeff-Formel:\n"

        Whileloop 0,P!

            si &Boucle : f!=HypGNext(N&,P!,s&,&Boucle-1,f!)

                d'autre     : f!=Hyp0(N&,P!,s&)

            endif

            Imprimer " x =";&Boucle;": P(x)= ";format$("%g",f!);
            imprimer tab(44);HypGFormel(N&,P!,s&,&Boucle)

        Endwhile

        Waitinput
        cls
        imprimer "\n Verteilungsfunktion en supplément:\n"
        Addf!=0

        Whileloop 0,P!

            si &Boucle : f!=HypGNext(N&,P!,s&,&Boucle-1,f!)

                d'autre     : f!=Hyp0(N&,P!,s&)

            endif

            Addf!=Addf!+f!
            Imprimer " x_max = ";&Boucle;": P(x<=";&Boucle;")= ";format$("%g",Addf!)

        Endwhile

        waitinput
        Goto "Lup"