The code used to plot Figure 15.7 and an associated function that implements the method of Brooks, Pope and Marcolini for turbulent boundary layer trailing edge noise can be downloaded below. The BPM function was adapted from a version programmed by Remy Johnston.
BPM example