@CarlitosVillalbaGalea You can do
solO1:=subs(sol3,1-C1(t)-C2(t)-C3(t));
allvalues(solO1); #Length of output exceeds limit of 1000000
#But it seems to be OK by checking with
eval(%,params); #Assuming that params has been defined already

I replaced all floating point numbers by 3 or -3 depending on their signs. As an example, if the number 1.23456*10^12 appeared in your system it would have been replaced by 3, and -1.23456*10^12 by -3. That is what I called crude!
I have no answer to your second question.

I replaced all floating point numbers by 3 or -3 depending on their signs. As an example, if the number 1.23456*10^12 appeared in your system it would have been replaced by 3, and -1.23456*10^12 by -3. That is what I called crude!
I have no answer to your second question.

I know nothing about MapleSim, but have you looked at
?DAE
in Maple itself?

Since
restart;
factor(expand(3 - 4*cos(z) + cos(2*z))) ;
# returns 2*(cos(z)-1)^2
the inequality holds for all real z.
Thus we get an even more impressive version:
restart:
assume(z>0, beta>0,t::real);
b:= 3*beta*Pi*(3-4*cos(t)+cos(2*t))/(16*z):
-b+sqrt(b^2);

Since
restart;
factor(expand(3 - 4*cos(z) + cos(2*z))) ;
# returns 2*(cos(z)-1)^2
the inequality holds for all real z.
Thus we get an even more impressive version:
restart:
assume(z>0, beta>0,t::real);
b:= 3*beta*Pi*(3-4*cos(t)+cos(2*t))/(16*z):
-b+sqrt(b^2);

@Algorithmus I already submitted a report in the form of an SCR (= Software Change Request).

@Algorithmus I already submitted a report in the form of an SCR (= Software Change Request).

@Algorithmus Clearly a bug!

@Algorithmus Clearly a bug!

@Carl Love In Maple ln represents a one-valued function: The principal branch. So the identity is an identity modulo 2*Pi*I:
a:=-1+I: b:=-1-I;
evalc(ln(a)-ln(b)-ln(a/b));
#Result 2*Pi*I
Since the logarithm appears in exp this is no problem.

@MartinHepperle On this page at the top the menu item 'More' has 'Submit Software Change Request'.

You don't want randomness?
If you did you could use RandomTools as in this example, which also doesn't perturb 0:

Q:=proc( p,{ span::positive := 0.1 ,
             digits::posint:=4,
             number::posint:=1,
             method::name:=logarithmic})
    local q,i;
    q:=RandomTools:-Generate(
       float(range=abs(p)*(1-span)..abs(p)*(1+span),
       :-digits=digits,:-method=method),makeproc=true);
    signum(p)*[seq(q(),i=1..number)]
end proc;
Q(10);
Q(-.09,number=8,digits=5,method=uniform);

@Carl Love Thanks!

@Carl Love Thanks!