Hello!
I'm trying to solve numerically an ODE system with piecewise. And this piecewise is very important for this task.
This system describes behavior of a pulley with friction. There are some constants: m, c, g, mu and J. Values of this constants are not important.
> sys := m*a(t) = piecewise(a(t) < a0, F0*time, a(t) >= a0, 0), v(t) = diff(x(t), t), a(t) = diff(v(t), t);
> m := 5; F0 := 10; a0 := 5;
> initialconditions := x(0) = 0, v(0) = 0;