fixing another bug. Jacobi should be defined as a matrix.
the local equilibrium reaction system is running, but not producing right results yet.