2:- use_module(library(mcintyre)).    3
    4:- if(current_predicate(use_rendering/1)).    5:- use_rendering(c3).    6:- use_rendering(graphviz).    7:- endif.    8
    9:- mc.   10
   11:- begin_lpad.   12
   13a(X):binomial(X,20,0.5).
   14
   15:- end_lpad.

?- mc_expectation(a(X),1000,X,E).