On Wed, December 20, 2006 22:25, Riccardo (Jack) Lucchetti wrote: > x += (uniform()>p) sorry, that should be x += (uniform()<p) Riccardo (Jack) Lucchetti Dipartimento di Economia Facoltà di Economia "G. Fuà" Ancona