PDEs for calculation of Nodal Temperatures by GSM

Considering a Steady State 2D Heat Transfer across square walls, this is a general code that solves for the node temperature values for a square wall with specified boundary temperatures. Also, different contour plots indicate change of temperature at different nodes across the square grid.