Finite Difference Method to solve Heat Diffusion Equation in Two Dimensions.

Heat diffusion equation of the form Ut=a(Uxx+Uyy) is solved numerically. All units are arbitrary.

This code employs finite difference scheme to solve 2-D heat equation. A heated patch at the center of the computation domain of arbitrary value 1000 is the initial condition. Bottom wall is initialized at 100 arbitrary units and is the boundary condition. As the algorithm marches in time, heat diffusion is illustrated using a movie function at every 50th time step. Code also indicates, if solution reaches steady state within predetermined number of iterations. All the units are arbitrary.