mu = 0.99;



data_set = zeros(10,2);

for i = 1:10

    [x y] = poissonsim(mu,100);

    data_set(i,:) = [mean(x) sqrt(var(x))];

end