def solve(self, basis_kwargs, boundary_points, coefs_array, nodes, problem,<EOL>**solver_options):
raise NotImplementedError<EOL>
Solve a boundary value problem using the collocation method.

Parameters
----------
basis_kwargs : dict
    Dictionary of keyword arguments used to build basis functions.
coefs_array : numpy.ndarray
    Array of coefficients for basis functions defining the initial
    condition.
problem : bvp.TwoPointBVPLike
    A two-point boundary value problem (BVP) to solve.
solver_options : dict
    Dictionary of options to pass to the non-linear equation solver.

Return
------
solution: solutions.SolutionLike
    An instance of the SolutionLike class representing the solution to
    the two-point boundary value problem (BVP)

Notes
-----