Total: 1
This work proposes a higher-order iterative framework for solving matrix equations, inspired by the structure and functionality of neural networks. A modification of the classical Jacobi method is introduced to compute higher-order coefficient matrices through matrix-matrix multiplications. The resulting method, termed the higher-order Jacobi method, structurally resembles a shallow linear network and allows direct computation of the inverse of the coefficient matrix. Building on this, an iterative scheme is developed that, once the network parameters are determined for a known system, enables efficient resolution of system variations without re-computing the coefficients. This iterative process naturally assumes the form of a deep recurrent neural network. The proposed approach moves beyond conventional Physics-Informed Neural Networks (PINNs) by providing an explicit, training-free definition of network parameters rooted in physical and mathematical formulations. Computational analysis demonstrates improved order of complexity, suggesting a promising direction for algorithmic efficiency in solving linear systems. This methodology opens avenues for interpretable and scalable solutions to physically motivated problems in computational science.