[Search for users] [Overall Top Noters] [List of all Conferences] [Download this site]

Conference rusure::math

Title:Mathematics at DEC
Moderator:RUSURE::EDP
Created:Mon Feb 03 1986
Last Modified:Fri Jun 06 1997
Last Successful Update:Fri Jun 06 1997
Number of topics:2083
Total number of notes:14613

1889.0. "How close to sqrt(7) can a rational number be?" by EVTSG8::ESANU () Mon Aug 22 1994 08:11

The following one is due to the Romanian mathematician
Radu Gologan (born 1952):

	Let m,n be integers.

	    m
	If  - < sqrt(7)
	    n
	           m    1
	then also  - + -- < sqrt(7)
	           n   mn

Generalizations?
T.RTitleUserPersonal
Name
DateLines
1889.1An infinite sequence approximation to sqrt(7)?VMSDEV::HALLYBFish have no concept of fireMon Aug 22 1994 10:4123
    (1) Choose m1,n1 (both > 0) so that m1/n1 < sqrt(7) and gcd(m1,n1) = 1
    
    Define m2,n2 by: m2 = m1�+1, n2 = m1n1, reduced to lowest terms, viz:
    
    	m2     m1     1     m1�     1     m1�+1
        --  =  -- + ---- = ---- + ---- = ------ then reduced to lowest terms.
        n2     n1   m1n1   m1n1   m1n1    m1n1
    
    Stipulating .0, we see this generalizes immediately to a sequence
    of (mN,nN) each of which satisfies (mN/nN) < sqrt(7).
    
    Question: does (mN,nN) --> sqrt(7) as N --> infinity, for *all*
    (m,n) satisfying (1)? If not, what criteria determines convergence
    to what value?
    
    In the most basic case, we have:
    
    (1,1) -> (2,1) -> (5,2) -> (26,10) => (13,5) -> (170,65) => (34,13)
     1.0      2.0      2.5       2.6 reduce 2.6     2.615+ reduce 2.615+
    
    -> (1156,442) => (34, 13)  A loop! Does this always happen?
    
      John
1889.2re .1: No loop! and a partial answerEVTSG8::ESANUTue Aug 23 1994 03:5326
First remark: the sequence defined in .1 cannot loop, as it is strictly
increasing: (m/n) + (1/(m*n)) > (m/n). The last pair in the example given in
.1 is (1157,442), not (1156,442).


Let us call "seed" the first term of the sequence of rationals defined in .1.


To the first question stated in .1:

		Does the sequence defined in .1 converge
		to sqrt(7) for all seeds?

I answer "no", as the proposition enunciated in .0 is equally true for sqrt(3)
(!), for instance, instead of sqrt(7): so build such a sequence for sqrt(3);
the built sequence is admissible for sqrt(7) also, as sqrt(3) < sqrt(7), but
it cannot converge to sqrt(7), as it never passes sqrt(3).


I have no answer for the second question in .1:
		
		What criteria determine convergence of the
		sequence defined in .1 to what value?

Does anyone have an idea?
Mihai.
1889.3SSAG::LARYLaughter &amp; hope &amp; a sock in the eyeTue Aug 23 1994 06:0619
>I answer "no", as the proposition enunciated in .0 is equally true for sqrt(3)
>(!), for instance, instead of sqrt(7)...

I don't think so; 1/1 < sqrt(3) but 1/1 + 1/(1*1) isn't.

There's a famous 19th century theorem by Liouville that limits how "good" a
rational approximation to an algebraic number can be; it says that for any
number x which is the root of an n-th degree irreducible equation and any
rational approximation p/q to x, 

	abs(x - p/q) > c/q^n    for some c > 0 which depends only on x

for square roots, n = 2, so from the structure of this you could certainly
surmise that there is some maximal fixed d, 0 < d < 1, such that 

	m/n < sqrt(3) ==> m/n + d/(mn) < sqrt(3)

But I haven't the foggiest idea of what d is, just as I haven't the foggiest
idea of how to prove the result in .0 from Liouville's theorem...
1889.4re. .3: You're right, but!...EVTSG8::ESANUTue Aug 23 1994 07:148
You're right, I gave a slightly wrong counterexample, though a type .0
proposition is true for sqrt(3) if you restrict it to  m > 1  and  n > 1 .

