#!/usr/bin/python

from numpy.random import poisson
from numpy.random import normal

npoints=20000
cp1 = 1000;
cp2 = 12000;
lambda1=1.0; 
lambda2=10.0;
lambda3=5.0

mean1=10.0;
std1=1.0;
mean2=100.0;
std2=4.0;
mean3=4.0;
std3=4.0;

for i in range(npoints): 
	if i<cp1: print poisson(lambda1)
	elif i<cp2: print poisson( lambda2 )
	else: print poisson( lambda3 )

	#if i<cp1: print normal( mean1, std1 )
	#elif i<cp2: print normal( mean2, std2 )
	#else: print normal( mean3, std3 )