function [prb_MF prb_Decorr prb_MMSE] = mud(snr_in_dB,S,B,Rinv,R,h)

E = 1; numoferr = 0;numoferr2 = 0;
SNR = 10^(snr_in_dB/10);         % signal to noise ratio
sgma = (E/sqrt(SNR));              % sigma, standard deviation of noise
sgma2 = (E/(2*sqrt(SNR)));
N=1000;

for i=1:N
% Multi User Detection    
    r = S*B(:,i);
    r = filter(h,r)+sgma*(randn(16,1)+i*randn(16,1));    
    MF = S'*r;
    decision = decide(MF);
    err_MF(i) = calculerr(decision,B(:,i));
    
    Decorr = Rinv*MF;
    decision_D = decide(Decorr);
    err_Dec(i) = calculerr(decision_D,B(:,i));
    
    MMSE = (pinv(R+sgma*eye(10,10)))*MF;
    decision_M = decide(MMSE);
    err_MMSE(i) = calculerr(decision_M,B(:,i));
end
prb_MF = mean(err_MF);
prb_Decorr = mean(err_Dec);
prb_MMSE = mean(err_MMSE);