> #
> # perm.mws (Permutationen abzaehlen)
> #
> with(combinat): # Zusatzpaket Kombinatorik laden
Warning, the protected name Chi has been redefined and unprotected

> n:=6; # n setzen

                                n := 6

> G:=permute(n): # Alle moeglichen Permutationen generieren
> zaehler:=0: # Zaehler zuruecksetzen
> i:=1: j:=2: k:=4: l:=5:
> for a from 1 to numbperm(n) do # Alle Permutationen durchlaufen
>     R:=G[a]; # a-te Permutation betrachten
>     if (R[i]<R[j] and R[i+1]>R[j+1] and R[k]<R[l] and R[k+1]>R[l+1])
>     then zaehler:=zaehler+1 end if; # gegebenenfalls Zaehler erhoehen
> end do:
> zaehler/numbperm(n); # Wahrscheinlichkeit berechnen und ausgeben

                                 1/9

>