Change the API for the initial guess

Currently the user is expected to pass along an initial guess for the coefficients used to combine the basis functions into an approximate solution to the BVP. Would be better if the user only needed to provide a list of arrays defining the initial guess for each endogenous variable. The solver could then optimally fit the basis functions to determine the initial guess for the coefficients.