Strictly speaking `order = len(coeffs) - 3` since it only has degree > 2.
But perhaps it's better to make `a_coeff` and `b_coeff` keyword arguments and allow creating a model without passing coefficients so that it can be fit? (But leave the fitting for a future release).