Sum over a dimension

Suppose I have a V-by-K A matrix, a K-by-K B matrix, and K-by-N C matrix. How should I sum over k from 1 to K specified summation(k=1..K)of product (A_vk*B_kk*C_kn).

Answer 1

[Roger Stafford]

Bd = diag(B);
D = zeros(V,N);
for iv = 1:V
  for in = 1:N
    D(iv,in) = sum((A(iv,1:K).').*Bd(1:K).*C(1:K,in));
  end
end