dxs = xcoord[1:n] - xcoord[0: n - 1]