modelDisplay('window') modelPrecision(6) modelCheck('MhHugginsmod.txt') modelData('MhHugginsdat.txt') modelCompile(3) modelInits('MhHugginsint1.txt', 1) modelInits('MhHugginsint2.txt', 2) modelInits('MhHugginsint3.txt', 2) modelGenInits() samplesSet(N) samplesSet(b) samplesSet(g) samplesSet(pie) modelUpdate(100000, TRUE) samplesStats("*") samplesHistory("*") samplesDensity("*") modelSaveLog('HugginsMhLog') samplesCoda(*, MhHugginscoda) modelQuit()