% scripts can be obtained from https://github.com/zhou-feng/matlab_scripts ...