BVP Shooting method with free parameters

You need to make sure you add the right amount of conditions so you don’t end up with a singular Jacobian. This is something that the least squares formalism would simplify, but it doesn’t exist yet.