Enforce Kinematic Constraint in Initial Value Problem ode45() MATLAB I have the following nonlinear system and accompanying ODE: My solution in MATLAB: