Partial Coherence Analysis in MATLAB

I have collected 3 acceleration vs time signals and 1 sound pressure level vs time signal. I have created a multi input,single output system, where the accleration signals are the input and sound pressure level signal is the output. I want to compute partial coherence against frequency between every input and output. Partial coherence is different from the typically calculated coherence. I have written a code for this on MATLAB but I am not sure whether it is correct or not theoretically. Therefore, I would highly appreciate if anyone can review my code. Thank you.

Moreover, I was wondering, if there is a code in MATLAB or Python, that I can use for partial coherence analysis for vibration analysis. Thank you.