[BHE] Using central differences to compute the jacobian for nonlinear solver of 'basic_newton'

Using central differences to compute the jacobian in the BHE process, so that the nonlinear solver of 'basic_newton' can be alternatively used.

The ctest has also been added, in which the 'basic_newton' is used to solve the original benchmark: BHE_1P.