2D heat conduction using ADI method for convection boundary

Code to solve 2D heat conduction equation using ADI method. There is convection at all boundaries. Symmetry gives other boundaries.

The code uses finite difference scheme and ADI method to solve for temperature profile of a square block. There is convection at all boundaries. The computational domain is the first quadrant with no flux boundary conditions at midplanes (of whole block) due to symmetry. The tri-diagonal matrices are solved using Thomas algorithm. Resulting temperature isolines are plotted for a particular time. The code can be modified to identify the steady state if needed.