t <- sort(runif(20, 0, pi))
x <- t**2 # ok
x <- c(sin(t[1:10]), cos(t[11:20])) # not ok

length(boxplot.stats(diff(x), coef=10)$out) > 0

splinefun(t, x)(t, 3) # obtains the jerk of x with respect to t