function z=capacity(M,corr,value,XPD,alpha,output,SNR)

%**************************************************************************
%This program works with "capacity_plot_main.m" file.Figs
%2.6,2.7,2.9,2.12 are determined using this program. It operates on
%file "capacity_plot.m"
%z=capacity_plot_main(M,corr,value,ric,k,XPD,alpha,output)
%M -> number of antennas (M x M) system
%corr -> 1 if with correlation, 0 if uncorrelated (for a 2x2 system only)
%value -> correlation coefficient value from 0 ->1
%XPD -> 1 if antenna XPD is to be investigated, 0 if not (for a 2x2 system
%only)
%alpha -> XPD value
%output -> defined by 'erg' and 'out' for ergodic capacity or outage
%capacity respectively
%EXAMPLE: capacity_plot_main(2,1,0.5,1,0.5,'out')-> determines the outage
%capacity for a 2x2 system with a correlation of 0.5 and an XPD of 0.5
%**************************************************************************

%vary SNR through 20 dB

temp2=[];
for i=1:length(SNR)
    temp1(i)=capacity_snr(SNR(i),M,corr,value,XPD,alpha,output);% change this file to suit requirements e.g. capacity_rician,
                                                                %    capacity_water(for water-filling) etc.
    temp2=[temp2 temp1(i)];
    temp1(i)=0;
end
z=temp2;