Consider also the .0 proposition for 2 instead of sqrt(3) or sqrt(7).
It works for all integer m and n, so this *is* a good reason for not all
seeds breeding sequences convergent to sqrt(7).

Mihai.
1889.5VMSDEV::HALLYBFish have no concept of fireTue Aug 23 1994 12:0418
    Thanks for pointing out the looping error. Serves me right for thinking
    while typing.
    
    (1,1) is readily seen to converge to 1+�, where � = (1+sqrt(5))/2,
    the ratio of even (or odd) Fibonacci numbers.
    
(1,1) -> (2,1) -> (5,2) -> (13,5) -> (34,13) -> (89,34) -> (233,89) ->
(610,233) -> (1597,610) -> (4181,1597) -> (10946,4181) -> (28657,10946)
    
    Since (1+�)� < 7 one can select integers m,n so that 1+� < m/n < sqrt(7)
    and since the sequence is monotonic increasing it is clear 1+� is not
    what EVERY possible m/n converges to.
    
    I conjecture that any given m,n run thru the process indefinitely will
    converge to (a + sqrt(b)/c), where a,b,c are integer functions of m,n.
    All we need is somebody with MAPLE to run a few different m,n ...
    
      John
1889.6HANNAH::OSMANsee HANNAH::IGLOO$:[OSMAN]ERIC.VT240Wed Aug 24 1994 12:436

The conjecture about m/n<sqrt(7) implying that adding 1/mn is still < sqrt(7),
is that a "famous unproved conjecture" or is there a published proof ?

/Eric
1889.7re. .5 & re. .3 - too long, I'm afraid...EVTSG8::ESANUWed Aug 24 1994 13:12177
re. .5 : This is a nice remark you made about the seed (1,1) breeding the
classical Fibonacci sequence! Could it be because of a coincidence? (see 6.2
hereafter).

Let me first recall some basic facts and detail your remark:


6.0. Let x(k) = m(k)/n(k) be the sequence of rationals built as you described it
in .1 . Then

	m(k+1)  = m(k)**2 + 1
	n(k+1) = m(k) * n(k)

(** means exponentiation). 


6.1. Let u(k) be the classical Fibonacci sequence:

	u(1) = 1
	u(2) = 1
	u(k+1) = u(k) + u(k-1)

u(k) is readily extended to all integer indices by

	u(0) =0
	u(-k) = u(n) * (-1)**(n-1)


6.2. It can be proved by induction that the seed (1,1) breeds the sequence

	x(k) = u(2k-1)/u(2k-3)

(use the identity   u(n)**2 - u(n-2)*u(n+2)  = (-1)**n ).

(:-  In fact, you get

	x(k+1) = (u(2k-3) * u(2k+1)) / (u(2k-1) * u(2k-3)) =
		= u(2k+1) / u(2k-1)
-:)

In particular   x(1) = u(1)/u(-1) = 1/1 = 1 , but the seed is
(m, n) = (u(1), u(2)) =(1,1) , the coincidence being that   u(2) = 1 = u(-1).

As, for a constant p,   u(n)/u(n-p) --> �**p , it results that
x(k) --> �**2 = 1 + � .


