Ordinary differential equation approximation of gamma distributed delay model