### Finding the primitive solutions of the diophantine equation ax^{2}+bxy+cy^{2}=n, where b^{2}-4ac < 0, a > 0, n > 0.

This is a somewhat different method to the one in in L.E. Dickson's *Introduction to the theory of numbers* and uses a substitution due to Lagrange in 1770.

See note.

The program is a BCMath version of a BC program `posrep(a,b,c,n,e)` .
E = 1 is verbose.

*Last modified 16th March 2015*

Return to main page