code: def _compute_residuals(self, coefs_array, basis_kwargs, boundary_points, nodes, problem):
          
          coefs_list = self._array_to_list(coefs_array, problem.number_odes)
          derivs, funcs = self._construct_approximation(basis_kwargs, coefs_list)
          resids = self._assess_approximation(boundary_points, derivs, funcs,
                                              nodes, problem)
          return resids
docs: Return collocation residuals.
  
          Parameters
          ----------
          coefs_array : numpy.ndarray
          basis_kwargs : dict
          problem : TwoPointBVPLike
  
          Returns
          -------
          resids : numpy.ndarray