The psi function for least trimmed mean
The analytic derivative of rho
1d array
ndarray
psi(z) = z for |z| <= c
psi(z) = 0 for |z| > c
TrimmedMean.psi