code/figure5.m