MATLAB code for Estimation of the optical communication link BER