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)
. */