Hi I want to solve system of euations: z*z+c=z 2*z=1 algsys([2*z=1, z*z+c=z], [c]); does not work I do : 1. eliminate z: (%i6) e1:z*z+c-z; (%o6) z^2-z+c (%i7) e2:2*z-1; (%o7) 2*z-1 (%i8) e3:eliminate([e1,e2],[z]); 2. solve e3 with respect to c (%i14) solve([e3[1]], [c]); (%o14) [c=1/4] Is it the best way ? Adam