code/figure10.m