This paper introduces a novel symmetric finite difference method (SFDM) that leverages the inherent symmetry properties of the Bratu equation to efficiently compute accurate solutions and bifurcation diagrams in high dimensions, surpassing the limitations of previous methods.
