Numerical solution of the heat equation with non‐linear, time derivative‐dependent source term