В вашем коде вы просто берете средние из двух последовательных чисел из p
и сохраняете их в mu_p
, что является совершенно неверным выводом, который вы получаете.
Для среднего значения всех чиселиз индекса [1 - n] я бы предложил следующее.
mu_p = cumsum(p) % built-in function for cumulative sum.
denom = 1:n_r % denominator for each index in mu_p = # terms till the index
mu_p = mu_p ./ denom % elementwise division
Пример
Ввод: p = [1,2,3,4,5,6,7,8,9];
Вывод: mu_p = 1.0000 1.5000 2.0000 2.5000 3.0000 3.5000 4.0000 4.5000 5.0000