MHB Akiva's question at Yahoo Answers regarding a Diophantine equation

  • Thread starter Thread starter MarkFL
  • Start date Start date
AI Thread Summary
The Diophantine equation X^2 + 1 = 2Y^2 has known solutions such as (1, 1), (7, 5), and (41, 29), with a recursive relationship established for generating further pairs. The discussion highlights the connection between this equation and Pell's equation, as well as its relationship to Pythagorean triples where one leg is one unit longer than the other. A recursive formula for x and y is proposed, leading to a closed-form solution that can generate an infinite number of pairs. The derived formulas allow for the calculation of any nth solution pair (x_n, y_n).
MarkFL
Gold Member
MHB
Messages
13,284
Reaction score
12
Here is the question:

X^2 + 1 = 2y^2 Diophantine equation?

I've been thinking about what whole numbers x and y satisfy:
x^2 + 1 = 2y^2

I know that (1, 1) works, as do (7, 5) and (41, 29), but I want a formula that gives me the next one, or any given one. What's the pattern?

The question is equivalent to asking which triangular numbers are twice other ones (the first three are (0, 0), (3, 2) (3rd is 6, 2nd is 3), and (20, 14)), and which right triangles with integer legs have one leg 1 unit longer than the other (first three triplets are (0, 1, 1), (3, 4, 5), (20, 21, 29)).

I'm not sure how to solve it. Which numbers x and y satisfy it? Is there an equation to find, say, the 20th x and y? Are there an infinite number of pairs that satisfy it, or just a few? Please help.

I have posted a link there to this topic so the OP can see my work.
 
Mathematics news on Phys.org
Hello Akiva,

The Diophantine equation you are considering is closely related to Pell's equation.

You have made the astute observation that the solutions to this equation are related to Pythagorean triples where the legs differ by 1. You may find http://www.mathhelpboards.com/f49/pell-sequence-2905/ to be of interest.

Okay, you have found the solutions:

$(x,y)=(1,1),\,(7,5),\,(41,29)$

At this point, we could hypothesize that $x$ and $y$ can be given recursively with:

$$S_{n+1}=6S_{n}-S_{n-1}$$

which leads to the characteristic equation:

$$r^2-6r+1=0$$

which has the roots:

$$r=3\pm2\sqrt{2}=\left(1\pm\sqrt{2} \right)^2$$

and so we know the closed form is:

$$S_n=k_1\left(1-\sqrt{2} \right)^{2n}+k_2\left(1+\sqrt{2} \right)^{2n}$$

Next, we may using the initial values for each sequence to determine the parameters.

For $x_n$, we may use:

$$x_1=k_1\left(1-\sqrt{2} \right)^{2}+k_2\left(1+\sqrt{2} \right)^{2}=1$$

$$x_2=k_1\left(1-\sqrt{2} \right)^{4}+k_2\left(1+\sqrt{2} \right)^{4}=7$$

This gives rise to the system:

$$3\left(k_1+k_2 \right)+2\sqrt{2}\left(k_2-k_1 \right)=1$$

$$17\left(k_1+k_2 \right)+12\sqrt{2}\left(k_2-k_1 \right)=7$$

Multiplying the first equation by $-6$, then adding, we get:

$$-k_1-k_2=1\,\therefore\,k_1=-1-k_2$$

Substituting into the first equation for $k_1$, there results:

$$-3+2\sqrt{2}\left(2k_2+1 \right)=1$$

Solving for $k_2$, we find:

$$k_2=\frac{\sqrt{2}-1}{2}$$

Thus:

$$k_1=-1-\frac{\sqrt{2}-1}{2}=-\frac{\sqrt{2}+1}{2}$$

and so we find:

$$x_n=-\frac{\sqrt{2}+1}{2}\left(1-\sqrt{2} \right)^{2n}+\frac{\sqrt{2}-1}{2}\left(1+\sqrt{2} \right)^{2n}$$

$$x_n=\frac{1}{2}\left(\left(\sqrt{2}-1 \right)\left(1+\sqrt{2} \right)^{2n}-\left(\sqrt{2}+1 \right)\left(1-\sqrt{2} \right)^{2n} \right)$$

For $y_n$, we may use:

$$y_1=k_1\left(1-\sqrt{2} \right)^{2}+k_2\left(1+\sqrt{2} \right)^{2}=1$$

$$y_2=k_1\left(1-\sqrt{2} \right)^{4}+k_2\left(1+\sqrt{2} \right)^{4}=5$$

This gives rise to the system:

$$3\left(k_1+k_2 \right)+2\sqrt{2}\left(k_2-k_1 \right)=1$$

$$17\left(k_1+k_2 \right)+12\sqrt{2}\left(k_2-k_1 \right)=5$$

Multiplying the first equation by $-6$, then adding, we get:

$$-k_1-k_2=-1\,\therefore\,k_1=1-k_2$$

Substituting into the first equation for $k_1$, there results:

$$3+2\sqrt{2}\left(2k_2-1 \right)=1$$

Solving for $k_2$, we find:

$$k_2=\frac{\sqrt{2}-1}{2\sqrt{2}}$$

Thus:

$$k_1=1-\frac{\sqrt{2}-1}{2\sqrt{2}}=\frac{\sqrt{2}+1}{2\sqrt{2}}$$

and so we find:

$$y_n=\frac{\sqrt{2}+1}{2\sqrt{2}}\left(1-\sqrt{2} \right)^{2n}+\frac{\sqrt{2}-1}{2\sqrt{2}}\left(1+\sqrt{2} \right)^{2n}$$

$$y_n=\frac{1}{2\sqrt{2}}\left(\left(\sqrt{2}+1 \right)\left(1-\sqrt{2} \right)^{2n}+\left(\sqrt{2}-1 \right)\left(1+\sqrt{2} \right)^{2n} \right)$$

So, in summary, we find an infinite number of solution pairs, given:

i) recursively:

$$x_{n+1}=6x_{n}-x_{n-1}$$ where $$x_1=1,\,x_2=7$$

$$y_{n+1}=6y_{n}-y_{n-1}$$ where $$y_1=1,\,x_2=5$$

ii) closed form:

$$x_n=\frac{1}{2}\left(\left(\sqrt{2}-1 \right)\left(1+\sqrt{2} \right)^{2n}-\left(\sqrt{2}+1 \right)\left(1-\sqrt{2} \right)^{2n} \right)$$

$$y_n=\frac{1}{2\sqrt{2}}\left(\left(\sqrt{2}+1 \right)\left(1-\sqrt{2} \right)^{2n}+\left(\sqrt{2}-1 \right)\left(1+\sqrt{2} \right)^{2n} \right)$$
 
Suppose ,instead of the usual x,y coordinate system with an I basis vector along the x -axis and a corresponding j basis vector along the y-axis we instead have a different pair of basis vectors ,call them e and f along their respective axes. I have seen that this is an important subject in maths My question is what physical applications does such a model apply to? I am asking here because I have devoted quite a lot of time in the past to understanding convectors and the dual...
Insights auto threads is broken atm, so I'm manually creating these for new Insight articles. In Dirac’s Principles of Quantum Mechanics published in 1930 he introduced a “convenient notation” he referred to as a “delta function” which he treated as a continuum analog to the discrete Kronecker delta. The Kronecker delta is simply the indexed components of the identity operator in matrix algebra Source: https://www.physicsforums.com/insights/what-exactly-is-diracs-delta-function/ by...
Back
Top