MATLAB code in my research