Branch made to modify the first half of the integrand of zeta in file ver2.0/latm/SRC_1setinput/integrands.f90. Only kept the spiMatElem(ix,icp,ic) of line 1475 to check convergence in response zeta_bulk (41).