Hi Sha,
It would look something like this
p_values = 1:10; % Example range of p values
T_pr = @(p, i) p * i; % Example function for T_pr
P_i = ones(1, 10) / 10; % Example probabilities summing to 1
E_Tpr = zeros(size(p_values));
for idx = 1:length(p_values)
p = p_values(idx);
E_Tpr(idx) = sum(T_pr(p, 1:10) .* P_i); % Summation
end
[~, min_idx] = min(E_Tpr);
P_optimal = p_values(min_idx);
disp(['Optimal P*: ', num2str(P_optimal)]);
