How to plot orbit of binary star and calculate its orbital elements?

I have a set of dates, position angles ($\theta$) and angular separations ($\rho$) for visual binary star. For example:

1994.3132  57.0    0.097
1996.1783  73.4    0.051
2006.1938 182.0    0.126
2009.2597 192.9    0.152
etc.


What are algorithms to calculate orbital elements and to plot the orbit of the binary star?

Are there available codes (or programs) for this calculations?

-
Would scicomp.stackexchange.com be a better home for this question? –  Qmechanic Jan 8 '14 at 22:46
The link has been changed to the following link: astro.gsu.edu/~gudehus/binary.html –  user1527422 May 8 at 14:30

If it's a binary, it's fairly simple (in comparison with a ternary system), because the stars of the binary orbit about their common barycenter in Kepler ellipses.

An orbit simulator, see here; replace the central star by the barycenter of the binaries.

If you don't find a ready-to-use software, try to solve it numerically by simulating the orbit of the binaries according to Kepler's law, and vary mass, distance, excentricity assumptions until they match to the observations. Use optimization methods, e.g. hill climbing algorithms, or gradient methods.

-
Or, see for example: sourceforge.net/projects/phoebe –  Moriarty Jan 9 '14 at 3:23
Needs some generalization to non-eclipsing binaries. –  Gerald Jan 9 '14 at 4:09

The good software for orbit calculation is "Binary" (Binary Star Combined Solution Package) by CHARA team:

http://www.astro.gsu.edu/~gudehus/binary.html

-