We present a generalized expression for the transfer matrix of finite and infinite one-dimensional spin chains within a magnetic field with spin pair interaction J/r p , where r = 1, 2, . . . , nv is the distance between two spins, nv is the number of nearest neighbors reached by the interaction, and p ∈ [1,2]. With this generalized expression, we calculate the partition function, the Helmholtz free energy, and the specific heat for both finite and infinite ferromagnetic 1D Ising models within a zero external magnetic field. We focus on the temperature Tmax where specific heat reaches its maximum. We calculate J/(kBTmax) numerically for every values of nv ∈ {1, 2, . . . , 25}, which we interpolate and then extrapolate up to the critical temperature as a function of p, using a novel functional approach. Two different procedures are used to reach the infinite spin chain with an infinite interaction range: increasing the chain size and the interaction range by the same amount, and increasing the interaction range for the infinite chain. As we expected, both extrapolations lead to the same critical temperature, although by two different concurrent curves. Our critical temperatures as a function of p fall within the upper and lower bounds reported in the literature and show a better coincidence with many existing approximations for p close to 1 than for the p values near 2. We report an averaged inverse critical temperature J/(kBTc) = 0.532 for the one-dimensional spin chain with p = 2. It is worth mentioning that the well-known cases for near (original Ising model) and next-near neighbor interactions are recovered doing nv = 1 and nv = 2, respectively.