Some drafts for figure 2
Important addition: function fPCA.m (operates like radial_dist.m, but for times = [50 200])