Singular BVP must be on interval [0, b] with 0 < b. Singular term must be square matrix of size neqn.