Could it be that you missed specifying values for the various constants such as h2, h3 etc, which will be needed for a numerical solve?
