Create diffusion1d.m