My MATLAB code for EM solver by IE MOM