Considering that Eq. (4) is a parabolic partial differential equationin one space variable x and time t, it was numerically solved
by the pdepe function in Matlab, R2010a (The MathWorks, Inc.,USA) with the aforementioned initial and boundary conditions