Modeling of timber harvesting operation by system dynamics and the productivity estimation punction