Compute the value of the right-hand side of the system of ODEs.

        Parameters
        ----------
        basis_funcs : list(function)
        nodes : numpy.ndarray
        problem : TwoPointBVPLike

        Returns
        -------
        evaluated_rhs : list(float)

def _evaluate_rhs(cls, funcs, nodes, problem):
        """
        Compute the value of the right-hand side of the system of ODEs.

        Parameters
        ----------
        basis_funcs : list(function)
        nodes : numpy.ndarray
        problem : TwoPointBVPLike

        Returns
        -------
        evaluated_rhs : list(float)

        """
        evald_funcs = cls._evaluate_functions(funcs, nodes)
        evald_rhs = problem.rhs(nodes, *evald_funcs, **problem.params)
        return evald_rhs