I assume this is a simple problem, but since I don't have a strong maths background, I have trouble to implement this in matlab. Could anyone help me with it? Many thanks!
The concentration of a molecule A can be written as a convolution of of the diffusion equation solution with continuous sources. The concentration of A at position r at time t is
or simply speaking
, where f is the function of the continuous point source and h is the diffusion function.
I want to plot the concentration of A versus time t and also I need to use this equation to produce A(r,t) for a set of ODEs.