6.3. The set of sequences v(n) , recurrently defined (like Fibonacci's):

        	 p-1
	v(n+p) = Sum (a(i) * v(i))
	         i=0

is a vector space. If all the roots   {x(k} / 0 <= k <= p-1}  of the p-degree
polynomial

	       p-1
	P(x) = Sum (a(i) * x**i)
	       i=0

are distinct, then the sequences

	v(k,n) = x(k)**n    (0 <= k <= p-1)

form a vector basis for this vector space.


6.4. I like your conjecture about an arbitrary seed (m,n) breeding a sequence
convergent to some << (a + sqrt(b)/c), where a,b,c are integer functions of m,n
>> . Could it be that for an arbitrary seed (m,n) the sequence
x(k) = m(k) / n(k)   depends always on recurrent sequences like above (6.3)?

Feeling against: the sequences x(n) converge so rapidly, that some of them could
very well have transcendental limits.

I smell continued fractions...


6.5. The problem of the reducing to lower terms:

If no reducing to lower terms of m(k) and n(k) was made while building the
sequence x(k) = m(k) / n(k), then a product formula can be quickly calculated
for x(k), and the infinite product obtained is (rapidly) convergent. Also in
this case the initial n (from the seed (m, n) ) is only a scaling factor,
otherwise the formula for x(k) depends on m only (like in the Fibonacci case,
obtained for the seed (1, 1), where n = u(1) = 1 seemed to be there by chance).

But there necessarily *is* such a reducing to lower terms, at least concerning
the factor 2:

Let Zq be the ring of the modulo q integers.

Consider the mapping Fq: Zq x Zq -> Zq x Zq , defined by

	Fq(m, n) = (m**2 + 1, m * n)

Then	Fq(1, 1) = (0, 1)
	Fq(0, 1) = (1, 0)
	Fq(1, 0) = (0, 0)
	Fq(0, 0) = (1, 0)

and it loops, so while building the sequence x(k) = m(k) / n(k), there are
reducing to lower terms for the factor 2 every two terms of the sequence x(k),
starting at latest for k = 4.

Having a look about reducing to lower terms for the factors 3 and 4, I obtained
the inverse conclusion: in Im(F3) and in Im(F4), the pair (0, 0) is a dead
start, that is you never get to it unless you start with it. So, while building
the sequence x(k), there is never any reduction to lower terms for the factors 3
and 4, unless there is such a reduction in the passage from x(1) to x(2).

I couldn't get any farther, perhaps an algebraist could be of help?


===============================================================================
re. .3 : Let me recall some other Liouville-related results:
     
6.6. The real number x is << approximable by rationals of order n >> if there is
a constant K(x), depending only on x, such that

	abs(p/q - x) < K(x)/q**n

THEOREM. Any rational is approximable to order 1, and to no higher order.

THEOREM. Any irrational is approximable to order 2.

More precisely:
THEOREM (Hurwitz). In this case (irrational x) the constant K(x) can be taken as
1/sqrt(5) for all x.

THEOREM (Hurwitz). In this case (irrational x) the constant K(x) = 1/sqrt(5) is
the best possible one: there is no K(x) < 1/sqrt(5) which fits all irrational x .

THEOREM. A quadratic irrational (= quadratic surd) is approximable to order 2
and no higher order.

The theorem you mentioned:
THEOREM (Liouville). A real algebraic number of degree n (root of an integer
polynomial of degree n) is not approximable to any order greater than n.

COROLLARY. If a real number is approximable to any integer order, then it is
transcendent.


6.7. The case of the problem .0 (for sqrt(7)).

The constant c that you mentioned:

>	abs(x - p/q) > c/q^n    for some c > 0 which depends only on x

appears in the proof of Liouville's theorem as  c = 1/M, where

	M =     max   abs(f'(y))
            abs(y-x) <= 1

and f is the polynomial of which x is a root.

In the case of x = sqrt(7) we have:

	f(y) = y**2 - 7
	M = 2 * (sqrt(7) + 1)

We need   c/(q**2) > 1/(p*q)   in order to prove .0, i.e.   p/q > 1/c = M ,
which is false in this case.

Anyway, the Liouville inequality is valid only for p/q lying in an interval
(x-d, x+d), so... 


References:

G.H. Hardy, E.M. Wright, An Introduction to the Theory of Numbers, Oxford
University Press  << The Holy Writ ! >> 


We have strayed - so pleasantly - far away from the problem .0 , or haven't we?

1889.8re. .6: Disappointed?EVTSG8::ESANUWed Aug 24 1994 13:2711
Does everybody wishes to know?...

Well, the answer is that the problem 1889.0 was published by Radu Gologan at
the beginning of the '80s (or late '70s) in the Romanian mathematical magazine
for high-school students, "Gazeta Matematica (seria B)". Knowing that they
would publish only elementary problems could help you to solve it in an
elementary way - I can produce my elementary solution any time. But I don't
know how far it can be generalized, and then I am puzzled by the limits of the
sequences defined in .1 by John, and do Liouville-type results help here?

Mihai.
1889.9Typing errors in .7EVTSG8::ESANUThu Aug 25 1994 05:4023
For those who have been patient and merciful enough to read 1889.7, there
were some typing errors:

------------------------------------------------------------------------
**Correction of 6.3:**
Replace                                /             By
------------------------------------------------------------------------
                 p-1                   /            p-1                 
        v(n+p) = Sum (a(i) * v(i))     /   v(n+p) = Sum (a(i) * v(n+i)) 
                 i=0                   /            i=0                 
------------------------------------------------------------------------
------------------------------------------------------------------------
**Correction of 6.5:**
Replace                                /             By
------------------------------------------------------------------------
Then    Fq(1, 1) = (0, 1)              /   Then    F2(1, 1) = (0, 1)   
        Fq(0, 1) = (1, 0)              /           F2(0, 1) = (1, 0)   
        Fq(1, 0) = (0, 0)              /           F2(1, 0) = (0, 0)   
        Fq(0, 0) = (1, 0)              /           F2(0, 0) = (1, 0)   
------------------------------------------------------------------------

Sorry!
Mihai.
1889.10HANNAH::OSMANsee HANNAH::IGLOO$:[OSMAN]ERIC.VT240Mon Aug 29 1994 15:124
I'm still enough confused about .7 to not know whether it answer my question
asking if the sqrt(7)  thing is a known provable thing.

1889.11Yes, 1889.0 is provableEVTSG8::ESANUAu temps pour moiTue Aug 30 1994 04:5213
I would very much like to see non-elementary proofs of 1889.0. I have mine,
and it is elementary, I have even a slight generalization for it (regarding
other numbers than sqrt(7)), but I realize that I do not know all the numbers
which can replace sqrt(7) in 1889.0.

Also, consider the sequence defined in 1889.1: it is convergent, but what is
its limit? It is not clear to me that it is sqrt(7)...

This problem (1889.0) seems to be rather deep, it seems to be close to a lot
of fundamental ideas.

Mihai.

1889.12HANNAH::OSMANsee HANNAH::IGLOO$:[OSMAN]ERIC.VT240Wed Aug 31 1994 11:5167
let v7 mean sqrt(7) to save me typing.

v7 = 2 + 1/x.

x = 1/(v7-2).  Multiplying top and bottom by v7+2:

x = (v7+2)/3.

Since v7 is between 2 and 3, v7+2 is therefore between 4 and 5.  So x is
between 4/3 and 5/3, so:

x = 1 + 1/y.

y = 1/(x-1) = 1/((v7+2)/3 - 1) = 3/(v7-1).   [ I multiplied top and bot by 3 ].

Multiplying top and bottom by v7+1:

y = 3(v7+1)/6 = (v7+1)/2.

Since v7 is between 2 and 3, v7+1 is between 3 and 4.  So y is between 3/2
and 4/2, so:

y = 1 + 1/z.

z = 1/(y-1) = 1/((v7+1)/2)-1) = 2/(v7+1-2) = 2/(v7-1) = 2(v7+1)/6 = (v7+1)/3.

