Changed Takahashi routine for interpolation of D_AB, rather than
interpolating constants.