¹small matlab script in the appendix