>> [x,y,z]=solve('x+y=r*z,y=x/y+1,z=y/x+1','x,y,z') x = r ry = 1/2+1/2*(4*r+1)^(1/2) 1/2-1/2*(4*r+1)^(1/2)z = 1/2/r*(2*r+1+(4*r+1)^(1/2)) 1/2/r*(2*r+1-(4*r+1)^(1/2))>> ezplot('z=1/2/r*(2*r+1+(4*r+1)^(1/2))')