Summary1. The feeding functional response is one of the most widespread mathematical frameworks in Ecology, Marine Biology, Freshwater Biology, Microbiology and related scientific fields describing the resource-dependent uptake of a consumer. Since the exact knowledge of its parameters is crucial in order to predict, for example, the efficiency of biocontrol agents, population dynamics, food web structure and subsequently biodiversity, a trustful parameter estimation is of utmost importance for 1 . CC-BY-NC-ND 4.0 International license It is made available under a (which was not peer-reviewed) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity.The copyright holder for this preprint . http://dx.doi.org/10.1101/201632 doi: bioRxiv preprint first posted online Oct. 13, 2017; scientists using this framework. Classical approaches for estimating functional response parameters lack flexibility and can often only serve as approximation for a correct parameter estimation. Moreover, they do not allow to incorporate side effects such as resource growth or background mortality. Both call for a new method to be established solving these problems.2. Here, we combined ordinary differential equation models (ODE models), that were numerically solved using computer simulations, with an iterative maximum likelihood fitting approach. We compared our method to classical approaches of fitting functional responses, using data both with and without additional resource growth and mortality.3. We found that for classical functional response models, like the often used type II and type III functional response, the established fitting methods are reliable. 4. Our method will enable researchers from different scientific fields that are measuring functional responses to estimate parameters correctly. These estimates will enable community ecologists to parameterize their models more precisely, allowing for a deeper understanding of complex ecological systems, and will increase the quality of ecological prediction models.