n=1000000;
print("number of samples is "); print(n);
sprintf("n= %d \n",n)

lambda=3.1

print("    k    P(X>k) ");

for (k in 30:50){
count=0;

for (i in 1:n) {
samp=rpois(15,lambda)
total=sum(samp);
if (total<=k) count<-count+1;
}
prob<-count/n;
print(k);
print(prob);

sprintf("k=%d prob=%f ",k,prob);
# print(ppois(k,15*lambda))

}