display('log')
check('MSmod.txt')
data('MSdat.txt')
compile(3)
inits(1,'MSint1.txt')
inits(2,'MSint2.txt')
inits(3,'MSint3.txt')
gen.inits()
refresh(1)
set(mu)
set(sd)
set(psi)
display('window')
update(100000)
stats(*)
history(*)
density(*)
coda(*,MScoda)