This project is all about discretizing 2D steady state heat conduction equation using finite difference method by employing a mesh which was developed to have growth rates in both direction based on user input like seen in ANSYS solver. This can also be run without the use of growth rate meshes. A mesh convergence study was also done to understand the importance of grids in solving such problems.