code: def _evaluate_rhs(cls, funcs, nodes, problem):
          
          evald_funcs = cls._evaluate_functions(funcs, nodes)
          evald_rhs = problem.rhs(nodes, *evald_funcs, **problem.params)
          return evald_rhs
docs: 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)