code/figure8.m