Hi, i have to solve the 2D heat equation:
∂T/∂t = α∇^2 T = α(∂^2T/∂x^2 + ∂^2T/∂y^2)
It is given that at the 4 boundaries the T is 0. I've already written some code, but it doesn't works. Can someone review my code and help me? The error occurs at the end, when the new T are put in a matrix. The code has been added as an attachment.