Given a generic expression of a conic:
$$Ax^2 + Bxy + Cy^2 + Dx + Ey + F=0$$
is there a way to write an expression for one of the branches as a function of the coefficients? I tried using the quadratic formula to get an expression for $y$: $$y=\frac{-(Bx+E)\pm \sqrt{(Bx+E)^2 - 4(C)(Ax^2 + Dx + F)}}{2C}$$
but this doesn't always work. Consider:
$$xy=1$$
Here, $A=0, B=1, C=0, D=0, E=0, F=0$, so $y=\frac{\cdot}{0}$, which isn't particularly helpful. In other cases it is not as bad, but still not what I'm looking for. E.g.
$$x^2 - y^2 - 1=0$$
Using the formula above, we get $$y=\pm \sqrt{x^2-1}$$ which seems nice, but $y=\sqrt{x^2-1}$ it is actually one half of each branch rather than one entire branch, as can be seen here.
http://www.wolframalpha.com/input/?i=plot%28x^2+-+y^2+-+1%3D0%29
http://www.wolframalpha.com/input/?i=plot%28y%3Dsqrt%28x^2-1%29%29
I am trying to draw one of these branches, so I need an ordered set of points along a predefined "grid" of either of the variables. Is it possible to do this?


I think you are going to have a more satisfactory experience doing this: find the center $\vec{P} = (x_0, y_0)$ of your hyperbola. Find the eigenvectors of the matrix $$ \left( \begin{array}{cc} A & B/2 \\ B/2 & C \end{array} \right) $$ and normalize and choose order and $\pm$ so that with basis $\vec{u},\vec{v}$ you can then write your branch as $$ g(t) = \vec{P} + \vec{u} \cosh t + \vec{v} \sinh t. $$ Note that while $\vec{u},\vec{v}$ are perpendicular to each other they are not particularly of any length, in the way I write it above. If you prefer an orthonormal basis you then just put one constant scalar multiplication in front of the $\cosh$ term and one in front of the $\sinh$ term. It will all work out if you actually do have a hyperbola, which happens when $B^2 > 4 A C.$
I did teach this about 20 years ago, a whole section on translations and rotations of conic sections for a engineering calculus course. The memories are dim, but it does appear that it is best to find an orthonormal basis for the above matrix first, then find the center expressed in that basis, and end up with constant scalar coefficients of the hyperbolic trig terms.