Gee, at end of the alphabet, let's wrap to a...

z is between (2+1)/3 and (3+1)/3 or 1 and 4/3.

z = 1 + 1/a.

a = 1/(z-1) = 1/((v7+1/3)-1) = 3/(v7-2) = 3(v7+2)/3 = v7+2.

Taking this all together,

v7 = 2+	1
	------------
	1 +	1
		-------------------
		1 +	1
			---------------
			1 +	1
				----------------
				v7 + 2

This means that the continued fraction denominators are this repeating pattern:

	v7 = [2 1 1 1 4 1 1 1 4 1 1 1 4 1 1 1 ....]

which is sort of interesting because of several things.  First, that what
looked "irrational" becomes just as regular as a repeating decimal when viewed
in this form.  Even transcendental numbers, such as e, confess to regularity.
e = [2 1 2 1 1 4 1 1 6 1 1 8 1 1 10 1 1 12 1 ....] although others, such
as pi don't.  (Does that make pi more transcendental than e ??)

Anyway, the other interesting thing about this continued fraction for 7 is
all the 1's.  Or, if not the 1's themselves, perhaps the fact that the
length before repetition is 4 denominators.  (Is length of repetition
related to length of repetition in decimal representation of 1/7 ?)

Does this representation of v7 help with the original puzzle ?  I thought
maybe it would if we could message the continued fraction into a form where
1/(AxB) is continually added on, but I don't see it yet...

/Eric
1889.13re. .12: pi is not more transcendental than eEVTSG8::ESANUAu temps pour moiTue Sep 06 1994 12:46128
                <<< RUSURE::NOTES1:[NOTES$LIBRARY]MATH.NOTE;7 >>>
                            -< Mathematics at DEC >-
