For comparison with other CAS choose from: Derive Macsyma Maple Mathematica Reduce
[ alpha * c1 - beta * c1**2 - gamma*c1*c2 + epsilon*c3, _
-gamma*c1*c2 + (epsilon+theta)*c3 -eta *c2, _
gamma*c1*c2 + eta*c2 - (epsilon+theta) * c3], _
[c3, c2, c1])
solve( _
>> Error detected within library code: system does not have a finite number of solutions You are being returned to the top level of the interpreter.