How to plot this euation?

Here is the code what i wrote but it does not give proper picture

:

clear all;
  clc;
  No=10^-6;
  Po=10;
  Rb=16000;
  Tb=1/Rb;
  Eb=Po.*Tb;
  Da=5;
  Da_db=10^(Da/10);
  Nc=32;
  Mc=4;
  L=8;
  K=1;
  M=3;
  j=0;

    Meu1=sqrt(1./(1+(No./(2.*Eb))+(2/(3.*Da_db.*Nc)).*((1+(Mc./5))*L*K-1)));

    % fact1=factorial(M-1+j);
    % fact2=factorial(j);
    % fact3=factorial(M-1);

    n=M-1+j;
    fact=factorial(n)/(factorial(j).*factorial(n-j));

    sum=0; 
    for i=1:M;
        for j=0:2;
         s(i)=fact.*(((1+Meu1)/2)^j);
         sum=sum+s(i);
         ber(i)=(((1-Meu1)/2)^M).*s(i)
        end
    end

    semilogy(Eb/No,ber)
    axis([0 10 10e-12 10e-1])

    xlabel('Eb/No,dB');
    ylabel('BER');
    grid on