On Wed, Jul 3, 2013 at 9:54 AM, Thilina Rathnayake <[email protected]> wrote: > > > Putting n = 0, 1, .... will result in solutions in the four classes. > > Putting n = 0 in the second one will give (40, 11), which corresponds to > a fundamental solution for on of the classes. Putting n = 0 in the fourth > one will give (12, 3) > which also corresponds to one of the fundamental solution we have found. > > Putting n= 0 in the first one and third one does not yield the same > fundamental > solutions found by diop_pell() (they simplifies to (-40, -11) and (-12, -3) > respectively). > I used the LMM algorithm discussed in the paper, so sometimes fundamental > solutions > found by different algorithms for each class may differ (LMM returns minimal > positive > solutions so this makes sense). > > Here is the paper, > > http://www.jpr2718.org/pell.pdf > > x**2 - 13*y**2 = 27. diop_pell(13, 27) is discussed at the last > paragraph of page 14.
I see! Yes, so I think that the only part missing is how to get from the general equivalence class, for example (12, 3), to the general solution in terms of "n" as Mathematica returns. Do you know how to implement that? Ondrej -- You received this message because you are subscribed to the Google Groups "sympy" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. To post to this group, send email to [email protected]. Visit this group at http://groups.google.com/group/sympy. For more options, visit https://groups.google.com/groups/opt_out.
