I have been asked to solve a diffusion equation with a source term using finite differences method. I need to numerically integrate the following equation either in MATLAB or C++. The equation is dT/dt = d2T/dx2 + S(x) The form of S(x) is some function given by a Gaussian profile. Could anyone have a solution to the problem! Thanks!