diffusion constant:
points:
dx:
dt: