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...