ASKSAGE: Sage Q&A Forum - Individual question feedhttp://ask.sagemath.org/questions/Q&A Forum for SageenCopyright Sage, 2010. Some rights reserved under creative commons license.Sat, 18 Oct 2014 10:28:28 -0500Solving 3rd degree Diophantine equation in Sagehttp://ask.sagemath.org/question/24518/solving-3rd-degree-diophantine-equation-in-sage/Is there a way to find all integer solutions to $y^2=x^3-x/25+9/125$ using Sage? I tried elliptic curves but its command integral_points() won't work as the curve seems not to be an integral model.Sat, 18 Oct 2014 09:17:43 -0500http://ask.sagemath.org/question/24518/solving-3rd-degree-diophantine-equation-in-sage/Answer by vdelecroix for <p>Is there a way to find all integer solutions to $y^2=x^3-x/25+9/125$ using Sage? I tried elliptic curves but its command integral_points() won't work as the curve seems not to be an integral model.</p>
http://ask.sagemath.org/question/24518/solving-3rd-degree-diophantine-equation-in-sage/?answer=24520#post-id-24520Hello,
It seems that there is no integral solution but I did not find a Sage command to directly answers that.
First of all, **EllipticCurve** complains because the equations is not defined with integral coefficients (there are rational though). But you can still enumerate the points defined over the rationals. The curve has a group structure, and in that case you are lucky since it is rank 1, i.e. isomorphic to Z (deduced from the **E.rank()** and **E.torsion_order()** commands below).
sage: E = EllipticCurve([0,0,0,-1/25,9/125])
sage: E.rank()
1
sage: E.torsion_order()
1
You can then look at the solutions
sage: p = E.gens()[0]
sage: p
(4/25 : 33/125 : 1)
sage: 2*p
(-34319/108900 : -8297279/35937000 : 1)
sage: 3*p
(82995246844/66933451225 : -24050077094058147/17316687833675875 : 1)
And you can check that these triple are solutions of the equation:
sage: x,y,z=p
sage: y**2 == x**3 - 1/25*x + 9/125
True
sage: x,y,z = 2*p
sage: y**2 == x**3 - 1/25*x + 9/125
True
And the theory of elliptic curves tells you that these are all the rational solutions.
In order to prove that there is no integral solution, you need to prove that the denominator gets larger and larger... but I did not know how it can be automatized.
VincentSat, 18 Oct 2014 10:28:28 -0500http://ask.sagemath.org/question/24518/solving-3rd-degree-diophantine-equation-in-sage/?answer=24520#post-id-24520