why is this in `interpolate.py` and not in `polyint`? Same question lagrange, actually. A minor detail, sure.