Оператор:
E(ii,i) = randomSample([1 0],[(p(ii,i)), 1 - p(ii,i)])
определяет локальную функцию E
и не является операцией присваивания матрице E
.Используйте
E[ii,i] = randomSample([1, 0],[p[ii,i], 1 - p[ii,i]])
(я исправил дополнительные ошибки в вашем коде, поэтому, пожалуйста, проверьте различия)
и для его запуска вы также должны написать:
p = ones(n,periods)*0.5