A dynamically consistent discretization method for Goodwin model