derivatives for lead dispersion
Add first and second derivative calculation to the lead dispersion function. The default behavior is not changed, as the derivatives are only calculated when an additional flag (deriv) is set to the order of the derivative.