interpolate implemented but not tested
OldSickPoor code and paper downloaded