================================================================================
Note 1889.13     How close to sqrt(7) can a rational number be?         13 of 13
EVTSG8::ESANU "Au temps pour moi"                   118 lines   6-SEP-1994 08:23
               -< re. .12: pi is not more transcendental than e >-
--------------------------------------------------------------------------------
Sorry for answering so late. I have also been working on your
interesting questions, Eric:

> This means that the continued fraction denominators are this repeating pattern:
>
>	v7 = [2 1 1 1 4 1 1 1 4 1 1 1 4 1 1 1 ....]
>
> which is sort of interesting because of several things.  First, that what
> looked "irrational" becomes just as regular as a repeating decimal when viewed
> in this form.  Even transcendental numbers, such as e, confess to regularity.
> e = [2 1 2 1 1 4 1 1 6 1 1 8 1 1 10 1 1 12 1 ....] although others, such
<<<< Q1 and Q2
> as pi don't.  (Does that make pi more transcendental than e ??)
> 
> Anyway, the other interesting thing about this continued fraction for 7 is
> all the 1's.  Or, if not the 1's themselves, perhaps the fact that the
<<<< Q3
> length before repetition is 4 denominators.  (Is length of repetition
> related to length of repetition in decimal representation of 1/7 ?)

<<<< Q4
> Does this representation of v7 help with the original puzzle ? 

(here v7 = sqrt(7))

-------------------------------------------------------------------
Q1. Even pi has some very regular expressions, as soon as we accept
non-regular continued fractions:

	     4  1**2 2**2 3**2 4**2
	pi = -  ---- ---- ---- ----
	     1+  3  + 5  + 7  + 9  +...


-------------------------------------------------------------------
Q2. So is pi 'more' transcendental than e?

The set of the numbers which have a periodic decimal representation
(starting after a finite number of decimals) is enumerable (it is Q,
the set of the rational numbers).

The set of the numbers which are the limit of a periodic regular
continued fraction (starting after a finite number of partial
quotients) is also enumerable.

Why is it so? Because it is the union of an enumerable set of (disjoint)
enumerable sets:
                            -------------
Let x = [a(1),a(2),...,a(m),b(1),...,b(n)] be the continued fraction
convergent to x and the barred part represent its periodic sequence of
partial quotients.

Then our set is =  U  {all sequences of digits of length m} x
                m,n in N
		     x {all sequences of digits of length n} =

                = U  (S(9)^m x S(9)^n)
                m,n in N

where S(9) = {0,1,2,3,4,5,6,7,8,9}

which is an enumerable set of finite sets, so it is enumerable.

Another proof of the enumarability of this set: it is also known since
Lagrange that this set here above (the set of the numbers which are
the limit of a periodic regular continued fraction, starting after a
finite number of partial quotients) is the set of the quadratic surds,
i.e. the numbers of the form

	A + sqrt(B)
        -----------
            C

(A,B,C integers,  B >= 0,  C <> 0).


Conclusion: this set being enumerable and dense in R, its
complementary has the power of the continuum and is dense in R, same
situation as Q and R\Q or algebraic numbers and transcendental
numbers. So, at a first glance, this confers to pi a kind of 'more
transcendental' flavour than e has, but even pi is part of a 'regular'
enumerable set (see 1 above)! So, finally, *in this sense*, really
every number is 'more transcendental' than any other number...


-------------------------------------------------------------------
Q3. Link between the periodicity of the decimal representation of 1/7
and the periodicity of the representation of sqrt(7) as regular
continued fraction.

I haven't succeeded in finding one... except the following loose one:

- the periodicity of the decimal representation of 1/q (q positive
integer) <= q

- the periodicity of the regular continued fraction convergent to
sqrt(q) <= card {(a,b) / a,b integers, a**2 + b**2 <= q**2}

(the number of points of integer coordinates in the origin-centered
disk of radius q).

Notation used: '<=' for 'less or equal'

Does anyone know better?


