poissrnd(a) генерирует случайное значение Пуассона с a в качестве среднего значения. По тому, как вы определяете a, ваш вывод представляет собой серию случайных значений, которые имеют различное среднее значение. Как известно, дисперсия распределения Пуассона равна его среднему значению. Это означает, что вы вычисляете дисперсию ряда значений, которые не представляют распределение Пуассона с одним параметром (они более похожи на значения с равномерным распределением!).
Вместо этого сгенерируйте a с помощью
a = 10000 + zeros(1,10000);