Abhinav , you do know that MonteCarlo is generation of random set of #s (based on area under the curve integration ).

If you have the equations of outage capacity then we I can look at it and tell you how to write , since I do not have communication toolbox

Dear Abhinav 

Look at this simple matlab program to determine MIMO channel capacity with Monte Carlo  realization.  , i hope that can be help you to solve your problem.  

clear

clc

SNR_dB=[0:12];

SNR_linear=10.^(SNR_dB/10);

nT=3; nR=3; % 3x3 Tx and Rx antennas

n=min(nT,nR); % determine the rank of channel matrix

I = eye(n);

N_trial= 100; % Monte Carlo (Number of realization)

C = zeros(N_trial,length(SNR_dB));

for k=1:N_trial;

H = randn(nT,nR);

HH=abs(H*H');

for l=1:length(SNR_dB)

C(k,l) = C(k,l)+log2(real(det(I+(SNR_linear(l)/nT)*HH)));

end

end

C_m=mean(C);

% Plots

figure

plot(SNR_dB,C_m,'r-o');

hold on,

xlabel('SNR[dB]'); ylabel('bps/Hz');

grid on

title('MIMO channel capacity nT=nR=3')

I agree with Diyari Hassan's answer

https://www.osapublishing.org/oe/viewmedia.cfm?uri=oe-24-1-219&seq=0

yet