solving system of equations



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