Dioxins and aflatoxins contamination along the dairy supply chain. We used a mixed integer nonlinear programming approach to maximize the performance of the sampling in terms of reducing the risk of the potential disability adjusted life years (DALYs) in the population. Decision variables are the number of samples collected and analyzed at each stage of the food chain (feed mills, dairy farms, milk trucks, and dairy processing plants) for each chemical, given a predefined budget. The model was applied to the case of monitoring for aflatoxin B1/M1(AFB1/M1) and dioxins in a hypothetical Dutch dairy supply chain, and results were calculated for various contamination scenarios defined in terms of contamination fraction and concentrations.