Script for figure 9
blog posts