code/figure3.m