diff(x,t) uses mathrm round the ds