Something like this:
m = 2 ; % i
n = 3 ; % j
C = rand(m,1) ;
k = rand(m,n) ;
M = rand(m,n) ;
iwant = zeros(m,1) ;
for i = 1:m
p = 1 ;
for j = 1:n
p = p*(1-k(i,j))^M(i,j) ;
end
iwant(i) = C(i)*p ;
end
Please cross check the code, I gave it quickly.
