MATLAB implementation of quasi-anharmonic analysis