def _compute_residuals(self, coefs_array, basis_kwargs, boundary_points, nodes, problem):
coefs_list = self._array_to_list(coefs_array, problem.number_odes)<EOL>derivs, funcs = self._construct_approximation(basis_kwargs, coefs_list)<EOL>resids = self._assess_approximation(boundary_points, derivs, funcs,<EOL>nodes, problem)<EOL>return resids<EOL>
Return collocation residuals.

Parameters
----------
coefs_array : numpy.ndarray
basis_kwargs : dict
problem : TwoPointBVPLike

Returns
-------
resids : numpy.ndarray