Just in case somebody is going to proceed with this thread: The article from which @Laiju apparently copied the algorithm can be found here (with much better readable formulas).
In principle, a SAS program implementing the algorithm would consist of only a few lines of code. Most of the terms in the formulas seem to be clearly defined. Moreover, the article contains a table of calculated values (Table 2, p. 245) for comparison. However, I was unable to replicate these values. I even tried to solve the equation for unclear parameters numerically (given a few values from Table 2; without sophisticated tools such as SAS/OR, though), but didn't obtain satisfactory results.
... View more