% ----------------------------------------------------
% Step One - Evaluation of the required constant terms
% ----------------------------------------------------

% (Sigma_m)^2 [sm2]
ebno=linspace(0,30,200);
erx0=1;
erxMUI=ones(1,5);
fc=1e11;
tm=5e-9;
B=4e9;
pulse=waveform(fc,tm,B);
Rb=20e6;
gamma_r=1;

sm2 = cp0901_sm2_OOK(pulse,fc);
  
% MUI energy summation [EMUI]

EMUI = sum(erxMUI./erx0);
  
% -------------------------------------
% Step Two - Evaluation of SIR and SNRn
% -------------------------------------

SIR = (gamma_r)/(sm2*Rb*EMUI);

EBN0 = 10.^(ebno/10);
SNRn = EBN0.*2;

% -----------------------------------
% Step Three - Performance evaluation
% -----------------------------------

SNRref = 1./((1./SNRn)+(1/SIR));

BER = 0.5.*erfc(sqrt(SNRref./2));
semilogy(ebno,BER);
axis([0 30  10e-6 1e0]);