Solver for the heat equation in an square domain using FVM. It includes:

• A solver with a source term at the middle of the domain and no-flux BC
• A solver with a source term at the middle of the domain, no-flux BC and variable coefficients.
• A solver with a source term at the middle of the domain and non-homogeneous BC