

if evalResults(end,sodaPar.iterCol)==sodaPar.nSamples+sodaPar.nOffspring

    subplotscreen(1,2,1,'figureNumber',1)
    subplotscreen(1,2,2,'figureNumber',2)
    
end


figure(1)
plotgelmanrubin(sodaPar,critGelRub)

figure(2)
marghist(sodaPar,evalResults)

drawnow
