Create diffusion_sim.m