In this manuscript we present B-POP, a MATLAB toolbox for bi-level optimization through multiparametric programming. It features i) bi-level programming solvers for linear and quadratic programming problems, and their mixed-integer counter-parts, ii) a versatile problem generator capable of creating random bi-level problems of arbitrary size, and iii) a library of bi-level programming test problems. The features of B-POP are demonstrated through detailed computational studies showing the capabilities and the scalability of the embedded algorithms. Moreover, two applications, i) a supply chain planning problem, and ii) a hierarchical model predictive control of a reactor system are chosen to show the applicability of bi-level programming and B-POP. (C) 2018 Elsevier Ltd. All rights reserved.