Seems like this could be possible with IntervalConstraintProgramming + ModelingToolkit and setting a large enough IntervalBox to cover all feasible values, unless I’ve misunderstood how IntervalConstraintProgramming works. But it would require ModelingToolkit v3, which conflicts with Symbolics.
I am Weaving solutions to mathematical problems from here (some plotting, algebra, diffeqs and a bit of optimization to come) into a pdf for students at my university, and trying to showcase how cool Julia is in this domain. So ideally I need a solution that doesn’t involve uninstalling and reinstalling key packages throughout.