Matlab code for computing the co-spectra of wave perturbations (see my thesis (Pages 54-56))