/*da interpolate on three body recombination by Steve Cota */
double da(double z);