On 10/9/06, Stavros Macrakis <macrakis at gmail.com> wrote: > On the other hand, > > solve([a*x+b*y=1,a^2=a,b^2=b],[a,b,x,y]) > > nicely gets you > > [[a = 1, b = 0, x = 1, y = %r5], > [a = 0, b = 1, x = %r6, y = 1], > [a = 1, b = 1, x = %r7, y = 1-%r7]] I stand corrected. best Robert