let c = h *. dt *. sqrt(k)