display('log')
check('model.txt')
data('data.txt')
compile(1)
inits(1, 'inits.txt')
gen.inits()
update(500)
set(mu)
set(tau)
set(sigma)
update(10000)
stats(*)
history(*)
density(*)
save('output')