modified MATLAB script CPML 1D program