Newbie trying to solve a 3D moving heat source problem with crank-Nicholson for a month

Revised in 3D moving heat source problem with crank-Nicholson ODE solver algorithm.