edf function

Compute empirical distribution function