Whether the solver Settings are reasonable

[nate]

@ss780060872

for n in range(num_steps):
    ...
    solve(F == 0, T, bcs)
    plot(T)
    ...

Your code has inefficiencies. See, e.g. Out of memory error in frequency for loop - #4 by nate.