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