file 02-11-2014 fix error in simulation, pbeta =1