Solves the 2D heat equation by simulating it. I used clne to set the boundary conditions. Suggestions for better methods are welcome.