I need to solve a system of differential equations and I have a function that returns dx/dt and du/dt, but I need to add the constraint that u is in the interval [0,1]. How can