code/figure7E.m