matlab script for figures