Is there a way to bound the values of a variable when using scipy.integrate.solve_ivp in python?

Question

I want to solve an IVP in python with two variables, x and u, but I need the values of u to be between 0 and 1. Right now it is giving me a solution with negative values for