Solution of diffusion equation

The diffusion equation as derived earlier can be written as

(37.4)   (37.5)   (37.6)   (37.7)   (37.6) (37.7)   (37.8) (37.9)   (37.10)   (37.11) (37.12) (37.13) (37.14)

Now the Matlab functions can be written as

Function pdefun

Function icfun

Function bcfun