MathGroup Archive 1999

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: Re: a tricky limit

  • To: mathgroup at smc.vnet.net
  • Subject: [mg15415] Re: [mg15387] Re: a tricky limit
  • From: "Richard Finley" <rfinley at medicine.umsmed.edu>
  • Date: Wed, 13 Jan 1999 20:57:29 -0500
  • Sender: owner-wri-mathgroup at wolfram.com

I think there may be a more direct way to see the limit of this
function... .if you multiply the product terms out you will see that
they approach 2/n  as x->1  (or 2/(n+1) as Paul found which is the same
thing since n->Infinity) ....since the limit (if it exists) cannot
depend on how x approaches 1  from below, let x = (n-1)/n which gives
1/(1-x) = n .  Finally, then  the whole function limit is 2 n/n as n->
Infinity or limit ->2.  That is,  the limit is exactly 2 (although I
haven't strictly proven that the limit  exists).  This is why it
appeared to be near 2!  

regards,  R Finley


>>> Paul Abbott <paul at physics.uwa.edu.au> 01/12 2:14 AM >>>
Arnold Knopfmacher wrote:

> I wish to obtain a numerical estimate (say 8 decimal digits) of the
> limit as x  tends to 1 from below of  the function
> h[x]=(Product[(1-fm[x]/(m+1)),{m,2,Infinity}])/(1-x) where
> fm[x]=x^(m-m/d) and d is the smallest divisor of m that is greater 
than
> 1. The problem is that when I replace Infinity by say 1000 as the upper
> limit  of the product, the function blows up near 1. Visual inspection
> of the graph of h[x] for 0<x<0.9 say, suggests that the limit should
> have a value around 2.1. Can anyone help?

With the definitions,

In[1]:= h[n_][x_] := Product[1 - f[m][x]/(m + 1), {m, 2, n}]/(1 - x)
In[2]:= f[m_][x_] := x^(m - m/Divisors[m][[2]])

series expansion of truncated products

In[3]:= h[10][x] + O[x, 1]
Out[3]=
    -2        1153
----------- + ---- + O[-1 + x]
11 (-1 + x)   1155

In[4]:= h[20][x] + O[x, 1]
Out[4]=
    -2        40020454
----------- + -------- + O[-1 + x]
21 (-1 + x)   33948915

indicates the leading (problem) term and the low-order approximations to
the limit you are after.  You can easily show that the leading
(problem) term, i.e.,

In[5]:= Product[m/(m + 1), {m, 2, n}] Out[5]=
  2
-----
1 + n

vanishes as n tends to Infinity.  

Approximations to the limit you are after can be found by determining
the coefficient of x in the Product with x replaced by 1-x, e.g.,

In[6]:= Coefficient[Product[1 - f[m][1 - x]/(m + 1), {m, 2, 10}], x]
Out[6]=
1153
----
1155

Alternatively, simple series rearrangement gives

In[7]:= 2/(1 + 10) Sum[1 - 1/Divisors[n][[2]], {n, 2, 10}] Out[7]= 1153
----
1155

Hence the number you are after is the same as

	Lim[2/(1 + p) Sum[1 - 1/Divisors[n][[2]], {n, 2, p}],p->Infinity]

I am not aware of a way of computing this expression exactly.  However,
you can easily find numerical approximations:

In[8]:= Timing[N[2/(1+10000) Sum[1 - 1/Divisors[n][[2]], {n,2,10000}]]]
Out[8]= {25.9 Second, 1.33959}

In[9]:= Timing[N[2/(1+20000) Sum[1 - 1/Divisors[n][[2]], {n,2,20000}]]]
Out[9]= {65.85 Second, 1.33972}

Cheers,
	Paul

____________________________________________________________________
Paul Abbott                                   Phone: +61-8-9380-2734
Department of Physics                           Fax: +61-8-9380-1014
The University of Western Australia            Nedlands WA  6907      
mailto:paul at physics.uwa.edu.au  AUSTRALIA                      
http://www.physics.uwa.edu.au/~paul

            God IS a weakly left-handed dice player
____________________________________________________________________



  • Prev by Date: getting the name of a symbol
  • Next by Date: Re: bug in the "Calculus`Limit`"
  • Previous by thread: Re: a tricky limit
  • Next by thread: Re: a tricky limit