display('log')
check('Test/Surgical_model.txt')
data('Test/Surgical_data.txt')
compile(2)
inits(1, 'Test/Surgical_inits.txt')
inits(2, 'Test/Surgical_inits2.txt')
update(500)
set(p[])
set(alpha)
set(beta)
update(1000)
stats(*)
history(*)
density(*)
autoC(*)
quantiles(*)
save('SurgicalLog')