Thanks Barton and Adam, I will learn the links. And about the solve function in Maxima, it seems always give simple parameterization first to x (= %r1) and y (= %r2), then solve for z. For example if I change the degree of x and y to 1, and z to 3.. it still gives x = %r1, y= %r2 and then solves z for 3 roots. Thank you ~bowo~