Initial guessing to bvp4c

I am working on a 4th order non-linear variable coefficient homogeneous ODE bvp. I am having issues getting a solution using bvp4c. This could be one of many things. Not having a solution within the boundaries I am providing, due to the non-linearity matlab cannot resolve to a solution, or the initial guesses for the functions are so far off it can't work with them. The error I am getting is

Error using bvp4c (line 251)
Unable to solve the collocation equations -- a singular Jacobian encountered.

I was wondering if there is a better way to provide initial guesses. I would like to give a vector of points that the solution should be near. Is that possible?

btw, I am unable to post example an example code.

The input for a initial guess is as follows

solinit = bvpinit(xinit,yinit)

Typically the xinit is a vector of points that range in the boundary of the problem (x-axis).

As for yinit, in a 2-point bvp, you will typically see [# #], and the numbers coorespond to the BC.

How can I use yinit to apply mulitple point conditions? For my problem I have 4 BC, do I need something that looks like this [# # # #]?

I guess I'm confused about how the yinit works. I've read on mathworks and stuff but there is limited information.

Thanks for the help in advance.