moving jacobian to a new subroutine
version converges for neq=2, nvar=1, for mu_e var