You can plot the solutions from the solve command by treating them as parametric equations for a surface and a line. Note that one of the solutions listed by solve is a subset of another of the solutions.
p=parametric_plot3d((r1,-r1,r2),(r1,-4,4),(r2,-2,2))