The psi function for Tukey’s biweight estimator
The analytic derivative of rho
1d array
ndarray
psi(z) = z*(1 - (z/c)**2)**2 for |z| <= R
psi(z) = 0 for |z| > R
TukeyBiweight.psi