-------------------------------------------------------------------
Q4. Well, these ideas might help (I really don't know!). The original
puzzle does also have an *elementary* solution.

Mihai.


-------------------------------------------------------------------
REFERENCES:

* G.H. Hardy, E.M. Wright, An introduction to the Theory of Numbers,
Oxford University Press        << The Holy Writ! >>

* Claude Brezinski, History of Continued Fractions and Pade'
Approximants, Springer-Verlag, 12 SCM
1889.14HANNAH::OSMANsee HANNAH::IGLOO$:[OSMAN]ERIC.VT240Tue Oct 04 1994 15:5221

Wording the base note a bit differently:

	For all m/n (integer m,n) less than sqrt(7), it is always more
	than 1/(mn) less.

Question:

	What's the best m/n we can find ?  That is, find m,n (integers)
	that are best fits to this equation:

		m/n + 1/(mn) = sqrt(7).

As I ask this, I'm not sure if a qualifier is needed.  For example, perhaps
it's relatively easy to find m/n *larger* than sqrt(7) but by less than
1/(mn) ?  If not, then no qualifier is needed.  If so, then we need to say:

	FOR M/N LESS THAN SQRT(7) and m,n (integers), what m,n best fits...

/Eric
1889.15Rewording the rewording? Rewarding!EVTSG8::ESANUAu temps pour moiMon Oct 10 1994 12:1810
Eric, could you please reword your rewording? Pay attention to:

>	What's the best m/n we can find ?  That is, find m,n (integers)
>	that are best fits to this equation:
>
>		m/n + 1/(mn) = sqrt(7).

LHS is rational, while RHS is irrational...

Mihai.
1889.16AUSSIE::GARSONachtentachtig kacheltjesMon Oct 10 1994 23:1614
re .15
    
>LHS is rational, while RHS is irrational...
    
    I'm sure Eric knows that. That's why he asks for "best fits" and not
    solutions.
    
    I guess what he's getting at, to give an example, is that given
    m,n < 1000, any m/n will differ from sqrt(7) by at least 1E-6, so he
    wants to know how to find the m and n that get as close as possible to
    sqrt(7). Presumably brute force is not the answer.
    
    I thought there was a sizable body of theory around nearest rational
    approximants (but I don't know any of it).
1889.17Rapidly convergent rational sequences have transcendental limitsEVTSG8::ESANUAu temps pour moiThu Oct 13 1994 07:5547
re. .14, re. .16:


The best rational approximants to any real number are the partial quotients
of its expansion as a continued fraction. Clear and complete presentations
of the subject are made in the following classical books:

[1] G.H. Hardy, E.M. Wright, An Introduction  to the Theory of Numbers,
Oxford University Press

[2] Ivan Niven, herbert S. Zuckerman, Hugh L. Montgomery, An Introduction 
to the Theory of Numbers, John Wiley & Sons, inc.



Beware: The sequence m(k+1)/n(k+1), defined by

	m(0)/n(0) < sqrt(7)

	m(1)/n(1) = m(0)/n(0) + 1/(m(0)*n(0)) < sqrt(7)

	...

	m(k+1)/n(k+1) = m(k)/n(k) + 1/(m(k)*n(k)) < sqrt(7)

might well *not* converge to sqrt(7).

I even suspect some of these sequences to have a transcendental limit, as
some of them converge extremely rapidly. my feeling is based on the
following grounds (see the Liouville-type results in 6.6, note 1889.7):

iv) The slowest convergent sequences of rationals have rational limits.

iii) On the third place are the rational approximants to quadratic surds.

ii) On the second place are, generally, rational sequences convergent to
algebraic numbers.

i) The quickest convergent rational sequences have transcendental limits.



Perhaps in the best position to solve the problem .0 are 15 years old kids
who have not yet heard about continued fractions...

Mihai.

1889.18HANNAH::OSMANsee HANNAH::IGLOO$:[OSMAN]ERIC.VT240Fri Oct 21 1994 12:1644
As I calculated in .12, the continued fraction for v7 (square root of 7)
is

v7 = [2 1 1 1 4 1 1 1 4 1 1 1 4 1 1 1 ....]

Suppose we lop off everything after the first 4.  That approximation gives:

v7 ~ 2+1/(1+1/(1+1/(1+1/4))) = 2 + 9/14 = 37/14

As reminded by Mihai in .17, 37/14 is a "best" approximation, which means
that no fraction with the same or smaller denominator can be any closer
than v7.

	(by the way, my DCL says 37/14 * 37/14 = 6.98469)

However, it's not clear to me that this approximation is best when measured
according to see what happens when 1/mn is added on.

What I'm worried about is, suppose C/D is closer to v7 than a/b.  But if C and
D are larger integers than a/b, then 1/CD is much smaller than 1/ab.  Hence
isn't it possible that

	a/b + 1/ab is closer to v7 than C/D + 1/CD, even though

	C/D is closer to v7 than a/b ?


