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.