Initialize the matrix and vector to calculate the equilibrium concentrations
Addresses #12 (closed). This PR also changes how numElements is calculated, moving to the convention in the Fortran code where it is one less than the value in the input file (see #46).