If so, then the continued fraction might not give us the best approximation
I'm looking for.

Perhaps another way to ask what I'm looking for:

	What is the limit of a/b + 1/ab for integers a,b as a/b approaches
	v7 from below ?

If you tell me the limit is v7, then give some support to your claim by
finding a,b such that

	a/b + 1/ab + c = v7

where a,b are integers and c is a positive real number less than
1E-20.  (reciprocal of 10 to the 20th power)

/Eric
1889.19AUSSIE::GARSONachtentachtig kacheltjesSun Oct 23 1994 03:287
re .18
    
>However, it's not clear to me that this approximation is best when measured
>according to see what happens when 1/mn is added on.
    
    Yes. I overlooked that when I posted .16. I don't know the answer to
    your question.
1889.20A generalization of .0EVTSG8::ESANUAu temps pour moiTue Nov 08 1994 03:3415
Let  k,m,n  be positive integers. The following
proposition is true for:

(i)    p = 4 * k  and  k, m, n > 0
(ii)   p = 4 * k + 3  and  k, m, n > 0
(iii)  p = 4 * k + 3  and  k = 0;  m, n > 1

	    m
	If  - < sqrt(p)
	    n
	           m     1
	then also  - + ----- < sqrt(p)
	           n   m * n

Mihai.
1889.21Proof of .20 (hence of .0)BATVX1::ESANUAu temps pour moiFri Dec 16 1994 09:1049
        m     1
(1)	- + ----- < sqrt(p)
	n   m * n

is equivalent to

                 1
(2)	m� + 2 + -- < p * n�
                 m�

For  m = 1  this means

        2
(3)	- < sqrt(p)
	n

which is true in the cases (i) and (ii), the case (iii) supposing m <> 1:

(i)    p = 4 * k  and  k, m, n > 0
(ii)   p = 4 * k + 3  and  k, m, n > 0
(iii)  p = 4 * k + 3  and  k = 0;  m, n > 1

(The case (iii) was introduced exactly for this purpose).

For  m > 1  we have  0 < 1/m� < 1 , so, as by our assumption  m� < p * n� , 
for proving (2) we only need to show that the following cases are
impossible:

(3)	m� + 1 = p * n�
(4)	m� + 2 = p * n�

which is straghtforward in the cases (i), (ii), (iii), by divisibility
(odd/even) considerations.

--- --- ---

I am not satisfied with this proof, though it has the merit of being
elementary: it does not penetrate the essence of an assertion of the type:

            m                 m     1
(5)	If  - < x  then also  - + ----- < x
            n                 n   m * n


For which x is this true, and for which x is this false?

It is clear to me that (5) is very deep, as the ideas stirred by .0 showed.

Mihai.
1889.22RTL::GILBERTMon Jan 02 1995 14:1010
Let S be the set of x > 0 such that:

	m               m     1
	- < x  implies  - + ----- < x, for any m and n positive integers.
	n               n   m * n

The numbers 1 and (3+sqrt(5))/2 are the two smallest elements of S.
The proof is long, tedious, and left to the reader.

(3+sqrt(5))/2 = 2.6180...
1889.23HANNAH::OSMANsee HANNAH::IGLOO$:[OSMAN]ERIC.VT240Tue Jan 03 1995 10:3529
Hmmm.  If (3+v5)/2 is the smallest S, then (1+v5)/2 has a counterexample.

In other words, we can find

	m/n < (1+v5)/2

such that

	m/n + 1/(mn) >= (1+v5)/2

First, let

	g = (1+v5)/2

so my hand doesn't get tired keeping typing the expr.  (g for golden ratio?)

Well, v5 is about 2.2 something, so

	g = 3.2/2

Try m/n = 3/2.

Question is, is 3/2 + 1/6 larger (or equal) than (to) g ?  I don't think so,
since 1/6 is .1666... which is smaller than .2...

So this m/n isn't a counterexample for g.  Can any of you find one ?

/Eric
1889.24you found one!IOSG::TEFNUT::carlinDick Carlin IOSG ReadingTue Jan 03 1995 11:109
>Question is, is 3/2 + 1/6 larger (or equal) than (to) g ?  I don't think so,
>since 1/6 is .1666... which is smaller than .2...

