MATLAB function for calculating the transfer function of the Big Muff Pi tonestack

Below is a gist containing the function used to calculate the transfer function of the Big Muff Pi tone stack. An example of how the code is be used can be found in this gist, which was used to produce the above image.