Shamirs algorithm

To get the accumulation points, solve the inequalities

|q*b(1) - p*b(2)| < delta
|r*b(1) - p*b(3)| < delta

with delta = SQRT(b(1)/2).

An existing Java IP-Solver is used. The IP-Solver does the computation recursively, so it might raise an StackOverflowException. In this case, choose a smaller width. This value determines the size of the interval, in which the values for p are searched for in one step.

source

DES