Help finding multiple feedback (MFB) high-pass filter design formulas?
Since you mention scripting, it might be better to start from the bottom, up: find the roots of the lowpass prototype, which can be transformed and scaled, then applied to whatever transfer function you have.
Currently, what you want is to implement a transfer function for each topology you find; that's going to hurt in the long run. But, if you say you have the lowpass transfer function then you can apply the RC-CR transformations (as mentioned in the comments). This involves replacing each R with x/C, and each C with x/R, where x=1/\$\omega_0\$. So, the lowpass can be transformed into highpass like this:

For the lowpass (bottom half), V(o11) is the MFB topology, while V(o12) is the Laplace expression of the transfer function. For the highpass (upper half), V(o21) is the MFB topology, while V(o22) is the Laplace t.f.. The values for the components are slightly bogus, except C2 which has been calculated for a 10 kHz corner frequency (for better alignment when plotting). It's fairly obvious the ones containing lp represent the lowpass and the ones containing hp represent the highpass.
The resulting transfer function is compared with the textbook t.f., seen in E4. The traces match, but the Laplace expressions have been slightly shifted to avoid overlapping (like the phases do), for better viewing.
From the math side, the lowpass and the highpass t.f. are these:
$$\begin{align}
H_{LP}(s)&amp;=-\frac{\frac{1}{R_1R_3C_1C_2}}{s^2+\frac{1}{C_1}\left(\frac{1}{R_1}+\frac{1}{R_2}+\frac{1}{R_3}\right)s+\frac{1}{R_2R_3C_1C_2}} \tag{1} \\
H_{HP}(s)&amp;=-\frac{\frac{C_1}{C_2}s^2}{s^2+\frac{1}{R_2}\left(\frac{C_1}{C_2C_3}+\frac{1}{C_2}+\frac{1}{C_3}\right)s+\frac{1}{R_1R_2C_2C_3}} \tag{2}
\end{align}$$
Apply the substitutions (I'm writing \$\omega_0\$ to see that it's there, but consider it 1, for a clearer picture):
$$\begin{align}
\omega_0&amp;=\sqrt{\frac{1}{R_2R_3C_1C_2}} \\
C'_1&amp;=\frac{1}{\omega_0R_1},\;C'_2=\frac{1}{\omega_0R_2},\;C'_3=\frac{1}{\omega_0R_3},\;R'_1=\frac{1}{\omega_0C_1},\;R'_2=\frac{1}{\omega_0C_2} \\
H_{HP}(s)&amp;\stackrel{C'_1,C'_2,C'_3,R'_1,R'_2}{=}H_{LP}(\frac1s) \\
&amp;\Rightarrow \\
H_{HP}(s)&amp;=-\frac{C'_1C'_3R'_1R'_2}{\frac{1}{s^2}+R'_1\left(C'_1+C'_2+C'_3\right)\frac1s+C'_2C'_3R'_1R'_2} \\
&amp;=-\frac{C'_1C'_3R'_1R'_2s^2}{1+R'_1(C'_1+C'_2+C'_3)s+C'_2C'_3R'_1R'_2s^2} \\
&amp;=-\frac{\frac{C'_1C'_3R'_1R'_2}{C'_2C'_3R'_1R'_2}s^2}{s^2+\frac{R'_1(C'_1+C'_2+C'_3)}{C'_2C'_3R'_1R'_2}s+\frac{1}{C'_2C'_3R'_1R'_2}} \\
&amp;=-\frac{\frac{C_1}{C_2}s^2}{s^2+\frac{1}{R'_2}\left(\frac{C'_1}{C'_2C'_3}+\frac{1}{C'_2}+\frac{1}{C'_3}\right)+\frac{1}{R'_1R'_2C'_2C'_3}}\,==\,H_{HP}(s)
\end{align}$$
If you need this in terms of the generic transfer function (omitting the '):
$$\begin{align}
H(s)&amp;=\frac{Ks^2}{s^2+\frac{\omega_0}{Q}s+\omega_0^2} \\
&amp;\Rightarrow \\
\omega_0&amp;=\sqrt{\frac{1}{R_1R_2C_2C_3}} \tag{3}\\
\frac{\omega_0}{Q}&amp;=\frac{1}{R_2}\left(\frac{C_1}{C_2C_3}+\frac{1}{C_2}+\frac{1}{C_3}\right) \tag{4}\\
&amp;\Rightarrow \\
Q&amp;=\sqrt{\frac{R_2}{R_1}}\frac{\sqrt{C_2C_3}}{C_1+C_2+C_3}
\end{align}$$
TLDR: you can use the design equations from the lowpass MFB to generate the highpass MFB.