@staticmethod
def _rearrange_ode(x, y, coeff_b, fx):
    """Rearrange coefficients for scipy solver.

        Parameters
        ----------
        x : np.ndarray(N,)
            Points from desired domain
        y : np.ndarray(order, N)
            Initial guess for the function values and its derivatives
        coeff_b : list[number] or np.ndarray(Order + 1)
            Coefficients for each differential part on ODE
        fx : np.ndarray(N,)
            Non-homogeneous term at given x

        Returns
        -------
        np.ndarray(N,)
            proper expr for the right side of the ODE equation
        """
    result = fx
    for i, b in enumerate(coeff_b[:-1]):
        result -= b * y[i]
    return result / coeff_b[-1]
