read AppsWZ:


for i from 1 to 5 do
print(`In a Hidden-Markov Model where there is fair die`):
print(`and a loaded die (with Prob(head)=1/3, Prob(Tail)=2/3)`):
print(` and the outcome was`, 10*i , `Heads and `, 10*i , `tails`):
print(`The most reasonable(Bayesian-Laplace style not ML) estimate`):
print(`for the probability of using the fair die is`):
lu:=ExpProbFirstDie([1/2,1/2],[1/3,2/3],[10*i,10*i]):
print(lu):
print(`that in floating point is`, evalf(lu)):
od:

print(`This took`, time(), `seconds `):
quit: