For comparison with other CAS choose from: Axiom Derive Macsyma Maple Mathematica
solve(
( 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));
{{c1=arbcomplex(12), c2=(c1*( - c1*beta*epsilon - c1*beta*theta + alpha*epsilon + alpha*theta))/(c1*gamma*theta - epsilon*eta), c3=(c1* 2 ( - c1 *beta*gamma + c1*alpha*gamma - c1*beta*eta + alpha*eta) )/(c1*gamma*theta - epsilon*eta)}}