Thus we have an approximation to a differential equation in the function OilSize[t] which we proceed to solve.

Input := 

eqDe = OilSize'[t] == pfit[OilSize[t]];