where \(B_0\)  is  a scaling factor chosen to make \(\frac{B}{B_0}\) dimensionless.

Python code: Padé fit to \(c_0\left(B\right)\)

define Python function for Eq. \ref{eq:expand_y_32}: