Implicit form - transient case

Objective: Crank Nicolson method implementation `(delT)/(dt) = alpha((del^2T)/(delx^2) + (del^2T)/(dely^2))` `(T_(i,j)^(n+1) - T_(i,j)^n)/(Deltat) = alpha(1/2(T_(i+1,j)^(n+1) - 2T_(i,j)^(n+1) + T_(i-1,j)^(n+1) + T_(i,j+1)^n - 2T_(i,j)^n + T_(i,j-1)^n))/(Deltax^2) + alpha(1/2(T_(i+1,j)^(n+1) - 2T_(i,j)^(n+1) + T_(i-1,j)^(n+1) + T_(i,j+1)^n…