How to compute Ergodic channel capcity as a function of number transmit antenna

How can I compute the channel capacity C as a function of K (K changes from 2 to 16) averaged over 1000 iterations based on the equation. C = E[log2(det(I_M + roh*H'DH))] where D = (1/K)I_K, (I_K is an identity matrix of dimension KxK) roh = 10dB, I_M is an identity matrix = I_K. C was supposed to be a row vector of size 1x16 but I was getting a 1x1 dimensional matrix.

clear; close all; clc;
M = 2:16;
K = 2:16;
N = 1000;
roh = 10;

for ii = K
    I_M = eye(ii);
    D = (1/ii) * I_M; 
    H = 1/sqrt(2) * (randn(ii,ii) + 1i * randn(ii, ii));
    C = mean(log2(det(I_M + roh * H'* D * H)));

      C_sum = 0;
      for k = 1:N
          C_sum = C_sum + C;
      end
  end

Answer 1

[alhossin aljadai]

Slamoalycom bro Abdullahi, if its bossible I need to contact u regarding a matlab code which Im orking on these days,but stuk with it. my email is [e-mail] thank you in advance.