Question is, is 3/2 + 1/6 larger (or equal) than (to) g ?  I think so,
since 1/6 is .1666... which is greater than .2/2...

Dick

1889.25HANNAH::OSMANsee HANNAH::IGLOO$:[OSMAN]ERIC.VT240Tue Jan 03 1995 14:334
whoops, yeah, sorry, i forgot the /2 ...

/Eric
1889.26re .22: Outline of the proof?EVTSG8::ESANUAu temps pour moiWed Jun 21 1995 07:4014
> Let S be the set of x > 0 such that:
>
>	m               m     1
>	- < x  implies  - + ----- < x, for any m and n positive integers.
>	n               n   m * n
>
> The numbers 1 and (3+sqrt(5))/2 are the two smallest elements of S.
> The proof is long, tedious, and left to the reader.

Can you give us the outline of the proof? Also, how does the set S look
like? What are its connected components, for instance?

Thank you,
Mihai.
1889.27RTL::GILBERTWed Jun 21 1995 22:1343
> Can you give us the outline of the proof?

Let S be the set of x > 0 such that:

	m               m     1
	- < x  implies  - + ----- < x, for any m and n positive integers.
	n               n   m * n

The numbers 1 and (3+sqrt(5))/2 are the two smallest elements of S.


Suppose some 0 < x < 1 is in S.  Let m = 1, and let n be the smallest positive
integer such that 1/n < x.  Then m/n + 1/(m*n) = 2/n > x, so x is not in S.
To see that 2/n > x, note the following:

	if       x > 1/2, then n = 2, and 2/n = 1 > x;
	if 1/2 > x > 1/3, then n = 3, and 2/n = 2/3 > 1/2 > x;
	if 1/3 > x > 1/4, then n = 4, and 2/n = 2/4 > 1/3 > x;
	...
	if 1/(j-1) > x > 1/j for an integer j >= 3, then take m=1 and n=j,
	so that m/n + 1/(m*n) = 2/j > 1/(j-1) > x.

Claim: 1 is in S.  This is equivalent to showing that there are no integers
m and n such that m/n < 1 < m/n + 1/(m*n); proof is left to the reader.
Hint: 1-m/n > 1/n.

To show that (3+sqrt(5))/2 is the second-smallest element of S, it's
necessary to show that for any 1 < x < (3+sqrt(5))/2, it's possible to find
m and n such that m/n < x < m/n + 1/(m*n).  To show that (3+sqrt(5))/2 is
in S requires showing that there are no such m and n.  In both cases, the
fractions approximating (3+sqrt(5))/2 were taken from the partials of the
continued fraction expansion of (3+sqrt(5))/2.


> Also, what does the set S look like?  What are its connected components,
> for instance?

That's a good question.  S doesn't have any 'connected components' -- it's
just a set of isolated points.  Suppose otherwise, that all x in some range
(x0,x1) were in S.  Pick any rational m/n in this range.  Then S contains
no points in the range (m/n,m/n + 1/(m*n)) -- but we assumed S contained
(x0,x1), so we have a contradiction.

1889.28re. .27EVTSG8::ESANUAu temps pour moiThu Jun 22 1995 09:046
Thank you a lot. Can you find all the elements of S? (Must be tough).

Mihai.

P.S. For the record, an isolated point of a set (the set defined by the point)
*is* a connected component of that set.
1889.29.21 seems to have slight inaccuracyHANNAH::OSMANsee HANNAH::IGLOO$:[OSMAN]ERIC.VT240Tue Oct 10 1995 15:5723
The following is excerpted from .21:


>       2
>(3)	- < sqrt(p)
>	n
>
>which is true in the cases (i) and (ii), the case (iii) supposing m <> 1:
>
>(i)    p = 4 * k  and  k, m, n > 0



For n=1 and k=1, we have p=4, so the statement becomes

	2 < sqrt(4)

which of course isn't exactly right.

Probably doesn't change subsequent discussion, but this is math, right ? So
such inaccuracies are worth mentioning ?

/Eric
1889.30AUSSIE::GARSONachtentachtig kacheltjesTue Oct 10 1995 19:565
    re .29
    
    I didn't check .21 or .29 but, yes, a proof that has a strict inequality
    which should actually include equality as well can doubtless be used to
    prove all sorts of contradictions.