I have settled this question here; in brief, a FORTRAN implementation of algorithms for solving this problem have been published in the Collected Algorithms of the ACM.
J. M. isn't a mathematician
- 4.1k
- 31
- 50