In this analysis, we first performed a critical review of one-compartment models used to describe metal toxicokinetics in invertebrates and found mathematical or conceptual errors in almost all published studies. In some publications, the models used do not represent the exact solution of the underlying one-compartment differential equations; others use unrealistic assumptions about constant background metal concentration and/or zero metal concentration in uncontaminated medium. Herein we present exact solutions of two differential-equation models, one describing simple two-stage toxicokinetics (metal toxicokinetic follows the experimental phases: the uptake phase and the decontamination phase) and another that can be applied for more complex three-stage patterns (toxicokinetic pattern does not follow two phases determined by an experimenter). Using two case studies for carabids exposed via food, based on previously published data, we discuss and compare our models to those originally used to analyze the data. Our conclusion is that when metal toxicokinetic follows a one-compartment model, the exact solution of a set of differential equations should be used. The proposed models allow assimilation and elimination rates to change between toxicokinetic stages, and the three-stage model is flexible enough to fit patterns that are more complex than the classic two-stage model can handle.