let t_j j = ((float j) *. dt)