diff_ek function

Differentiate x and y and compute change amplitude via an algorithm proposed by Engbert & Kliegl (2003)