Non-negativity constraint violated at t0. All indices specified in ODESET(''NonNegative'',idx) must satisfy 1 <= idx <= neq