1:- use_module(library(mcintyre)).    2
    3:- if(current_predicate(use_rendering/1)).    4:- use_rendering(c3).    5:- endif.    6:- mc.    7:- begin_lpad.    8
    9m(S,N,P):multinomial(S,N,P).
   10
   11query(N,P):-
   12    m(S,N,P),
   13    S=[N1,N2|_],
   14    N1>N2.
   15    
   16:- end_lpad.   17
   18prob_query(N,P,Prob):-
   19  mc_sample(query(N,P),1000,Prob)