1:- use_module(library(mcintyre)).    2
    3:- mc.    4
    5:- begin_lpad.    6
    7flip(h,P,N):P; flip(t,P,N):1-P.
    8bias(X): beta(X,2,5).
    9
   10next_toss(T):-
   11    bias(Bias),
   12	flip(T,Bias,0).
   13
   14previous_tosses(Tosses):-
   15    length(Tosses,NumTosses),
   16    tosses(Tosses,NumTosses).
   17
   18tosses([],0).
   19tosses([H|T],NT):-
   20    bias(Bias),
   21    flip(H,Bias,NT),
   22    NT1 is NT-1,
   23    tosses(T,NT1).
   24
   25:- end_lpad.

?-mc_lw_sample(next_toss(h), previous_tosses([h,t,h,h,h,h,h,t,h,t,h,h,h,t,h,h,h,h,h,h]), 1000,P). ?-mc_lw_expectation(bias(B), previous_tosses([h,t,h,h,h,h,h,t,h,t,h,h,h,t,h,h,h,h,h,h]), 1000,B,E). */